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1. Introduction 


To aid engineers that design power supply systems, special analytical tools 
are being introduced. These tools include integration routines that start with 
the description of a power supply in state-equation form and yield analytical 
results that include transient performance, steady-state equilibrium stability, 
and the response of the power supply to spurious inputs (audio susceptability) • 

The Department of Electrical Engineering at the University of South Florida 
has worked to develop two analysis tools that can be used with the state-equation 
analysis package. 

The first tool uses a computer program that works with the SUPER*SCEPTRE 
circuit analysis program and prints the state equations for an electrical net- 
work. The state equations developed automatically by the computer program are 
used to develop an algorithm for reducing the number of state variables required 
to describe an electrical network. In this way a second tool is obtained in 
which the order of the network is reduced and a simpler terminal model is obtained. 
The price of the reduction in order is a reduction in the accuracy with which the 
model represents the system. 

At the conclusion of the first year, the computer program that converts the 
circuit description of an electrical network to a set of state equations was 
completed. 

In the second year the computer program was modified to accept a larger 
class of circuits including circuits with mutual inductance. The work on re- 
ducing the number of state variables required to describe an electrical network 
was completed and is described in this report. 

The first task that was completed was the development of a computer pro- 
gram that would print the state-equations for an electrical network. The ana- 
lytical background for the state equations of an electrical network is presented 
in Chapter 2. From the analysis in Chapter 2, the generality of the approach 
can be seen together with the reason for the presentation of data as matrices 
which are inverted in the state equations. The analytical basis for reducing 
the number of state variables using eigenvalue sensitivities concludes 
Chapter 2. 
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In Chapter 3 a description of the steps required to get a set of state 
equations is given. Chapter 3 is a user’s manual for the program to PROd uce 
^tate Equations (PROSE program) . The user begins with an electrical network.. 

By following the directions in Chapter 3, he can obtain the state equations for 
the network. Additional features of the program are also described in 
Chapter 3. 

The description of the programs as written in FORTRAN is given, in Chapter' 

4 and a complete listing of the program appears in the Appendix. 

Several examples are given in Chapter 5 that illustrate the use of the PROSE 
program in obtaining the state equations for an electrical circuit. 

A summary of the work accomplished under this grant appears together with 
recommendations for continued research in Chapter 6. 
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2. State Equations 

The standard form of the state equations for an electrical or mechanical 
system is 

» 

X •= f(x,u,t) 
y = s(x,u,t) 

where x is an n x 1 state vector, 
y is an m X 1 output vector, 
u is a p X 1 input vector, 
t is time, 

f is the state transition function and 
g is output function* 

When the system is linear, the equations can be n-nritten as 
X ■ = Ax + Bu 
y = Cx + Du 

where A is n x n transition matrix 
B is n X p 
C is m X n 
D is m X p 

The matrices A, B, C, and D may be functions of time. In certain nonlinear 
problems the matrix representation is possible. The elements of the matrix 
are given as an algebraic expression in terms of circuit components such as 
resistors,, inductors and capacitors. The value of a component may depend upon 
a voltage or current (the element may be nonlinear) . 

The state equations for an electrical circuit can be derived by the 
following steps: 

The circuit diagram is the starting point. The model for each 
circuit device is obtain as a collection of two terminal elements. 

The final circuit diagram contains 

a. capacitors 

b. resistors 

c. Inductors 

d. independent voltage sources 

e. dependent voltage sources 

f. independent current sources 

g. dependent current sources 
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2. The graph of the circuit is used to obtain a tree with the tree 
branches selected in the following order: 

a. voltage sources 

b. capacitors 

c. resistors 

d. inductors 

All voltage sources are selected before any capacitors are selected; 
capacitors are selected to add new nodes to the tree before any resistors are 
added; etc. The circuit elements which are not tree branches ate called links 

3. The components are partitioned and listed in Table 1. The groups of 
components in the tree are listed across 'the top and the groups of 
components that are links are listed along the side. 

TABLE 1 



C4 

R5 

L6 

E7 

Cl . 

F14 

0 

0 

F17 

R2 

F24 

F25 

0 

F27 

L3 

F34 

■ F35 

F36 

F37 

J8 

F84 

F85 

F86 

F87 


where 

Cl is a diagonal matrix of capacitor links 
R2 - is a diagonal matrix of resistor links 
L3 is a matrix. of inductor links 
J8 is a column matrix of current sources and 

C4 . is a' diagonla matrix of capacitor tree branches 
R5 is a 'diagonla matrix of resistor tree branches 
L6 is a- matrix of inductor tree branches 
E7 is a column matrix of voltage sources and 
The symbols are matrices that shown the connections between the 1 
links and the j tree branches. 

-ff i' number of links in class p 

ij I j 1 5, Q .5. number of tree branches in class j 
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where 


fpq 


1 if tree branch q is in the loop formed by placing link p 
back in the tree and if q is in the same direction in the 
loop as p. 

“1 if tree branch q is in the opposite direction in the loop 
defined by p. 

0 if tree branch q is not in the loop defined by p. 


The state variables are selected as the column matrix 


13 

V4 • 


Inductor Current State Variables 


If mutual Inductances are allowed, the equations describing the inductor 
currents are given as 

L3 13 + L36 16 = V3 (1) 

L63 13 + L6 16 = V6 (2) 

where L3 is the inductance matrix for the inductor links 

L36 
L63 

L6 is the inductance matrix for tree-branch inductor. 

Using the connection equations 

16 = 13 + J8 


} 


are the inductance matrices for mutual inductance between 
link inductors and tree-branch inductors and 


'36 ■ "se 

which may be differentiated to give 


16 = F 


T 

36 


13 




Equation (4) may be substituted into Equation (1) to get 

• [L3 + L36 F-^ ] 13 = V3 - L36 F.^ J8 

36 86 


(3) 

(4) 

(5) 


The connection equations relating link voltages to tree branch voltages givea 
V3 = -F^^ V4 - Fg^ V5 - F^g V6 - F^^ E7 (6) 

The inductor vcltage matrix is removed by the following steps beginning with 
Equation (2) which is repeated next. 

V6 = L63 13 + L6 16 


= L63 13 + L6 (F^g 13 + F^,J J8) 


86 


( 2 ) 

(7) 


The equation for 13 becomes 


[L3 + L36 F^g + F36 L63 + F36 L6 F 


36 


] 13 


= -F35 V4 - F35 V5 - F37 E7 


( 8 ) 
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To reduce this equation so that only the state variables 13 and V4 to- 
gether with independent sources E7, Ey, J8, and Jx appear, the resistor volt- 
ages and currents must be obtained as a function of the state variables. 


Resistor Links and Tree B ranches 

The resistor voltages and currents can be expressed by the following 
matrix equations. 


R2 

0 


’12 


~V2 

0 

G5 


V5 


I5_ 


0 -E 
T 

LP35O 


24 


13 




0 -F 
T 

l’^25 “ 


25 


12 

V5 


0 

1^85 ° 


1 

'J8' 


E7 


Equation (9) can be solved as 


'l2~ 


R2 

-1 

1 

0 

1 

>ri 

1 


’ 13* 

4* 

0 

-E 

27 


-js' 

V5^ 




.^35 0 . 


V4 


F ^ 
L 85 

0 


E7 


If the matrix 


partitioned as 


then 






R2 


25' 


-^25 


-1 


xs 


ml 

m3 


m2 

m4 




and 


(9) 


( 10 ) 


*R2 F * 


mlR2 - m2F^_ ralF^c + w2G5 
25 25' 


U 0 

:^23 


m3R2 - m4F'?,. m3F.,_ + m4G5 

25 25 J 


0 U 

L J 


R2 

^25 

«-l 

R2ml 4- Foc™3 

25 

R2m2 + F25B14 


U' 

0 

T 

G5 


T 

T 




-^25 

-F23ml + G5m3 

-F23itt2 + G5m4 


0 

u 

- 

_ 


- 






( 11 ) 


( 12 ) 


where U is a unit matrix of the appropriate size. The solution for the matrix 
-2 

MR is obtained from Equation (11) as 
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f T 1 -1 

ml » 1 r 2 + 

m2 « -mlF25R5 
T 

m3 = R5F25ral 

m4 “ — E.5F 25 | 

Beginning With Equation (12) yields another expression for 
and as 

M - [^^5^2^25 + G5 ] 

m3 = m4F2^2 

m2 ® — G2F22^^ 

ml = G2 [u - G2F23m4F25] 


Solving Equation (10) for V5 yields 
V5 = m4P25l3 - ~ mSF^^E? + m^Fg^JS 

Equation (15) is substituted into Equation (8) and solved' to give 
the solution for state variable 13 as 


-® 35 **’’ 35^3 ■ -[ 


*34 ' * 35"®24 


13 -Mj_ 

-*35’'‘*85^8 ■ *36“*86f 


]V4- [ 


*37 - *35"**27 


] 


E7 


( 13 ) 


(14) 


(15) 


( 16 ) 


Canacitor Voltages 

The voltage state equation ca-n be obtained beginning xd-th the capacitor 
equations 


ClVl = II’ 
C4V4 == 14 


(17) 


The solution .begins with the second equation in (17) 

C4V4 = 14 = fJiI + F^Jl2 + F3^I3 + FgJj8 

The variables II and 12 must be eliminated in the solution. Begxnning wxth 
11 

Cl”^rl == VI = 


(19) 
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( 20 ) 


== + ¥^^12 + F^Jl3 + J8>] 

“1 —1 

Denoting Cl by SI and C4 by S4 and solving gives 

II . M/l {-F^^ S4 [fJi2 + + fJj81 


where M. = [SI + F,,S4F,h 
s 14 14 

Solving Equation (10) for 12 given 


( 21 ) 


12 = -ml[F24V4 + F2 qE 7] + m2 [F^J 13 +.EggJ8] (22) 

When Equations (21) and (22) are substituted into Equation (18) and the result 
simplified, the result is 

V4 - S4{[[-[fJ 341 + [FjJbZFj^ + 

* 'hK' h4^^^2l “IV'' 

+ Fj^^S4 - U 44 I IF 3 J + rj mZFg^l J3 

+ h4“ - VI 

Fiyi ?} C 23 ) 

Combining Equation (16) and Equation (23) gives the state equations. 


13 


1 

0 

1 L 

-1 

~^35”4F35 F^^mSF^^ 

“13" 

-V4- 


- Q C4_ 


JFaa-F, Im "V , S J [F Jm 2F,W. ]-[u, , -F, Vf, ,S, ]F JmlF,, ], 

V4. 


Lt « 44 -=' l 4 “ fh 4 S 4 "’' 8 l +^ 2 l '" 2 r 85 l 


“^37 ^35 ^ ^27 



' J8" 


- E7_ 


-“ 6^86 -^ 36«^36 


-F.J M ^ F._ 
14 s 17 



J8 ■ 


_ E7 . 


(24) 
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The state equations are produced as the output of the preprocessor computer 
program known as SCEPTRE I. These state equations are the input to the computer 
program SCEPTRE II which performs the integration required to obtain a transient 
solution. 

Notice that Equation (24) appears in matrix form even though each element 
may depend on time or on another voltage or current that can be obtained at 
times when the state variables are known. 

The state equations are ca] culated by processing the output of SCEPTRE I 
to obtain matrix equations in the form similar to Equation (24) . \'Jhile this 
equation is non-linear, it can be linearized arquM a quiescent value as specified 
by initial conditions for the state variables 

The linearized state equations are 




A 




+ E 



(26) 


The eigenvalues of this equation are eigenvalues of the constant matrix, A. 

In the representation of any physical system a compromise is made in 
choosing the order of the mathmatical model. The order may be order of the 
differential equation that describes the system or may be the number of poles 
in the transfer function for a linear system. In any physical system that 
is modeled with lumped parameters, a decision must be made on how many energy 
storage elements (capacitors, inductors) are included in the model. More 
energy storage elements increase accuracy and computational difficulty. The 
number of independent energy storage elements is equal to the number of state 
variables required . 

As the number of state variables is increased the accuracy can be in- 
creased. However, in choosing the equivalent circuit for a physical system 
such as a power supply, some energy storage elements may he selected that 
do not affect the accuracy of the model significantly. A simpler description 
that did not reduce accuracy would give a better model. 

In a state variable description of a physical system, when an independent 
energy storage element is discarded, the number of state variables is reduced 
by one. Energy storage elements that are incidental to the accurate description 
of a physical system usually contribute to only the shortest time constants 
or the largest eigenvalues of the A matrix. 
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The eigenvalues of the A. matrix are obtained using a subroutine. The 
sensitivity is calculated for each parameter by increasing the parameter by 
IZ and calculating the eigenvalue. For example the sensitivity to a resistor 
is given by 


•R P AR 


(100) ^ 

-p 


(27) 


where p is the eigenvalue, 

Ap is the change in the eigenvalue 
AR is the change in resistor and ■ 

R is. the resistor value. 

The sensitivities can be used to reduce the number of state variables. 
First the eigenvalue with the largest magnitude is selected. The sensitivities 
of this eigenvalue to each parameter is scanned. The circuit parameter 
(resistance, inductance, or capacitance) to which the largest eigenvalue has • 
a sensitivity with a magnitude within 3 per cent of unity is selected. If 
the selected circuit parameter does not affect other eigenvalues, it is . 
removed from the network. It ‘is set equal to zero if the sensitivity is 
negative; if the sensitivity is positive, the circuit parameter is set to 
infinity. 

The desire is to force the eigenvalue to negative infinity. From the 
equation- 


^ = y sP 

p R 


AR 

R 


-I- I sP ^ J ^ 

R ^'i “ C ^ L 


L L 


(28) 


it can be inferred that the per cent change in an eigenvalue can be maximized 
by selecting a circuit element, 5, with the largest sensitivity and maximizing 
oP' AC 




(29) 


for that term. 

In every case the eigenvalue will be moved away from the origin toward 
infinity. If it is moved to infinity, a state variable is eliminated by the 
removal of a circuit element. The correlation of eigenvalue going to infinity 
and state variable disappearing can be observed from the characteristic equation 

ISI - aI = + a + ... + a,S + a^ (30) 

' ‘ n-1 10 


In factored form the caharacteristic equation becomes 

K ... (-S + i). 

Pi P2 Pn 


(31) 
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Obviously one term disappears the largest eigenvalue p^. 

If not moved to infinity, the eigenvalue is made larger in magnitude 
and the process can be repeated so that other circuit elements are deleted 
and the eigenvalue is moved to infinity. 

By considering the characteristic equation again it can be seen that 

|SI - a| = s“ + a + ... + a- (32) 

' ' n-1 0 

where 



In words, the sum of the pole sensitivities to a parameter is less 
than or equal to 1, For this reason in most cases of interest if the sensitivity 
of an eigenvalue is nearly equal to one, the sensitivity of other eigenvalues 
will be very small. Thus, the circuit element is deleted will only cause one 
eigenvalue to have an appreciable change in magnitude. 

The sensitivity calculations are made in the PROSE program for the 
eigenvalue with the maximum absolute value. If the magnitude of the sensitivity 
is within three per cent of unity for a circuit parameter, the circuit parameter 
can be removed from the circuit. In the examples in Chapter 5 the computer 
output shows how the circuit elements can be removed. 
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3. User's Manual for PROSE 


3.1 Introduction 

The user's manual for the computer program to PROd uce ^tate Equations is 
given in this section. The PROSE program can be used to derive the state 
equations for any electrical network that can be analyzed with the SCEPTRE* 
circuit analysis program and compute eigenvalue sensitivities with the 
following restrictions: 

1. Each circuit component can be modeled by two-terminal elements 
which can be represented as either a voltage source, a current 
source, a resistor, a capacitor, or an inductor. 

2. No FORTRAN Subroutines are allowed. 

3. No linear, secondary, or primary dependent sources are permitted. 

(In place of the SCEPTRE dependent sources, voltage or current 
sources can be defined by expressions which can give a dependent 
source) . 

The PROSE program can he' used by following these steps: 

Step 1. A circuit diagram is drawn for the circuit for which the 
state equations are desired. Each component in the original circuit has 
been represented as a collection of two terminal elements in the circuit 
diagram. 

Step 2. The circuit diagram is translated into a- listing for tje 
SUPER*SCEPTRE circuit analysis program.* 

Step 3. Input sources are requested under OUTPUTS subheading. 

Step 4. After the RUN CONTROLS heading, the instruction EXECUTE 
SETUP PHASE ONLY is inserted. 

Step 5. The control cards are added and the complete program is entered. 

The PROSE program responds by outputting the state equations and a state 
transition matrix for a linear set of equations. In addition the PROSE program 
computes a set of "sensitivities” of the eigenvalues of the transition matrix 
•that result from a change in each of the circuit elements. 


* The format for the computer input can be found in one of the following: 

1. Extended SCEPTRE Volume I. User's Manual . AD-A009-594, 

2. James C. Bowers and Stephen R. Sedore, SCEPTRE: A Computer Program 
For Circuit Analysis and Systems Analysis , Prentice-Hall, Inc., 
Englewood, Cliffs, New Jersey 1971 

3. User's Manual for SUPER*SCEPTRE , AD-AOll-348 
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3.2 


Example 

To fix ideas consider as an example obtaining the state equations - 
for the circuit in Figure 1. 



Fig. 1. Boost Power Supply Average 
Model TD Interval Model 


The input description for PROSE use is: 

Line Number 


7 /job card for 240 k byte region 1 
//STEPl EXEC PROSE ^ 
//INPUT DD * 3 
CIRCUIT DESCRIPTION 4 
BOOST POWER SUPPLY 5 
AVERAGE MODEL 6 
TD INTERVAL 7 
ELEMENTS ,8 
EGjG-l=l 9 
RL, 1-2=1 10 

Lj2*"G=2 11 

C,3-G=3 12 
RC, 4-3=4 13 
R,4-G=5 14 
OUTPUTS 15 
EG 16 
RUN CONTROLS 17 
EXECUTE SETUP PHASE ONLY 18 
STOP TIME=5 19 
END . 20, 


Fig. 2. PROSE Input Listiiig for Circuit of Figure 1 
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The PROSE system output is condensed below. 

DL= (-IL*RL+EG) / (L) 

DC=-VC/((R+RC))/(C) 

STATE VARIABLE NAMES 
IL 
VC 

MATRIX A 

-S.OOOOOD-Ol 

0.0 

INDEPENDENT SOURCE NAMES 
EG 

MATRIX B 

-5.00000D-01 
-3.70370D-02 

EIGENVALUES OF MATRIX A 
REAL PART 
-5.0000000D-01 
-3.7037037D-02 

Fig. 3. Summarized PROSE output for input 
listing of figure 2. 

Figure 1 indicates the circuit whose state variable equations are 
computed. Each element is given a name and each node is labeled with 
an encircled character. 

Figure 2 shows the IBM OS job control language (lines 1 through 3) 
to execute the system. Lines 4 through 14 in Figure 2 describe the 
circuit in standard SCEPTRE forrmat.^ The OUTPUTS request (lines 15 & 16) 
alerts PROSE to regard EG as an independent source in its determination 
of the B matrix. 

The RUN CONTROLS (lines 17 through 19) include the special line: 


See EXTENDED SCEPTRE, Volume I. User's Manual, AD-A009-594 


IMAG. PART 
0.0 
0.0 


-5.00000D-01 

-3.7O370D-O2 
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EXECUTE SETUP PHASE ONLY. Both this card' and the STOP TIME ~ number 

dummy value are- required. The END card finishes the PROSE input description. 

The output of Figure 3 indicates the state variable equations for 
DL s^’^d DC as well as the state variable names (IL and VC) , the 

linear transition' matrix A, the independent sources (EG) , the input matrix 
B, and the eigenvalues for the transition matrix. 

The matrices are determined from the state variable equations and the 
element values in the form: X = AX + BU, where X is the vector of state 
variables and U is the vector of independent sources. 

Additional outputs include the sensitivities which will be further 
described later in this document. 

3.3.0 PROSE Language 

The PROSE system input language is the standard- SCEPTRE transient 
solution language. For those unfamiliar with the SCEPTRE language, the- 
reference of note 1 must be consulted. 

3.3.1 Descriptions 

PROSE will allow all standard SCEPTRE descriptions; i.e. CIRCUIT, 

MODEL, and RERUN. However,, the rerun description is ignored and model' 
description is used only if models under it are called in the circuit 
description. 

3.3.2 Elements 

All SCEPTRE elements (E,R,L,C,J) are allowed by PROSE. Element values 
may be constants, expressions or tables. PROSE will allow defined para- 
meters in expressions but will not allow FORTRAN function subroutines, 
equations, secondary’ dependent sources, linear dependent sources, and 
primary dependent current sources. Voltage and current source deriva- 
tives are also- allowed. 

Models may be called as elements as in standard' SCEPTRE, but model 
elements must also conform to the above rules. 

Since PROSE is transient analysis oriented, no AC sources nor ele- 
ments ^dth bounds are permitted. 
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3.3.3 Defined Parameters 

Standard SCEPTRE transient oriented real defined parameters may be utilized. 

3.3.4 Outputs 

At least one output request is necessary for proper execution. The 
OUTPUT requests of voltage or current sources implies chat these sources are 
independent’ sources and are used in the PROSE computation of the B matrix. 

3.3.5 Initial Conditions 

Initial conditions may be supplied and sometimes are required since 
PROSE does not compute a DC operation point prior to matrix computations. 

The initial conditions are used to set the system operating point when a 
linear set of state equations are computed. If the initial conditions are 
not supplied, all initial state variable values are assumed to be zero. 

3.3.6 Functions 

Only Tables are allowed under functions, PROSE does not permit equations. 

3.3.7 Run Controls 

PROSE requires the special run control: EXECUTE SETUP PHASE ONLY for 

execution of the state variab'le equation generator. 

A dummy STOP TIME is also required. 

For debugging purposes, PROSE allows the run control: WRITE SIMUL8 DATA 

All other run controls are illegal. 

3.3.8 End 

An END card is required at the end of the input to PROSE 
No other heading or subheading cards may be included. 

3.4 Features 

3.4.1 PROSE system output 

The output from PROSE begins with a FORTRAN routine, -Within the routine 
is a listing of the state variable equations. These equations are of the form: 
DL or DC = f (elements , I ,V ) 

J_i o 

dl dV 

where DL = — r— and DC = — ^ 
dt dt 

where I and V are the circuit state variables. Also, if any matrices appear 
- in the state equations, the matrix elements will be listed and their inverse 
elements (elements of the inverted matrix) will be specified by a trailing 
"I” in the state equations. 

Following the subroutine output, the linear circuit sensitivity analysis 
results are given. 
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The state variable names are listed, followed by the state transition 
matrix A. These are followed by the independent source names (if any) and 
the input matrix B (if any independent sources exist). 

The eigenvalues of matrix A are then displayed followed by the 
sensitivity analysis. This analysis is computed by increasing each element ■ 
value (except for independent sources) by 1% separately and recomputing the 
transition matrix eigenvalues. 

PROSE lists the sensitivities of each element to the maximum eigenvalue 
and also gives an indication of elements which may be removed (either short 
or open circuited) to eliminate a state variable from the system. 

3.4.2 Mechanical Systems 

2 

The standard SUPER*SCEPTRE mechanical description systems may be 
used by PROSE giving state variable equations for these systems. Combinations 
of circuit and mechanical descriptions may also be analyzed. 

3.^5 Limitations 

3.5.1 No equations may be specified 

3.5.2 No FORTRAN functions 

3.5.3 No linear or secondary or primary dependent sources 

3.6 Error Messages 

Error 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


RUN INITIAL CONDITIONS ONLY specified 
VECTOR EQUATIONS specified 
RUN IC VIA IMPLICIT specified 
RUN AC specified 

EY type elements exist (Linear dependent source) 

JX type elements exist (Linear dependent source) 

J8 type elements exist (Primary dependent source) 

JO type elements exist (Secondary dependent source) 

LOTaJER limit TABLE specified 

RUN MONTE CARLO specified 

RUN OPTIMIZATION specified 

RUN SENSITIVITY specified 

RUN WORST CASE specified 

EXECUTE SETUP PHASE ONLY not specified 

3 

SCEPTREl DTASAV dataset empty 
SCEPTREl PGMSAV dataset empty^ 

3 

SCEPTREl PGMSAV dataset incorrect 
No state variables exist 


See User's Manual for SUPER*SCEPTRE, AD-AOll-348 


- 17 - 



22 Pushdo™ stack overflow^ 

5 

23 Invalid symbol encountered 

24 Invalid use of XTABLE format 

25 Symbol divided incorrectly in pGjiSAV^ 

26 -■ More than 300 elements, tables, expressions, and- 

defined parameters 

27 Invalid symbol pushed into stack^ 

28 Invalid expression symbol pushed into stack 

29 Pushdovm stack overflow*^ (same as 22) 

30 A TABLE appeared incorrectly in an expression 

5 

31 Invalid use of subscript zero 

5 

32 A state variable cannot be located 

33 Stack KSHUV overflow - possible computational delay 

34 Passive element is variable - element set to zero 

35 .Matrix is too' large for display 

3 

Probably SCEPTREl syntax error 

4 

Possibly a computational delay exists 

5 

Please notify author; -those are system errors. 

All errors tjlthout note 5 are user errors. 

Note; References refer to AD catalog numbers at; 

National Technical Information Services 
U. S. Department of Commerce 
Springfield, VA 22161 


4.0 PROSE Principles of Operation 

4.1 Introduction 

The PROSE system operates in conjunction with the SUPER*SCEPTRE 
package'' . The- system utilizes a standard unmodified SUPEF*SCEPTRE 
system only as an equation generator, i.e. the execution phase (simu- 
lation)' of STJPER*SCEPTRE is not required. The PROSE system is a "post- 
processor" to the equation generator. Thus, the total PROSE system 
operates in six sequential job steps. 

4.2 Steps of the PROSE system 

The ■ first two steps of the operation of PROSE are the standard 
SUPER and SCEPTl steps of a SUPER*SCEPTRE job. The SUPER step perforins 
any processing required oh a standard SUPER*SCEPTRE input stream (i.e. 
mechanical, digital, transfer function, etc) to produce actual SCEPTRE 


6 


See previous reference AD-A011-348‘ 


- 18 - 



coding. The SCEPTl step Is also standard SCEPTRE input processing and 

7 

equation generation . . The state variable equations generated are placed 
in a FORTRAN subroutine (SIMUL8) and written on dataset PGMSAV. Also, 
data required for execution of the equations is placed on dataset DTASAV 
by SCEPTl. 

With this setup stage performed, the PROSE system begins processing 
on the PGMSAV and DTASAV files. This third step reads the DTASAV file’ 
and verifies its legality according to the criteria setup in the PROSE 
manual. 

After verification, the system searches the PGMSAV file for the 
state variable equations. The equations are pushed out in an equation 
stack .(NFUSH) with the variable name (left hand side of the equation) and 
a pointer to its equation. in NPUSH placed in array KLOCA. All expressions 
and state equations are thus saved in memory. During this process, the 
element names are also saved in array NSYM. 

The system next begins writing its FORTRAN routine. The routine is 
one which can compute the A and B matrices of- the linearized, state variable 
form: 

X = AX + BU 

where X is the vector of state variables and U is the vector of input sources. 

Prior- to the state equation generation, the system outputs the elements 
of any R,G,L, or S matrices which exist in the circuit. 

The state equations are now generated by the following process. The - 
symbol table (NSYM) is searched .sequentially for state variables, and the 
equation for the derivative of that variable is initiated. The process 
Involves location of the equation in the stack (NPUSH) , popping the- equation 
until a variable name is located, pushing the current equation pointer (onto 
KSHUV), and then recursively generating the equation for the current variable. 
At each location of a variable in an equation, either a set of parentheses 
are forced around the new equation (if the new equation has any sum terms in 
it), or not (if no sum terms exist); or, possibly, the variable is constant 
(elemenL value) or a state variable itself, in which case the variable name 
is entered in the equation. This process may form awkwardly long state 
equations but decreases as much as paracticable the number of parentheses 


^SCEPTRE SUPPORT II, VOLUIIE 11, Revised Mathematical Fromulation, AD-88-2-385 


- 19 - 



and forces only constant valued variables to appear in the equations. Each 
state variable equation is thus generated by this same procedure. ' 

The routine has now become a self-contained representation of the state 
equations of the input circuit. 

The next two steps of PROSE perform the compilation of the SIMUL8 
subroutine and linkage edit the SIMUL8 object with the object (load module 
in actuality) of the simulation portion of the PROSE system. ■ 

The simulation (GO step) reads the DTASAV file and begins the generation 
of the A and B matrices by alternately setting one and only one state variable 
to 1, all others are set to zero. Of course, all constant element values are 
set to their nominal values as specified in the input description. This A 
and B matrix generation is performed by subroutine EIGAN which also calls the 
IBM Scientific Subroutines HSBG and ATEIG to compute the eigenvalues of the 
real A matrix. Finally, the simulation performs a sensitivity analysis by 
Increasing each constant non-input element value, respectively, by 1% and- 
recomputing new A and B matrices and eigenvalues of the new A matrix. 
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5. Examples 

.Three examples are -given which illustrate the use oi the PROSE program 
in obtaining state equations for electrical networks. For each example- the 
SCEPTRE listing is given together with the circuit diagram- The state 
equations are printed as the output. Notice that 

V 

DLX represents 


and 


dV 


DCX represents 


CX 


dt 


5.1 RLC Circuit 

A simple RLC circuit is used for the first example. The circuit 
diagram is shown together with the SCEPTRE' input in Figure 5.1a. The 
state equations produced as the output of the PROSE program is also 
listed in Figure 5.1b. 

In Figure 5.2 the A and.B matrices are given. Next -the eigenvalues 


are calculated together with 'the sensitivities of the maximum eigenvalue of 
one per cent change in each parameter. If the sensitivity of the circuit 
element is within 3% of unity, the circuit element can be removed from the 
network. How the circuit element can be removed (by setting it equal to zero 
or by setting it equal to infinity) is indicated in several of the examples 
that follow. 
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TH’S IS VERSION N OF SUPERS SCEPTRE UPDATED 2/12/75 

BX3-97W58I 

.............. 


CIRCUIT CcSCRIPTION 
EXAMPLE : 

elements 

E,G-1=1 

Rjl-2=1 

L, 2-3=1 
C, 3-0 = 1 
RL, 2-G=i 
OUTPUTS 
E 

RL‘f{ CCMRCLS 
STOP tlHE=i 

execute SETUP PEASE ONLY 
ENO 


END CF PRE-PROCESSOR 





Figure 5.1a. RLC Circuit Example. 
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o n o n f ■> o o 


S IMUL AT IC-v! HAS BEEN DELETED 

C<><><><><><><><><><><><> <><><><><><><><><><><><><> 

STATE VARIAELE EQUATIONS 

OL=t-VC-U>fR+E)/(L) 

CC={-VC/{RL) + IU/IC) 

ENC 


Figure 5.1b. State Equations for Figure 



NASA-LdWIS research CtNTER 


PRCSfc FROOPAf' 

UNIVERSITY OF SOUTH FLORIDA 


EXAMPLE : 


state variable names 


IL 

VC 


MATRIX A 


-UCCCCCD+CO -l.COOO''0+OD 
l.CCCCCC+CC -l.CCCCCD+CC 


INDEPENDENT SCLRCES 
E 


MATRIX 8 


1 . CCOCCC+00 
o « c 


EIGENVALUES CF MATRIX A 
REAL PART 


-i.CCCCOODC+00 
-1. CCCCCCCO+CC 


IMAG. PART 

1. 000000'^ D+jO 
-l.CCCCOOOC+00 


SENSITIVITY ANALYSIS 

MAXIMUM EIGENVALLc 2 MAGNITUDE = 1. AIA2136D+00 

EUBMSNT SENSITIVITY 


•C 

R 

RL 

L 


-C.AS6 
0 .2 50 
-0.248 
-C.4SE 



NO SENSITIVITIES WITHIN 3« OP 1 . 

C MAY BE PEMCVEC (I.E. SET TO ZERO IF SENSITIVITY 

SET TO INFINITY IF SENSITIVITY 


< ZERO 
> ZERO.) 


Figure 5.2. Matrices, Eigenvalues and Sensitivities for KLC, Circuit. 
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5 . 2 Mutual Inductance 

A simple transformer system is shown in Figure 5.3 together with the 
SUPER*SCEPXRE input. The state equations are given in Figure 5.4. 

Notice the generation of the inductor matrix equations, XML, and the 
references to the inverse elements of this matrix, XMLI. In Figure 5.5 
the A matrix, B matrix, eigenvalues and sensitivities are listed. 
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SUPcR'^SCE FTPS SYSTEF INPUT — QI/T5 VERSION — S/360 -- OS RELcASS 21.8- 
THIS IS Vc^>SIC^ N CP SUFER*SCCPTR= UPDATED Z/lZ/75 

dev=lcp=6 ev the' lKiversity gf scuth flcrica. 

PLEASE REPORT. ANY SOP 'TWARE 'ERRORS TO OR. JAMES. C. BOWERS 8L3-97A-2?B1 


CIRCUIT CSSCRIPTIGN- 

yUrUAL lADUCTANCE TEST .FOR PROSE 

ELEMENTS 

E, 1-2=1 

Rl, 2-3=1 

LI ,3-1=1 

L2,A-==i 

H,L1-L2=C.S5 

F2,A-5=i 

CUTPlTS 

E,VLl,VL2,ILltIL2rIR2 
RUN CONTROLS 
SfO? TIME=5 

EXECUTE SETUP PHASE ONLY 
END 


END CF FRE-FRCCESSCR 


Figure 5.3. 


Mutual Inductance Example. 
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o o o o ri o n o o 


SIMULA! lOi'i HAS SEEN DhLETEO 
C 

C<><>0<>0<>0<><><><><><><><><><><><><><>0<><><><><><><><><><>0<>0< 



EKC 


Figure 5.4. 


State Equations for Mutual Inductance Example. 
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NASA-LEhlS research CENTER 


PHCJE PROGRAM 


UNIVERSITY Cf SCUTH FI.ORIDA 


MUTUAL INCLCTANCh TEST FOR PROSE 


state variable ^A^'ES 


ILl 

IU2 


.''ATRIX A 


-l.C25tAC+Ci 9.743590+00 
5.7A359f, + CC -1.C2564D+C1 


^OEPE^.CE^T SGLPCES 


MATRIX 3 


1 .C:564C + 0j. 
7425SD + C3 


higexvalles cf matrix a 


REAL PART 


-2-CCCOOCOC+Ol 
-5. I2E2C510-CI 


SS-SSITIVITY analysis 

MAXIMUM. EIGE'JVALLc H 1 MAGNITUDE = 2. C CCCOOOD+01 

ELEMENT SENSITIVITY 

R1 C.5CC 

R2 O.SCC 

LL -9.069 

L2 -9.C69 

NO SENSITIVITIES WITHIN 3? OF 1 . 

LI MAY EE REMOVED (I.E. SET TO ZERO IF SENSITIVITY < ZERO 

SET TO INFINITY IF SENSITIVITT > ZERO.) 


Figure 5.5. Matrices, Eigenvalues and Sensitivities for Mutal Inductance Example. 
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0.0 



5.3 Integral Pulse Frequency Modulation Regulator 

The circuit model for -an integral pulse frequency modulation regu- 
lator* is given in Figure 5.5. The state equations follow in Figure 
5.7. The A matrix is shown in Figure 5.8 with the eigenvalues for the 
A matrix. Also sho\ra is the variation in • eigenvalues that accompanies 
a one per cent change in each parameter. 


* R. P. Iwens, Y, Yu, and J. Triner, "Time Domain Modeling and Stability 
• Analysis of an Integral Pulse Frequency Modulated DC' to DC' Power 
Converter", PESO *75 Record 
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SUPER*5C£PTRc SYSTEM INPUT 


0?/75 VERSION 


S/360 


OS RELEASE 21.8 


i A A A V » )» 4 A A « ip* «>!■ W»i5> A## A* A**#*#*)*^**** 


Tl'IS IS V:;SSIJ.‘J iN OE SUPcR’XSCEPTRE UP0AT£D 2/12/75 
DHVELGPdC 8Y Th£ UNIVERSITY OF SOUTH FLORIDA. 

PLEASE REPORT ANY SOF rW\RE ERRORS TO OR. JAMES C. BOWERS 813-074-2581 


CIRCUIT DESCRIPTION 

r.TEGPU PFM DC TO DC CONVERTER 

IV.ENS AND YU - PESC 75 OECORO P80-90 

NASA LEWIS RESEARCH CENTER 

ELEMENTS 

El ,G-1=3C 

PO, 1-2=3. i. 15 

LOt 2-3= 2 5 ‘jE“6 

CD,3-4=3CCE-6 . 

R5f4-G=C . J?7 
RL , 3-G=lC 
R1 ,3-5=2c.Tt3 
R2 j 5-G = 13. 5E3 
C2i 3-0=0.0226-6 
P3, 5-6= 105 3 
BA ,6-7= ICC :3 

EX.7A- ' = xl ( {EI-ILO-xpO-VPL)«PN) 
ER,G-7A=X2( (R2/(R1+R2)}*PKR) 

ECP, 7A-j=XA{0.0) 

C=0P=0. 3 
J0P,G-*-' = X3(I£.lP) 

Cl ,8-G=22CvE-.2 

DEFINED PARAMETERS 

PN=C.6b 

P.<R= IS 

CU7PUTS 

61 ,VRL( :C) ,VC1 (EC) 1 1 L0( T ) , VCO { VC) 

RU J CONTROLS 

STOP Ti:i:=ico:-6 

EXECUTE SETUP PHASE ONLY 
END 


END CF PRE-PROCESSOR 


2 2 
TS) W , 



C2 



Figure 5.6. Integral Pulse Frequency Modulation Regulator Model. 
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©CO 




no o o o n 


SIMULA! [CU‘ HAS BEEN DELETED 
C 

cooo <><><><><><><><><><><><><><><><><><><><><><><> O^XXXXXXXX 


STATE V^RIABLe cQUATIOfJS 

D.0={-VC2-ILn/t l./(ROn-( ( ;R2/(R1 + R21 5*PKR) l-{0.00)+EI)/(Ln» 

c::= ( 1 + (+vc?-vcu<-n ( R2 /( ri+r2) )>!<pkr) ) + { d.ooi )/ir 5 ) + { +vc2+n ( 

SR’/(R1+R2) )-'PKR) l+(0. JO) 1/ (RL} + {+ (+1. / (R2 )*(-( ( (P2/ (R1+R2) )«PKR) )- 
)-l ./(RU -.'(-VCE) }/((l./(R3)+l./(R21+L./(Rl))) + (MR2/t'’i + 

$S2 ) I-PKR) ) + (J.0ol)/tR2) + «-n (EI-ILO'-RO-(+VC2+H lR2/(Rl+»2) )-"PKR) > + 
i(2.CL) ) ) 4(0.00) >/(R4)-IL0) ) / (Cl) 

Di;2=(-( 4VC2-VC0+( ( (R2/(R1+R2) I :<PKR ) ) + (0.00) ) / ( R5 ) - ( +VC2+ ( ( ('’2/( 

S' L + S2) )^-PKR) ) + (C.O0) )/(RL)-(+VC2~( +1./' (Raji^l-l t ( R2/ ( R1 + R2 ) ) '^PKR) )- 
S( ).C0) ) -l./(Ri) v(-VC2 ) ) / ( (1 ./ (R3 )4l./ (R2 )4l ./ (Ri) ) J )/(Rl )4ll.U)/( 
SC2) 

Di;C= (+VC2-VCO+{ ((R2/(R1+R2) )!-PKR) ) + ( 0. 00) ) /( R5 ) / (CO ) 


COOO :><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>< 

c 

ei(C 


Figure 5.7 


State Equations for Figure 5.6 



!ASi-L'W13 C.’iTCR 


P=!Oii PR J3P&M 


UMIViKSiry CF south FL0°I0A 


m=3R..'.L pf;i r*; n ctivcrte’^ 

;»£•', 5 A.NO YU * PtSC 73 ■’SCORO P90-90 
HASA LrVitS RESEARCH CcNTgR 


STATE /AJIABL; \AM£S 


ILO 

VCL- 

VC2 

VCU 


MATRIX A 


R.CICCCC+C. 

0.0 

-4.!;coocD+o3 

0.0 

<..SA545C+0o 

0.0 

5.94B630+C7 

-5.903190+09 

4 . 545 45C+07 

0.0 

-4.948660+08 

5.903190+03 

C.C 

C.3 

4.329000+04 

-4.329000 + 04 


lNU£Oc.,CiST ShURCCS 


El 


MATP i X B 


>D,0C*0 J 
“2 .YSRSSC+Oi 
0.0 
0.0 


El GE'iV.\U,2£ OF RATRIX A 


RsAL PART IMAG. PART 


C.O O.I> 

-5.T«90825U+''8 ''.0 

-3.AB229790+02 3.623A0A6C+03 

-3. AH22979U402 -3,62390460+03 


SSRSITIVITY A-H.LYSIS 


'AX I-'+U > 
LE.MENT 

EIGEN/SLUS 

SENSITIvn 

EX 

0.0 

2 R 

:.c 

c: 

J.C 

C2 

-■).990 

CO 

-C.CCO 

R3 

.c:o 

R3 

-3.C39 

R5 

-0.9B3 

RL 

-J.COB 

?2 

-) .009 

R1 

-0.CC3 

R4 

O.C 

LD 

C.CCO 

JOP 

0.0 


2 KAGNITUne = 5.94908200+08 

«*«■ CtEIlENT MAY PE REMOVED BY SETTING TO 2 frc). 

*+>* ELEMENT MAY BE REMOVED BY SETTING TO 7FR0. 

Figure 5.8.' Matrices, Eigenvalues 
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Sensitivities for .Figure 5.6. 



6. Summary and Conclusions 


The- work during the two years of the NASA Grant NSG-3096 has resulted 
in a computer program which will convert a SCEPTRE listing for an electrical 
or mechanical system into a set of state equations for the network. Additional 
features of the program to PROd uce ^tate Equations include subroutines to: 

1. linearize the state equations about an operating point as 
specified by the initial conditions 

2. obtain the eigenvalues for the linearized state equations 

3. obtain the sensitivities of the eigenvalue with the maximum 
magnitude 

4. use the sensitivity calculations, to indicate which circuit 
element can be removed from the network without adversely 
affecting the systems dynamic operation 

At the end of the first year the PROSE program had limitations on the type of 
networks it would' accept. In the second year several of these limitations- 
were removed at a great cost in programming manpower. The most significant 
improvement is the ability PROSE now has to accept networks with mutual 
inductance. 

Additional research and development could profitably be spent in applying 
the PROSE program to the design of automatic control systems- Using SUPER*SCEPTRS 
an electromechanical description of the plant to be controlled would be reduced 
to the state equations for the systems. The state equations could be used in 
an algorithm to automatically compute the control matrix. 
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Appendix 


This appendix gives the Fortran listing for the computer program that 
PROd uces ^tate Equations. 
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FORTR/SN IV G1 RELEASE 2.0 


MAIN 


DATS = 78178 


1A/4A/45 


PAGE 0001 


CCCl 


CCC2 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PP.CSSl PRl 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPfU 

PRi 

PRl 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRl 


COPYRIGHT 1977 BY JCHN E. O'REILLY, JR. PRl 

10/18/77 ■ PRl 

CCPYRICFT 1976 EY JOHN E. O'REILLY, JR. PRl 

05/13/76 . PRl 

PRl 

PRl 

TRICE.\T£NMAL -VERS ION OF P-R-O-S-c PRl 

PRl 

* " PR’ 

^ * <*//////;////////////////////////////// . PRl 

* X. # # * ///////////////////////////////////// PRl 

pRi 

///////////////////////////////////// PRl 

*///////////////////////////////////// PRl 

# v si! ■!« PRi 

/ n / It n j n / n t u u uu / ! n n n/ u f pri 

**>:'»« n/tnun/n ininiinnniitntiin, pri 

PRl 

/////////// ///////////////////////////////'/////////// PRl 

( n n ! ut m u t u n t ! n ! 1 1 n n n ! I n ! I u t ! t ! I ! n m u pri 

PRl 

///////////////////////////////////////////////////// PRl 

.///////////////////////////////////////////////////// PRl 

PRl 

7//////////////////////////////////////////////////// . PRl 

I n ! I i > ! n ! ) I rn n t ! I n ! 1 1 n n N > n } ! n ! I ) n n t > n / N pri 

PRl 

//////////////■/////////////////////////////////////// PR 1 

/■////////// ////////////////////////////////////////// PRl 

PRl 

PRl 

BLOCK DATA . PRl 


<><><><><><><><> <><><><><><><><><><><><><><><><><><><><><><><><><>PR1 


PRl 

PRl 

P-R-O-S-E - GEAERA’TE STATE EOLATIONS FROM SCtPTRE PRl 

BUILT PGKSAV AND DTASAV FILES.’ THE ROUTINES PRl 

COMPUTE THE STATE VARIABLE EOLATIONS AND A AND B PRl 

MATRICES. PRl 

PRl 

PRl 

OOO^rXXXXXXXXXXXXXXXXXXXJOOOOOOOOOOPRi 
IMPLICIT REAL^P lA-H ,0-ZI PRl 


CCCl 
OCOZ 
0003 
0G04 
CCC5 
CCC6 
OU07 
0CC8 
CCC9 
0010 
0011 
CC12 
0C13 
00 lA 
0015 
CC16 

0017 

0018 
CC19 
0C2C 
0021 
002 2 
0023 
002A 
0025 
CC26 
0027 
0023 
CC29 
CC3C 

0031 

0032 
CC33 
003A 
0035 
GC36 

0037 

0038 
0U39 
OOAC 

0041 

0042 
C'G43 
CCAA 

0045 

0046 
CC47 
0046 
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FORTRA^ IV G1 release 2.0 


RLK DATA 


DATE = 78i.7a 


14/A4/45 


PAGE 0002 



C 




PRl 

CCAS 


C 




PRl 

CC5C 


C 

VARIABLE 

DEFINITIONS - NASA 

PRl 

0051 


C 

NSYH 

- 

ELEMENTS SYMBOL STORAGE ARRAY 

PRL 

005 2 


C 

KSYM 


•NURBE'’ CF SYMBOLS IN NSYM 

PRi 

CC53 


C 

MP^AY 

- 

INPUT AND OUTPUT 80 BYTE ARRAY 

PRi 

0054 


C 

KRAY 


’CURRENT loc;ation pointer in NRAY 

PRl 

0055 


c 

NLEN 

- 

ARRAY CF lengths OF EACH TYPE OF ELEMENT 

PRl 

CC56 


c 

NPT 

- 

POINTERS TO ELEMENT TYPES IN NSYR. 

PRl 

CC57 


c 

NT 


TYPE OF ELEMENT VALU£(=0 IS NOT £J TYPE) 

PRl 

005R 


c 



• 

PRl 

CC55 

0CC3 


COMMON 

/NASA/ KSYM(6t300) , KSYR , NRAY ( 80 ) ,KRAY , NLEN( 8 ) » NPT ( 8 ) , 

NT (300PR1 

0C6C 



1 » 



PRi 

0061 


c 



PRl 

CC62 


c • 




PRl 

CC63 


c 

variable 

DEFINITIONS - PUSHED 

PRl 

0064 


c 

KLOCA 

- 

array of left HAND SIDE ENCODED NAMES AND , 

PRl 

0065 


c 



LOCATION OF EOUATIGN IN NPUSH AND 

PRl 

CC66 


c 



SUR. TERMS FLAGS 

PRi 

0067 


c 

KLOC 

- 

LENGTH OF KLOCA 

PRI 

006 8 


c 

NPUSH 

- 

ECUATICN ARRAY 

PRl 

CC69 


c 

KPUSH 


LENGTH CF EOUATICN ARRAY NPUSH 

ORl 

0C70 


c 

KPMAX 

- 

MAXIMUR length OF NPLSH(LCOO) 

PRi 

00 71 


c 



• 

PRl 

0072 

0CC4 


CDPMON 

/PLShEO/ KLOCA (3 »400 )* KLOC, NPUSH 11000 ) r KPUSH, KPM AX 

PRl 

OC73 


c 



. 

PRl 

0074 


c 




PRi 

0075 


r 

VARIACLE 

DEFINITIONS - SHGR 

PRl 

CC76 


c 

NSH’JV 

- 

PUSH OCHA STACKtFIFC TYPE) 

PRl 

CC77 


c 

KShUV 

- 

CURRENT POINTER INTO NPUSH STACK 

PRl 

0078 


c 

KST • 

- 

START CF CURRENT EQUATION IN NPUSH 

PRl 

CC7S 


c 

KiiMD 

" 

END CF CURRENT ECUATICN IN NPUSH 

PRl 

ccao 


c 

KSU8 

- 

current VARIABLE'S SUBSCRIPT VALUE 

PRl 

0081 


c 

KEXTD 

- 

CURRENT VARIABLE’S EXTENDED VALUE 

PRl 

00B2 


' c 

KFAR 

- 

PARENTHESIS FLAG 

PRl 

0083 


c 




PRl 

0C84 

CC05 


CCPHCN 

/SHGR/ NSHUV(500),KSHUV,KST,KcND,KSUB,KEXTD,KPAR 

PRl 

0G35 


c 




PRl 

CC66 


c 




PRl 

0087 


c 

variable 

DEFINITIONS - CNTRL 

PRl 

0088 


c 

LIST 

- 

ccbug flag 

PRl 

0039 


c 

LI 

- 

RESERVED WORD 

PRl 

CC9C 


c 

CTROL 

- 

CONTROLS FROM DTASAV 

PRl 

0091 


c 

KXK 

- 

XML, XMRC, AND XMS ARRAY SIZES 

PRl 

0092 


c 



• 

PRl 

CC93 

0CC6 


C ORMON 

/CNTRL/ LI ST, Li , CTROL (160) ,KXM(3) 

PRl 

CC94 


c 




PRl 

0095 


c 




PRl 

0096 
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FORTRAN IV 31 RSLSASF 2.0 


DLK DATA 


DATE = 73178 


1A/AV45 


PA3E 0003 


CCC7 • 


0CC3 • 


II 

% 


% 


O' 


c 

VARIABLE 

DEFINITIONS - FILES 

PRi 

0097 

c 

MPRINT 

- 

PRINTFP. FILE NUMBER 

PRl 

0098 

c 

NPGK 

“ 

SCEPTRE PROGRAM FILE NUMBER 

PRI 

CC99 

c 

NDTA 

" 

SCEPTRE DATA FILE NUMBER 

PRl 

OlOC 

c 




PRI 

0101 


CCNMCN 

/FILES/ NPP.IiNT,NPGM,NDTA,NEQU 

PRl 

01C2 

c 




PRi 

0103 

c 




PRl 

OlOA 

c 

VARIABLE 

DEFINITIONS - CHARS 

PRl 

0105 

c 

NSP 

- 

SPACE 

PRl 

01C6 

c 

'JLP 

“ 

LEFT PARNTHESIS 

PRl 

0107 

c 

NRP 

- 

RIGHT PARENTHESIS 

PRi 

0108 

c 

NCCM 

- 

COMMA 

PRi 

01C9 

c 

NtQ 

* 

EQUAL Sl’CN 

PRl 

CllC 

c 

.NO 

- 

ALPHA LETTER "0" 

PRi 

0111 

c 

NPLUS 

- 

PLUS SIGN 

PRi 

0112 

c 

NPINUS 

~ 

MINUS SIGN 

PRi 

0113 

c 

NSL 

- 

SLASH 

PRi 

OllA 

c 

NZcR 

- 

ZERO NUMERIC 

PRi 

0115 

c 

NXM(3)- ALPl-A 'L'f'R't'S' 

PRi 

0116 

c 

N'UP(IO) 


NUMERICS 

PRl 

0117 

c 




PRi 

0118 


CCKMO;<l 

/CHARS/ NSP, NLP tNRP,NCOM,NEQ,ND, NPLUS, NMINUSi 

,NSL,NZER,NXM(3PRl 

0119 


U ,NUK(10) 

PRi 

C12C 

c 




PRl 

0121 

c 




PRi 

012 2 

c 

VARIAGLE 

DEFINITIONS - HOLDIT 

PRl 

0123 

c 

HEAD 

“ 

HEADI.SC CARO BUFFER 

PRl 

C12A 

c 

KHEAD 

- 

NUMBER HEADING CAROS 

PRl 

0125 

c 

N CUM 1 

- 

RESERVED WORD 1 

PRi 

0126 

c 

ELTS 


ELEMENT VALUES 

PRi 

0127 

c 

KELTS 

- 

NUMBER CF ELEMENTS 

PRl 

0128 

c 

K£L . 

- 

NUMBER ELSMENTS-NUMBER £J ELEMENTS 

PRi 

0129 

c 

SCR 

“ 

SOURCE DERIVATIVE VALUES • 

PRl 

013C 

c 

NSOR 

- 

SOURCE DiiRIVATIVE NAMES 

PR 1 

0131 

c 

KSDR 


NU.MlsE° SOURCE DERIVATIVES 

PRl 

0132 

c 

NCUM3 

“ 

RESERVED WORD 3 

PRl 

0133 

c 

XI c 

- 

INITAL STATE VARIABLE VALUES 

PRi 

013A 

c 

KSV 

“ 

NUMBER STATE VARIABLES 

PRl 

0135 

c 

NCUMA 

“ 

RESERVED WORD A 

PRl 

0136 

c 

pak 

- 

CEFINEC parameter VALUES 

PRl 

0137 

c 

NPAR 

- 

0£FIN=D PARAMETER NAMES 

PRl 

0138 

c 

KDP 

- 

NUMBE"’ DEFI.NED PARAMETERS 

PRl 

0139 

c 

KCUM5 

“ 

R35£R'/ec WORD 5 

P.Rl 

OlAC 

c 

XTAB 

" 

TABULAR FUNCTICN VALUES 

PRi 

OlAl 

c 

NTAB 

“ 

TABLE NAMES 

PRl 

01A2 

c 

KTAB 

- 

NUMBER ENTRIES IN XTAB 

PRi 

01A3 

c 

KTA32 . 

- 

NUMBER CF TABLES 

PRl 

OlAA 
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C 

KEJ - ^UPBSR CF INPUT (INCEPcNDENT) SOURCES 

PRl 

0145 


C 


PRi. 

0146 


c 

. 

PRl 

0147 

0CC9 


COPKON /HOLCIT/ HEAD ( 9, 1 1 1 , KH5A0 t NOUK 1 , E L TS ( 300 ) , KE LTS , K£ 

L,SDRI50)PR1 

0148 



1 ,NSOR (6,5Vi ) ,KSCR,NDUM3,XIC(50 ),KSViNOUM4,PAR('100),NPAR( 6i 

100) .KOPfPRl 

0149 



2NDUM5,XTABI 500) ,ATAB(7,50J , KTAB ,KTAB2 ,KEJ 

PRl 

015C 


c 


PRl 

0151 


c 


■PRl 

0152 


, c 

SET CONSTANTS 

PRi 

0153 


c 

PU.SHED 

PRl 

0154 

0010 


DATA KLOC/0/ 

PRl 

0155 

cell 


DATA KPUSh/l/,KPMAX/1000/ 

PRl 

0156 


c 

FILES 

PRl 

01-57 

0012 


DATA NPRINT/6/,NPGM/l/,NDTA/2/iN£0U/3/ 

PRl 

0158 


c 

CFARS 

PRl 

0159 

0013 


DATA NSP/iH / , MP/1 H ( / , NRP/ 1 H ) / , NCOM/1 H, / 

PRl 

016C 

OO’IA 


DATA NEO/iH=/,ND/lHD/,NPLLS/lH+/fNMINUS/lH-/ 

PRl 

0161 

0015 


DATA NSL/lH//,NZER/lHa/ 

PRl 

0162 

ecu 


DATA .SXP/lHL.lFRjiHS/ 

PRl 

0163 

0017 


DATA .NLM/lHC,lHl,lH2,lH3,lH4,iH5,iH6,lH7tlHa,lH9/ 

PRl 

0164 


c 

FOLDIT • 

PRl 

0165 

GCia , 


DATA KTAe/0/,KTAB2/0/,i<SDR/0/,i<ELTS/0/.KDP/0/ 

PR 1 

0166 


c 


PRi 

0167 

ocio 


DATA KXM/3=i'C/ 

PRl 

0163 

CC2u 


END 

PRl 

0169 
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PAGE OCO 


C <5<><><>0<><>0<><><>0<><>'<><><><><><><><><><><><><><><>0<><>OPRI 



c 

eEGIN PAIA ROUTINE 

PRi 


c 


PRi 


c 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRi 


c 


PRI 

CCCl 


IMPLICIT REALMS lA-H, 0-2) 

PRi 

CCC2 


COPMC.V /NASA/ NSYM(6, 300) ,KSYf',NRAY(80) , KRAY 

,NLEN(3 ),NPT (8 ) ,NT {300PR1 



1 ) 

PRI 

CC03 


CCPMON /CNTRL/ L 1ST , L 1 , CT ROL ( 1 60 ) , K XM ( 3 ) 

PRi 

CCC4 


CONMON /PILES/ NPRINT,NPGP,NCTA,N£OU 

PRI 


c 


PRi 


c 

GET SYMEOLS 

PRI 

CCC5 


CALL SETUP 

PRI 


c 


PRi 


c 

• 

Pftl 


c 

GENERATE STATS EQUATIONS HEADER 

PRI 

CCC6 


CALL PTCUP 

PRi 


c 

BUILD,, EQUATIONS . 

PRi 

0007 


CALL eiLD 

PRI 


c 


PRi 


c 

• 

PRI 


c 

FIN ■ 

PRI 

CC03 


CALL FINIS 

PRI 

CCC9 


STOP 0 

PRI 

OCIC 


END 

PRI 


II 

G > 

■ fea 


017C 

cm 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 
0101 
0182 

0183 

0184 

0185 

0186 

0187 

0188 
0189 
019C 

0191 

0192 
C193 
0194 
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CCCi 


CCC2 

OOOJ 

ccc^ 

CCC5 

0006 

CC07 

CCCH 

OCO'} 

OCIO 

cell 

CC12 

•0013 

CC14 

CC13 

0016 

CC17 

ecu 

00 U 
CC2-0 
CC21 
CC 2 2 
0023 
CCi^ 
CC2: 
0026. 
0C27 
QC2e 
CC2-? 
CC30 
CC31 

i'ri7 

0033 

CC3A 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


1 


2 

3 


A 

5 


6 


7 


SLBROLTINt BACK {.\A1<E ,I\CHAR') ' i 

PRl 

<><><><><><><> OOOOOOOOOOOOOOOOOOOOOOOOOO PPi 

PRl 

CHAi>IG£ NAPE TO KCHAR IN MECHAMCAL FORMAT PRl 

CFANGES PP'l 

GP=TP CXXX=ZXXX AND CHANGE XXX ‘ PRl 

I=F EXCEPT ICW=TJ ' PRl 

V=V EXCEPT VCU=UJ PRl 

ELEMENTS PPl 

CU=M EU=U- JU=R LU=K RU=D PRl 

CVv=J EW=G, JVi=C PRl 

PRl 

<><><><><><><><><><><><><><><><><><><><><><>< ><><><><><><><><><><> PPl 

PRl 

DIMENSION NEU5), NELF(8)f NST(2), NSTF(A), NCHARI6), NAHE(6} PRi 

EQUIVALENCE ( NG ,NELF ! 7)) t I ND ,NELF ( 5) ) » (NT.NSTFO)) PRl 

6CUIVALENCE (NU,NUF{2n, (NCtNELF ( 8 ) ) , (Ni<,NSTF C A ) ) PRl 

DATA Nf L/1HC,1H=,1HJ,1HL,1HR/ PRl 

DATA nelf/:hm,ihu,ihr,ihk,ihd,ihj,ihg,ihc/ PRl 

DATA i>jST/lHI, IFV/ ’ PRl 

DATA .\STr/lhF,!FV,lPT,lHW/ PRl 

data NZ/1HZ/,NP/1HP/,NDLX1H / , PRl 

KST=1 ‘ PRl 

KN=0 • PRi 

DO 1 1=1,6 • PRi 

IFINAMEI I) .NE.NBL) KN=KN+1 PRl 

NCFARCI ) = NAME(U PRl 

■IF(XN.LE.l) RETURN PRl 

IF(NChAR( 1) .NE.NG.0R.NCHAR(2) .NE.NP) GC TO 3 PRl 

NCFAR(1)=NT PRl 

RETURN PRl 

DO A KTNS=1,A PRl 

KL0C=KN-KINS+1 ' PRl 

IFIKLCC.LS.l) RETURN PRi 

IFIMCHAR(KLCC) .EC.NU.OR.NCHAR(KLOC),EQ.NW) GO TO 5 PRl 

CONTINUE PRl 

RETURN PRl 

IF(.\CHAR(1 J.NE.ND) GO TO 6 PRl 

NCHAR(1I=NZ PRl 

KST=2 PRl 

DC 7 Kl=l,2 PRl 

IF1NCHAR(KST).NE.NS7(KU) GC TO 7 , PRl 

NCHAR (KST) = NSTF(KU • PRl 

KST=KST+1 . . PRl 

IF (NCriARIKSTJ.EC. NC. AND. NCHAR (KLOC) . EO.NW ) NCHAR {KST-U=NSrF (K 1+2) PR I 

GO TO 8 PRl 

CCNTINUc PRl 


018 5 
0196 

0187 

oue 

0188 
O20C 
0201 
0202 
0203 
02-:a 
02C5 
02G6 

0207 

0208 
G2C8 
0210 
0211 
0212 
0213 
021A 
0215 
C2U 

0217 

0218 
0218 
022C 
0221 
0222 
0223 
022A 
02 Z 5 
0226 

0227 

0228 
0228 
023C 

0231 

0232 

0233 
023A 

0235 

0236 

0237 

0238 
0238 
02AC 
02AI 
02A2 
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C035 

8 

K£L=5 

CC36 


KELF=0 

CC37 


IFI.NCHARIKLCC) .EC.NU) GO TO 9 

003S 


KEL=3 

0039 


KcLF=5 

CC4-, 

9 

OC 10 Kl=i,KEL 

0041 


JF(NCHAR(KST).NE.NEUK1) ) 00 TC 10 

0042 


NCl-AR (KST ) = NELF(K1+KSLF ) 

CC43 


GC TO 11 

CC44 

1C 

CUNT! NOE 

0045 


RETURN 

CC46 

11 

IF(KLCC.EC.KN) GO TO 13 

CC47 


KELF=KNj:KLOC ’ 

0043 


DO 12 Kl=i,SELF ' 

CC49 

12 

NChAR (K.L0C+Ki-1)=NCHAR(KL0C+Kl) 

CC50 

13 

NCHAR IKK)=NEL 

0C51 


GO TO 2 

CC52 


END 
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PRl 0243 
PRl C244 
PRl 0245 
PRl 0246 
PRl 0247 
PRl 024S 
P.Rl 0249 
PRl 0259 
PRl 0251 
PRi 0252 
PRl 0253 
PRl 0254 
PRl 0255 
PRI 0256 
PRi 0257 
PRl 0258 
PRi 0259 
PRl 0269 




FORTRAN IV GI REUSASH 2.0 


CCDOUT 


DATE = 78178 


1VA4/A5 


PAGE OCOi 


OCCi 


CCC2 


CCC3 

CCCA 

0005 

0CC6 

CCC7 

0003 

CCC9 

CCIC 

ecu 

0012 

CC12 

cc:a 

0013 


c 

c 

c 

c 

c 

c 

c 


c 


SUBROUTINE BCOOUT (KVAL) ' PRl 0261 

PRl- 0262 

<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>PR,l 0262 

PRl 0266 

PUT ThE BCD CHARS FOR IKVALJ INTO OUTPUT STRING PRl 0265 

PRl 0266 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRl 0267 

PRl 0268 

CCRRCN /CHARS/ NSPf NLP,NRP,NC0M,N50,N0tNPLUS»NMINUS,NSLf NZ£R,NXM( 3PR1 0269 


J,KNUNS(1C) 

PRl 

027C 


PRl 

0271 

PVAL=KVAL 

PRl 

0272 

NLcN=INT( ALCGIO (AMAXO (MVAL,1 ) ) J +1 

PRl 

0273 

M01V=iC«*{Nl.EN-l)’ 

PRl 

0276 

KFL=0 

PRl 

0273 

be 1 I=1 tMEV 

PRl 

0276 

MPTR=MVAL/HOI V+1 : 

PRl 

0277 

IF(MPTR,tC.l.ANC.KFL.EO.O> GO TO i 

PRl 

0278 

CALL FUTCHR (KNUKS ( NPTR ) ) 

PRl 

0279 

kfl=i 

PRi 

0230 

MVAL=MVAL- {MPTR-iJ*MOI V 

PRl 

02 81 

MCIV = .''C1V/10 

PRl 

C232 

RE TURN 

pr; 

0262 

END 

PRi 

0286 
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OCOl 


SUBROLTINE OILD 

PRl 

0235 


C 


PRl 

0236 


C 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOP'^l 

0287 


C 

. 

PRl 

0288 


C 

BUILD THc stats VARIABLE SQUATICNS 

PRl 

0287 


C 


PRl 

0290 


C 

<><><><><><><><><> <><><><><><><><><><><><><>00 OOOOOOOOOP71 

0291 

CCC2 


IMPLICIT RS4US(A-H,0-Z1 

PRl 

0292 

CC03 


COMMON /NASA/ NSYM t 6,300 ,KSYM,NRAY( 80) ,KRAY ,M = N (8 ) ,KPT (8 ) , 

NT (30GPR1 

0293 



1 ) 

PRl 

0294 

CCCA 


COMMON /PL-S.HEO/ XLDCA(3,400) ,KLOC,NPUSH(LGOO)»KPUSH,KP.V!AX 

PR 1 

0295 

0005 


C3MM0-I /CHARS/ NSP ,NLP ,NRP ,RCCN ,NcQ, ND , NPLUS, NM INUS ,NSL iNZER 

,NXM(3PRi 

0296 



1 ) ,NUM (10 ) ■ 

PRl 

0297 

CCC6 


CCPMCA /CATRL/ L 1ST , LI , CT ROL ( 160 ) , KXM ( 3 ) 

PRl 

029 8 

0007 


COMMDn /files/ .\PRINT,NPGM,NDTA,N£QU 

PRl 

0299 

CCC3 


CIMSNSICN MPAV16), LXM(6) 

PRl 

0300 

CCC9 


DATA PRAY(1)/1FC/ 

PRl 

03C1 

CCiC 


data LXM/1HX,1HM,1H ,1H / 

PRl 

03C2 

0011 


KFL = 0 

PRl 

0303 


C 

CLEAR CUT array ONCE ONLY 

PRl 

03C4 

CC12 


DO 1 1=1,72 

PRl 

C30S 

0013 

1 

MRAY( I )=NSP 

PRl 

0306 

0014 


K R AY= 7 

PRl 

0337 


C 

OUTPUT ANY MATRICES 

PRl 

03C8 

CCi5 


OC 4 1=1,3 

PRl 

03C9 

0016 


!F(KXM(U. = 0.C1 GO TO 4 

PRl 

0310 

CC17 


IF(KFL.EC.O) WRITE (NEQU,2) 

PRl 

0311 

ccis 


KFL = 1 

PRl 

0312 

0017 

2 

FORMAT (22HC MATRIX E CUA TI CNS , 58X , / , 1 HC ,79X ) 

PRl 

0313 

CC20 


WRITE (N£CU,3) NXK ( I) ,KXM( I) ,KXM ( I ) 

PRL 

0314 

CCZl 

3 

FORMAT (CXtIEHCIMENSION X R , A1 , IH ! , 13, 1H, , 13, 1H),52X) 

PRl 

0315 

0022 

4 

CONTINLE 

PRl 

0316 

0023 


DO 6 1=1,3 

PRl 

0317 

CC24 


I F (KxXl ) .EC.O ) GO TO 6 

PRl 

0318 

CC25 


JN=KXM(I ) 

PRl 

0319 

0026 



PRl 

0327 

CC27 


L,XM{3 )=NXM(I) 

PRl 

0321 

CC23 


OC 5 J=1 , JAK 

PRl 

0322 

002? 


NB=l+( J-D/JN 

PRl 

0323 

0030 


NA=J-tNe-U4JN 

PRl 

0324 

CC31 


I F INA.GT.N3) KV AL=100000* I + (N A-1 ) =t«JN+NS 

PRl 

0325 

CC32 


IF(NA.LE.AB) KVAL=:CC0004I + (M3-1 j-KUN+NA 

PRl 

0326 

0033 


CALL PUTNAM (LXM) 

PRL 

0327 

CC34 


CALL eCCCUT (NA) 

PRl 

0328 

CC35 


CALL PLTCHR INCCM) 

PRl 

0329 

0036 


CALL aCDOLT (NO) 

PRl 

0330 

0037 


CALL PUTCFR (NRP) 

PRl 

0331 

CC3S 


CALL PUTCHR (NEC) 

PRl 

0332 
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0C39 


CALL POP (KVAU 

PRl 

0333 

lC^C 

5 

CALL FUTLIN 

PR 1. 

0334 

CCAi 

6 

CONTI ^L'E 

PRl 

0335 

0042 


IFIKFL.EO.l) V^RITE (N£QU,7) 

PRl 

0336 

0043 

■7 

.FORMAT (1FC,79X,/,1HC,79X,/, 1HC,79X) 

PRl 

0337 

CC44 


'write (NEC.U,8) 

PR 1 

0338 

0045 

8 

FORMAT (3CHC STATE VARIABLE EQUAT I CAS ,50X ,/ , 1 HC,79X ) 

PRl 

0339 


C 

LCOP CN STATE VARIABLES 

PRl 

0340 

CC46 


DC 11 KSVPT=3,4 

PRl 

0341 

CC47 


KSV=I\LENIKSVPT1 

PRl 

0342 

0043 

• 

IF(KSV.EO.O) GO TO 11 

PRl 

0343 

CC49 


LSV=NPT(KSVPT)-1 

PRl 

0344 

0C5C 


DO 10 JSV=1 ,KSV 

PRl 

0345 


c • 

POINT TO STATE VARIABLE NAME 

PRl 

0346 

0051 


LJSV=LSV+JSV 

PR-i 

0347 

✓ 

c 

PUT THE NAME 

PRl 

0348 

CCS2 


DO 9 1=1,5 

PRl 

0349 

00S3 

9 

MRAYt I + 1)=MSYM( I,LJSV) 

PRl 

0350 

CC54 


CALL PUTNAM IMPAY) 

PRl 

03 51 


C 

AND AN ECLALS SIGN 

PRl 

0352 

0055 


NRAY(KRAYJ=NEfl 

PRl 

0353 

0056 


KRAY=KRAY+1. 

PRl 

0354 


c 

DECODE TP.E NAME 

PRl 

03 5 5 

0CS7 


CALL DECMAM ( NS YM t 1 , L JSV] , 0 , J LCC ,KVAL ) 

PRl 

0356 


c 

• MAKE HER LOOK LIKE STATE VARIABLE DERIV. 

PRl 

0357 

CC5S 


KVALSV=KVAL+10900 

PRl 

0358 


. c 

POP TbE ECUATICN FDR THIS DERIV. 

PRl 

0359 

0059 


CALL POP {KVALSV3 

PRL 

0360 


c 

FORCE TFE LINE OUT . 

PRl 

03 61 

CC6C 

10 

CALL PUTLIN 

PRi 

0362 

0C61 

11 

CONTINUE 

PRl 

0363 

0062 


RETURN 

PRl 

0364 

0063 


END 

. PRl 

0365 
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OCOi 


CC02 

CCC2 


0CC4 

0005 

0C06 

CCC7 

C003 

0CC9 

CCIC 

cci: 
0012 
CC13 
CCI A 

0015 

CC16 

0017 

ooia 

C C 1 5 
0C20 
0021 
CC22 

0023 
CC2 A 
CC25 

002o 

CC27 

0028 

CC29 

0030 

CC31 

CC32 
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SUBROUTINE CECNAM (NAME >K SUB ,KLQS , KVAL ) 


PRl 

0366 

c 



PRI 

0367 

c 

<><><><><><><><><><><><><><><><><><><><><><><><> OOOOOOOOOPRi 

0368 

c 



PRl 

0369 

c 

DECODE tMAHE) INTO SUBSCRIPT LOCATION AMO 

VALUE 

PRl 

0370 

c 



PRl 

03 71 

c 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRi 

0372 

c 



PRl 

0373 


IMPLICIT . REALMS (A-H,a-Z) 


PRl 

03 74 


COMMOM /NASA/ NSYR 16 ,300 1 , KSY M , NRAY (80 ) t KRAY, NLEN 1 8 ) 

,NPrI6),NTC300PRI 

0375 


1] 


PRl 

0376 


CCM.MOY /CFARS/ NSP , NLP ,'NR P, NCQM, NEU ,N0 ,NPLUS,NMINUS»NSL ,N2ER ,NKM( 3PR1 

0377 


1) ,i\UMIlC) ' 


PRl 

0378 


CONND'J /CNTRL/ LIST.LI ,CTROL( 160) ,KXK(3: 


PR 1 

0379 


CCMMO'l /FILES/ MPRINT,NPGM,NOTA ,N60U 


PRl 

0380 


DIMENSICM NANS(6), N£LTS!7), MTYPE(6) 


PRl 

03 81 


DIMENSION NXTAPUS) 


PRl 

0332 


DATA MrLTS/iHC, IhR, IHL , IHE , IHJ , IHP / 


PRI 

0383 


DATA A'rYP£/lKD,lhV,lHI,lHS,ll-G,lHX/ 


PRl 

0334 


DATA AN/1HN/,\M/1HK/ 


PRl 

0335 


DATA NXTAB/ IHT, IHA, IHB, IHL, IHE/ 


PRl 

0336 


ECUIVALEA’CE IN I , NTY PE ( 3 ) ) , (NP,N£LTS( 6 ) 1 


PRL 

0387 


EOUVALENCE {MXTABd) , MELTS {?)) , (NCHART, 

NXTAB(I) ) 

PRl 

0388 

c 



PRl 

0339 


KLQC=0 

, 

PRl 

0390 


KELTS=0 


PRl 

C3 91 


KTYP£=0 


PRl 

0392 


KEXTD=0 


PRl 

0393 

c 

LOCATE SY.MfiOL TYPE 


PRl 

0394 


DC 1 1=1,6 


PRl 

0395 


IF(NAHE(ll.EO.NTYPEm) GC TO 7 


PRl 

039-6 

1 

CCNT IaUE 


PP.l 

0397 


KST=l 


PRl 

0398 

c 

LOCATE SYMBCL ELEMENT 


PRl 

0395 

2 

DO 3 1=1,7 


PRl 

0400 


IF<NA.''E<KST).eC.N£LTS(I) ) GO TO 8 


PRl 

04CI 

3 

CONTI MUE 


PR 1 

04C2 

C 

MUTUAL TYPE7 


PRl 

0403 


IF(NAR’;(KST).Ne.NM) CO TO 4 


PRl 

04 34 


KELTS=9 


PRi 

C4C5 


GO T3 ? 


PRl 

0406 

C 

VCLTAC-'^/CURRENT TYPE 


PRi 

04u7 

4 

IF(KTYPe.£G.2.CR.KTYPS.e0.3) GO TO 12 


PR 1 

0408 

C 

fork SYMBCL VALLE 

. 

PRl 

04C9 

5 

KVAL = KLOC + lCGC’i‘KELTS+lCCOC*KTVP=+lCCCCC«KEXTD 

PRI 

0410 


IF(LlST.SQ.l) WRITE (NPRINT,6) NAME,KLOC, 

KELTS.KTYPE 

,KcXTD PRi 

0411 

6 

FORMAT <6H DECNA«= ,6 A1 ,6K KLCC=,I3,7H KELTS=,I3,7H KTYPE=,r3,7H XEPRl 

0412 


1XTD=, 13) 


PRl 

0413 


6 2 
y-O 2 

SI 

ci ^ 

^ a 
9 
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CC23 


RETURN 


PRl 


c 

GOTTA TYPE NAME 


PRI 


c 

XPRSSSION IS-SOECIAL 


PRl 

CC34 

7 

IF(I.cG.6) GO TC 17 


PRl 


c 

NOT AN ''X"-TY?6 


PRl 

CC35 


KST=2 


PRl 

CC36 


KTYP£=I 


PRl 

CC37 


GO TO 2 


PRl 


c 

GOTTA ELEMENT 


PRl 

OC35 

3 

K£LTS=I 


PRL 


C 

PARAMETER IS SPECIAL 


PR 1 

CC39 


IFIKELTS.EO.?) go TC 5 


PRl 


C 

LOCATE NANE IN SYMBOL TABLE 


PRl 

CC40 

9 

CO 11 1=1 ,KSYM 


PRl 

CC4i 


00 10 J=l,5 


PRl 

0042 


KSTJM=KST+J-1 


PRI 

C04J 


IF(.NAMC(KSTJM).NE.NSYMfJ,n) GO TO 11 


PRl 

CC44 

10 

CONTINUE 


PRl 


C 

FOUND NAME IN STABLE 


PR 1 

CC4 5 


KLOC=I 


PRl 

CC46 


GC TC 5 


PRl 

CC47 

11 

CONTINUE 


PRl 


c 

NOT IN STABLE OR V/I TYPE 


PRl 

CC43 

12 

KSTAR=2 


PRi 


C 

V4'/I# TYPE 


PRl 

0049 


IF (LPLMPlNAMEjKSTAR ,1 ,KSTYU .EC. I) GC 

TC 5 

PRi 


c 

■ VflC/l;/C TYPE - STATE VARIABLE DERIV, 


PRl 

CC50 


I F (NAPrIKSTARJ . EC.NC) GO TO 15 


PR 1 


c 

ViiP/IYP TYPE 


PRi 

0051 


IFINAME (KSTAR) , EQ.NP) GO TO 16 


PRl 


c 

FORM LCCAT-vICN CF SYMBOL 


PRl 

CC52 

13 

KL0C = .\PT( KSTYL)+KSUS-1 


PR 1 

0053 


IFIKSLB.EQ. iCOCCCC) KLOC=KLCC+I 


PRl 

COS'* 


KPL=KLOC 


PRl 

CC55 


IF {KPL.GE. 1000000 ) KPL=KPL-1000000 


PRI 


' C 

X-ELcMFNT TYPE 


PRi 

C056 


DC 14 KELTS=1, 5 


PRi 

CC57 


if(nsym{i,kpl).£g.n6lts{keltsi) go to 

5 

PRl 

:C59 

14 

CCNTINCS 


PR 1 

0C59 


K£LTS=0 


PRl 

CC6: 


GC TO 5 


PRl 


c 

STATS VARIABLES 


PRl 

0061 

IS 

KTYPS=1 


PRi 

00a2 


GO TO 13 


PRl 


c 

V2P, V3P, XIIP, OR XI5P 


PRl 

CC63 

16 

K“ XTD=KSTYL+3 


PRl 

0064 


KTYPE=0 


PRl 

CC65 


KLCC=KSUB 


PRl 


OAIA 

0A15 

0A16 

0417 

0418 

0419 
C42C 

0421 

0422 

0423 

0424 

0425 

0426 

0427 
042c 

0429 

0430 

0431 

0432 

0433 

0434 

0435 

0436 
04 37 

0438 

0439 

0440 

0441 
04 4 2 

0443 

0444 

0445 

0446 

0447 

0448 

0449 

0450 

0451 

0452 

0453 
04 54 

0455 

0456 

0457 

0458 

0459 

0460 
C461 
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0066 

CC67 

0068 

0069 

0070 

CC71 

0072 

CC73 

007^ 

0075 

0076 

0077 

0078 
0C79 

0030 

0C8I 


0082 

CC83 

oca^ 
CC £5 
0CS6 
0037 

008o 

0089 

CC5C 


0C5i 

0092 

0093 


0099 

0095 


date- = 78173 14/94/45 


PAGE 0003- 


0 

17 

C 

C 


0 


18 

C 


19 


C 

20 

0 

21 

0 

C 


0 

22 


C 

C 

C 


23 


GO TO 5 
"X"-TYPE 
K STAR =2 
"XM" TYPE 

IF(NAM?(KSTAR).cO.NM) go to 18 
"XI" TYPE’ 

IF(.'MME(KSTAR),N£.NI) GO TO 21 

KSTAR=3 

"XI#" TYPE 

IF(LPI,PP(NAP£,KSTAR,1,KSTYU ,EC,1 ) GO TO 5 

GO TO 16 

KLOC=KSUG 

"XM"-£L£P£NT TYPE 

DO 19 K£XT0=1,3 

IF(NA»1P(3).EQ.NXM(KEXTD)) GO TO 20 

CCATIAUe 

K£ XTD = 0 

KL0C=0 

CC TO 5 ' 

XM$I - TYPE 

IF:nAME(4),E0.NI.0R.NAHE( 5) .EC.NI) ’KTYPE=6 

GO TO *5 

X TYPEf NON SPECIAL - "XTABLE"? 

IF(LDCATE{NAMS tNXTAB,2 ,5) .eC.O) GO TO 22 

"XT" TYPE? tu 

IF(N'AP£(2}.EQ.ACHART> GO TO 23 

"XWStf" TYPE 

KSTAP=3 

KFUMP = LPUf'P(NAy£,XSTAR, 2, KSTYL ) 
KSTYL=KSTYL/iO 

J'^^KPb.^P.EQ.C.ANO.KSTYL.NE.O) GG TO 13 

KST=1 
GO TO 9 
"XTABLE" 

KELTS=8 
KTYPE=0 
GO TO 5 
ERRORS 

TABLE appeared IN XPRESSICN ILLEGALLY 

CALL ERROR (30J 

END 


PRl 0462 
PR1‘ 0453 
PRl 0964 
PRl 0965 
PRl 0466 
PRl 0467 
PRl 0468 
PRl 0469 
PRi 0470 
PRl 0471 
PR: 0472 
PRl C473 
PRl 0479 
PRl 0475 
PRl 0476 
PRl 0477 
PRl 0478 
PRl 0479 
PRl C43C 
PRl 0431 
PRl 0482 
PRl 0483 
PRL 0434 
PRl 0485 
PRl 0436 
PRl 0487 
PRi 0908 
PRl 0489 
PRl 0490 
'’Rl 0491 
PRl 0492 
PRl 0493 
PRI 0499 
PR: 0995 
PRl 0996 
PRl 0497 
PRl 045? 
PRl 0499 
PRi 0500 
PRl 0501 
PRi 0502 
PRL 0503 
PRl 0504 


Q O 


to 

JO 
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ERROR 
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PAGE OOCi 


OCOi 


SU6R0UTINE ERROR (NUM) . 

PRl 

0505 


C 


PRl 

C5C6 


C 

<><><><>00<;><><><><><>0<><><><><><><><>0<><><><><><><><> <><><>? Ri 

0507 


C 


PRl 

050« 


C 

PUT ERROR NUMBER 

PRl 

0509 


C 


PRl 

051C 


c 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRl 

0511 


c 


PRl 

C5L2 

CC02 


CCPMCN /FILES/ NPRINT,NPGM,NCTA,NEOU 

PRl 

C513 


c 


PRl 

0514 

0003 


DIMENSION N8AD(34) 

PRl 

0515 

0C04 


DATA NPACL/lfi/ 

PRl 

0516 

CCC5 


DATA N BAD/ 1,4, 5, 6, 7, 8, 17, 18,19,20,23, 24»27f 28, 30 f 3 1*32,34/ 

PRi 

0517 

CCCto 


WRITE <NPRINT,1) NUM 

PRl 

0516 

OCO/ 

* 

FORMAT I//,14H ERROR NUMBER ,12,20H SEE ERROR LIST ,//) 

PRl 

0519 

CCC8 


DC 2 1=1,N8ACL 

PRl 

052C 

CCC? 


IF CNUM.EO, N3A0 m ) GO TO 4 ‘ 

PRl 

0521 

oou 

2 

CONTINUE 

PRi 

0522 

OCll 


WRITE (NPRINT,3) 

PRl 

0523 

CC12 

3 

format (388 STATE VARIABLE PRODUCTION SUPPRESSED.) 

PRl 

C524 

CC13 


STOP.C 

PRi 

0525 

C014 

4 

WRITE (NPRINT,5) 

PRl 

0526 

C015 

5 

FORMAT (48H SEVERE ERROR - SENSITIVITY ANALYSIS SUPPRESSED.) 

PRl 

0527 

CC16 


■ STOP 8 • 

PRl 

0526 

0017 


END 

PRi 

0529 
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0001 


SUOROLTINE FINIE 

PRl 

0530 


• c 


PRl 

0531 


c 

<><><><><><><><><> <><><><><><><><><><>o<><>< 

XXXSOOOOOOOPR 1 

0532 


c 


PRl 

0533 


c 

FINISH LP THE STATE EQUATIONS PRODUCTION. 

PRl 

0534 


c 


PRl 

0535 


c 

OOOOOOOOOOOOOOOOOOOCXXXSOOOOOOOOOOPll 

0536 


c 


PRl 

0537 

0002 


COMMON /FILES/ NPR INT,NPGM, NDTA ,NEQU 

PRl 

0538 


c 

, 

PRl 

0539 

CCC3 


WRITE (NSOUil) 

PRl 

054C 

0004 ' 

1 

FORMAT ( lHC,79X,/,lHC,79N,/,lHCt35(2H<>) ,lH<r 

8X,/,lHC,79Xt/t9H PRl 

0541 



1 £N0,71X) 

PRl 

0542 

CCC5 


RETURN 

PRl 

0543 

0006 


END 

PRl 

0544 



cool 
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GtTEQS 
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PAGE OOOi 


0001 


CCC2 

CCC3 

OCOA 

CCC5 

0006 

0007 

OOOS 

0C09 

CCIC 

ecu 

0012 

CC13 

OOIA 

CC15 

0016 

CC17 

CC18 

0019 

0020 

0021 

0022 

CC23 

CC2A 

0025 

CC26 

0C27 

CC28 


0C29 

CC3C 


SUBROLTI.NE GSTfiOS 


C 

C 

C 

C 

c 

c 

c 

c 


ORI 

Pfti 


PRl 

GET EQUATIONS FROM SIMTR INTO (NPUSH) PCINTEC TO BY ' PRl 


(KLOCA) AND KLOC 


PRl 
PR 1 


<><><><><><><><><><><><><><><><><><><><>< XXXXXXXXXXX^OP*^ 

PRl 

UPLICIT PEAL#«lA-h,Q-Z) ■ 


C 0PM OK 
COMMON 
1 ) 

COMMON 

COMMON 

COMMON 


/PUSHED/ KLCCA(3 ,400) , KLOC, NPUSHdOOO ) tKPUSH,KPMAX PRi 

/NASA/ NSYM! 6, 300) ,KSYM,NRAY( 80) ,KRAY,KL£N(8) ,NPT (8) ,^4T (300 PRi 

PRl 

/CNTRL/ LIST,L1,CTR0L(160),KXM(3) ' • PRl 

/FILES/ NPRIMT,NPGM,KDTA,NEGU PRl 

/CHARS/ NS?,NLP,NRP,NC0M.,Nc0,.ND,N'PLUS,NMINUS,NSL,NZSR,NXMi(3PR1 

PR 1 
PRl 


l),NUM(10) 

DIMENSION NAME) 6), NAUGP(6), MCNBXI6) 
DATA NQ/lhO/,N'’/lHM/ 

DATA NC/IHC/ ,NP/1 HP/, NX/ lHX/,NI/i HI/ 
DATA NALGM/ih/,lriA,lHU,lHG,lHM ,1H// 
DATA NI0iV8X/iHI,lH0,IHN,lH ,1HX,IHM/ 


KSLBSV=0 

: GIT NEXT CARD 

1 NEOF=LREAD(NPGM.,0 ) 

2 IFlNRAYd ).EQ.NC) GO TO 1 

; CHECK FOR " = 

DO 3 KEQ=6, 50 

I F (NRAY (KEQ) .EC .NEC) GO TO 5 

3 CONTINUE 

IF(LOCATE( NRAY,NAUGM,13,6) .EC. C) GO TO 14 
IF(LCCATE(NRAY,NIONBX,13,6).EQ.O) GO TO 18 
: TEST "FNC" 

4 IF (NEOF.EQ. C) GO TO 1 

KL0C1=KLCC+1 
KLCCAd ,KL0C1 )=0 
KL0CA(2,KL0C1) =KPUSH 
RETURN 

C FCUMD 

5 IF(LIST.EC.l) WRITE (NPRINT,6) NRAY 

6 FORMAT (1X,6CA1) 

GO OM PARAMETER OR XPRESSION 

IF (NRAY (7 ) . EQ.N'P.OR. (NRAY (7 ) . EQ .NX. AND .NRAY( 8 ) .Nc.NI.ANO.NRAY( 8) .NPRi 
le.NM)) GC TC 7 PRl 

IGNORE "DO" LOOPS PR* 

1 F (NRAY(7 ) .EO.NC .ANO.NRAY( 8 ) .EO.NO) GO TO 11 PRl 

KRAY=7 RRi 


PRL 
PR 1 
PRl 
PRl 
PRl 
PR 1 
PRl 
PRl 
PRl 
PRl 
PRl 
PR 1 
PRl 
PRl 
PRl 
PRi 
PR 1 
PRl 
PRl 
PR 1 
PRi 
PRl 
RRl 
PRl 
PRl 


0545 

0546 

0547 

0548 

0549 

0550 
05 51 

0552 

0553 

0554 

0555 

0556 

0557 
0552 
0559 
0550 

0561 

0562 

0563 

0564 

0565 

0566 
0557 

0568 

0569 

0570 
05 7 1 
0572 
05 7 3 

0574 

0575 

0576 
05 77 

0578 

0579 

0580 

0531 

0532 

0533 
0584 
C585 
0536 
05 37 
0538 
053S 

0590 

0591 

0592 


- 5 . 0 - 
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CG3i 


KSUB=0 

PRl 

0593 


CC32 


Kt;Ql=KEO + l ' ■ , 

PR-1 

0594 



c 

IGNORE ZERO VALUES 

PR i 

0595 


0C33 


IF(.S'RAY(KEQ1).=0.NZER} GO TO 1 , 

PRl 

05 96 



c 

GET SYPDCL 

PRi 

0597 


CC3A 


lF(LPLLL(RRAY,KRAYtNAME,OjKSUE).eO.-l ) GO TO 12 

PRl 

0598 



c 

PUSH EQUATION 

PRi 

0599 


CC35 


KRAY = K>^C + i 

PRl 

06CC 


CC3£ 


IFlKSLDSV.NE.O.ANO.KSUB.EC.lOpOOOO) KSUB=KSU8SV 

■pRl 

0601 


0037 


KPSHSV=KPLSH 

PRl 

0602 


0C33 


CALL PUSH (NAMFfKSUB) 

PRl 

060 3 


CC3^ 


KSUBSV=0 

PRl 

C6C4 



c 

CHECK FOR NCTHlhG PUSHED 

PRl 

06C5 


OOAO 

• 

IFIKPSHSV.EO.KPUSH) KLOC=KLOC-1 

. P.Rl 

0606 



c 

GIT NEXT STUFF 

PRl 

0607 


CCAI 


GO TO 2 

PRl 

06C8 



c 

PARH. OR XPRESSION - SAVE MAKE IN SYPBOL TABLE 

PRi 

0609 


CCA2 

7 

KSTAR=‘' 

PRl 

C-6I0 


GC43 


KSY«=KSYy+l 

PRi 

C6n 


CC-?A 


IF (KSYY.GT. 300) GO TO 13 

PR 1 

0612 


00A3 


DO B 1=1,6 

PRl 

0613 


CC^6 

8 

NSYH( I,KSY.v )=NS P 

PRl 

0614 


CCA7 ■ 


00 9 1=1,6 , 

PRi 

0615 


00A3 


lF(.'i=.AY(KSTAR),EQ.NLP.DR.NRAVtKSTAR).EC,N£Q) CC TO 10 

PP,i 

0616 


OC^O 


NSYM( I.KSYH )=NRAY (KSTAR ) 

PRl 

0617 


CC3C 


KSTAR=KSTAR+1 

PRl 

0616 


0051 

9 

CONTI NL'E 

PRl 

0619 


0C52 

10 

K3AY=KP0+i 

PRl 

062'! 



c 

PUSH ECUATICN 

PRl 

0621 


CC53 


KPSH3V=KPLSH 

PR 1 

0622 


0054 


CALL PUSH ! NSYYI liKSYM) ,0) 

PRl 

0623 


0'C55 


IF IKPSHSV .rO.KPUSH) KL0C = KL0C-1 

PRl 

0624 


CC53 


GO TO 2 

PRi 

0625 



c 

DETERMNE IF ONE PA.S5 DO ICOP 

PRl 

0626 


0057 

u 

KSUBSV=NU(-iB(NRAY, 19,20,11) 

PRi 

0627 


G'05a • 


KS(JSV1=NUFSINR£Y,'2R»26, II ) 

PRl 

0628 


CCS? 


IFIKSLBSV.NE.KSBSVI.OR.KSUBSV.EQ,!) Ksuesv=0 

PRl 

0629 

> a 

0C60 


GO TO 1 

PRl 

0630 



c 

ERRORS 

PRl 

0631 • 



c 


PRl 

0632 



c 

bad symbol 

PRi 

0633 


0061 . 

12 

CALL ERROR (251 

PRl 

0634* 



C 

SYMBOL TAELE OVERFLOW 

PRl 

0635 


0062 

13 

CALL ERROR (261 

PRl 

0636 



c 

PARSE OUT XML, XMR, AND XMS LENGTHS 

PRl 

0637 


CC63 


IF(NRAY(31).NE.NCDK) GO TO 4 

PRl 

0638 

. 

CCtA 


DO l'5 1=1,3 

PRl 

0639 


0065 


IF(NRAY( 22) ,6Q;NXM( I) J GO TO 16 

PRl 

0640 
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C C 6 6 

■ 15 

CONTINUE 

0C67 


GO TO 4 

CC66 

16 

N£K=30 

CC69 

• 

NST=28 

0070 

17 

KXM( I )=NUMBINRAY,NST,N£N,II ) 

0C71 


IF(KXM(I).GT,7) CALL ERROR (35) 

CC72 


GO TO 4 

0073 

IS 

I F(NRAY( 2S) .NE.NCnf*) GO TC 4 

0C74 


CO 19 1=1,3 

CC75 


IF(NRAY(19).EQ.NXM(D) GO TO. 20 

CC76 

19 

continue 

0C77 


GO TO 4 

CC78 

20 

NEN=27 

CC79 


NST=25 

0030 


GO TO 17 

ocei 


Ei\C 


PRl 

PRl 

PRl 

PRl 

PRl 

PRi 

PRi 

PRi 

PRl 

PRl 

PRl 

PRl 

PRl 

PRl 

PRi 

PRl 


0641 

0642 

0643 

0644 

0645 

0646 

0647 

0648 

0649 
065C 

0651 

0652 

0653 
06 54 

0655 

0656 
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■OCOl 


CCC2 
. 0008 ' 


0004 

CCC5 

0CC6 

0007 

. CCC8 
0007 
0010 


ooi: 

0012 

0013 

CCIA 

ccis 

0016 


0017 

CCIO 

CCi? 

CC2C 

0021 

CC22 

CC22 

0024 


CC25 

0026 

0C27 

CC23 

CC29 

0030 


SUBROUTINE MIDUP 


C 

c 

c 

c 

c 

c 

c 


PRl 
PR 1- 


C 

c 


2 

3 


<><><><><><><><><><><><><><><> <> OOOOOOOOOOOOOOOOOP^j- 


HRITE n£GUG» M.'-TRICESt ANC HEACER 


PRl 
PRi 
PRl 

<><><><><><><> OOOOOOOOOOOOOOOOOOOO^X^^^^^^^^^^'^'^^ 

PRl 

IMPLICIT REAL*R U-HtO-Z) „ 

COMMON /HOLOIT/ HEAD! 9, 1 1) ,KHEAD t f^OUMl t 6LTSI3C0 ) i KELTS » K EL, SDR i 50 ) PRx 
1 , NS DR (6, 5 C ) ,KSrP,NDUH3, X ICC 50 ) , KSV , NOUM 4 , PAR ( ICC) ,NPAR< 6 ,1’00) ,KOP ,PR1 
2NCUM5,XTAE(500) ,NTAD(7,50 !,KTAB,KTAB2,K£J 
COMMON /CHARS/ N SP ,NLP ,NRP , RCC M , NEO , ND , NPLUS , NM INUS , NS L , N ZcR, NXH I 3PR x 
lltNUMdO) ' 

COMMON /FILES/ NPRINT,NPGM,NCTA,NEQU PRi 

/CNTRL/ LIST,L1,CTR0L(160) ,KXM(3) 

/NASA/ NSYM( 6,300) ,KSYM,NRAY( 80) ,KRAY,NLEN(8) ,NPT (3)>NT (300 PRl 

PRl 

/PUSHED/ KL0CA(3 ,400) ,KLDC,NPUSHCiOOO) ,KPUSH,KPMAX PRl 

EQUIVALENCE { TI ME ,C TROL( 1 J) • 

DATA NE/lh£/,NJ/lHJ/ 

PLT OcBLG 

IF(LIST.N£. 1) GO TO 19 

DEBUG facility PP^ 

WRITE (NPRINT,!) NPT,NLEN,KSYM 

FORMAT ( iSH POINTERS TO NAME S ,8 1 7 ,/ ,1 8 H LENGTHS OF NAMES ,8I7,/,7HPR1 


COMMON 
COMMON 
1 ) 
COMMON 


1 LIM!T=, 17) 

DC 2 I=1,KSVH 

WRITE (NPRINT, 3) I , ( NSYM ( J , I ) , J=l ,6 ) 
FORMAT ( IH ,I4,2H ,6Ai) 


4 

5 

6 
7 


8 

9 

10 


WRITE (NPRINT, 4) ( I ,NPUSH ( I > ,N PUSH { I ) , 1 = 1 ,K PUS H ) 

FORMAT (2C0(/, 1X,4 (219, 8X,A1,3X) ) ) 

WRITE (NPR INT , 5 ) ( I ,KLOCA ( 1 , I ) ,KLOCA( 2, I ) ,KLUCA( 3, I ) , 1= 1»KLQC) 
FORMAT { 2C0( /,! X ,3 (41 9,3X) ) ) 

WRITE (NPRINT, 6) KXM 

FORMAT (6F iiXML=,I3,7H «XMRG= , I 3, t>H #XMS=,r3) 

WRITE (NPRINT, 7) T I ME , KELTS , KS V , KCP, KS DR, KT AB 

FORMAT (12H OOOTI ME = , 1 PD 14. 7 , / , 7H KELTS=,I3,5H KSV=,I3,5H KOP=, 
113, 6h KSCR=,I3,6H KTAB=,I3,//) 

DO 10 1=1, KELTS 

I F(N SYM( 1 , 1 ) .NF.NE. AND. NSYMI 1 ,1 > . NE. NJ) GOTO 8 
•IF (NT (I ) .EQ.O) CO TO 10 
WRITc (NPRINT, 0) I.ELTSd) 

FORMAT (12H <><><>ELNO= , 13 ,7H ELVAL=,1 PD14.7) 

CONTINUE 


PRl 

PRl 

PRl 

PRl 

PRl 

PR'l 

PRi 

PRl 

PR I 

PRi 

PRl 

PRl 

PRl 

PRl 

PRl 

PRl 

PRl 

PRl 

PR 1 

PRl 

PRl 


0657 

065S 

0659 

U660 

0661 

0662 

Ob63 

C664 

0665 

0666 
0667 
0666 

0669 

0670 
06 71 
0672 
0573 

0674 

0675 

0676 
0577 
06 7 8 
0679 
0,630 
0681 
0662 
0633 
0684 
0635 
06 8 6 
0687 
06 6 8 

0689 

0690 

0691 
C692 
C693 
0694 
C695 
C696 

0697 

0698 

0699 

0700 

0701 
07C2 

0702 
0704 
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PAGE 0002 



C 

PUT DEFINED PARAMETER VALLES 

PRi 

0705 

crC3: 


IFIKOP.EO.O) GO TO 13 

PRl 

0706 

C032 


DC 11 1 = 1 ,KCP 

PRI 

0707 

CC23 

li 

WRITE (\PRUT,12) 1,PAR{I) 

PRl 

07CE 

0C3A 

12 

FORMAT ( 13H OOOD PNUM= 1 1 3 1 7H DPVAL=rlPD14.7) 

PRI 

0709 


'C 

PUT SOURCE CERIVATIVE VALUES 

PRl 

07 IC 

CC35 

X3 

IF(KSDR.eC.C) GC TO 16 

PRl 

0711 

0036 


DO 14 I=1,KS0R 

PRl 

0712 

0C37 

lA 

WRITE INPRINT, 15) I,SDR(I) ■ 

PRl 

0713 

CC38 

15 

format (14H <><><>SDRNUM=,I3,8H S DRV AL= , 1PD14 .7 ) 

PRl 

071A 


C 

PUT TABLE VALUES 

PRi 

0715 

0037 

16 

IF(KTAE.SC.C) GO TO 19 

PRl 

0716 

CC^O 


DC 17 I=1,KTA8 

PRl 

0717 

CCAl 

17 

WRITE {KPftIKT,18) I,XTAB(n 

PRl 

071S 

0CA2 

la 

FORMAT ( 14H <> <> <>TA3NUM= , I 3 , 8H TAB VA L= , 1 PD14 .7 ) 

PRi 

0719 


c 

PUT STATE variable EQUATION SEPARATOR 

PRl 

072C 

CCA3 

19 

WRITE (NE0U,20) 

PRl 

0721 

OOAA 

20 

FORMAT ( 1HC,79X,/,1HC ,35(2H<» ,31H< , 3X ,/ ,1HC ,79X ,/ ,1 HC,79X ) • 

PRl 

0722 

CCA3 


RETURN 

PRl 

0723 

CCA6 . 


END 

PRi 

0724 
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DATE = 78178 


14/44/A5 


PAGE 0001 


SUeROLTIN-e POP (KPCPQ) . PRl 

PRl 

<><><><><><><><><><><><><><><><><:><><><><><><><><><><><><><><><><> pRi 

PRl 

POP THE ECUATICA POINTED TO 8Y VALUE (KPOPO) PROP (NPUSH) PRl 

PR^ 

<><><><><><><>< ><><><><><><><><><><><><><><><><><><>0<><><><><>0PrI 

PRl 

implicit REAL>SE(A-H,0-Z) . • PRl 

COPMON /NASA/ NSYM6, 300 ) ,KSY,M ,NRAY( 80) ,KRAY,NLSN ( 8) ,NPT(8) ,NT(300PR1 
1) ■ PRl 

COMMON /PUSHED/ KLOCA(3t 4CO) ,KLCC,NPUSHtlOOO) fKPUSH,KPMAX PRl 

COMMON /SEGR/ NShUVI 500) ,KSHUVtKST,KENO ,KSUB,K£XTD,KPAR PRl 

CCPM.CN /FILES/ NPRINT,NPGM,NCTA,N£QU PRI 

COMMON /CHARS/ LSP,NLP,NftP,NCCM,N£Q,NC,NPLUS,NMINUS,NSL,N2ER,NXM( 3PR1 
DjNU.MIlC) PRl 

CCK-'CN /CNTRL/ LIST , L 1, CT ROL ( 160 ) t K XM ( 3 ) PRl 

DIMENSICiN NDVII3) PRi 

DATA NnVI/lHOTlHV,lHl/ PRl 


KSHUV=0 PRl 

KSU5N1=C PRl 

LOCATE THIS SY"OOL VALUE IN STACK PRl 

IFILOCTSTI ILOC,KPCPD) .NE.O) GO TO 20 PRi 

CHECK IF parens: around SUM TERMS REOUIREO PRl 

KPAR=KLOCA{ 3,IL0C) PRl 

IF(KPAR.E0.1.0R.NRAY(KRAY-1) .EO.NSL) CALL PUTCHR INLP) ■ PRl 

FCRM FCF LCCATERS PRl 

KST=KLOCA( 2,ILPC) PRl 

ILOCP=ILOC+i • PRl 

KEN0=KLCCA(2,ILCCP)-1 ■ PRi 

POINT TO VALUE - NOT NAME/CHAR PRl 

IF(MPUSH(KST).GT.C.ANO.NPUSHn<ST).LT.lCCCCGCOCO) GOTO? PRl 

NAME/CHAR FCUNC PRl 

IFUPAR..\6.C.G«.NPUSHIKST).NE.NPLUS) CALL- PUTCHR (NPUSHtKST)) PRl 

KST=KST + 1 • ■ "pRl 

END OF POPPING? ‘ ' PRl 

I F (KST. LE. KeNC) go to 3 ■ PRl 

YHS - °LT RIGHT FAREN IF LEFT CNE PUT PRl 

IFIKPAR.EC. X) CALL PUTCHR (NRP) PRl 

CHECK FOR MCRE ECUATIONS THIS ITEM PRl 

DO 6 ILCSV=2fKLCC PRl 

IFIKEND+l.EQ.KLQCAI 2,ILCSV) . ANO.KLOCAd ,ILCSV-1 ) .EO.KLOC'A (X , ILCSV ) PRl 
1 ) GO TO 19 . . • PRl 

CCNTINUc ORl 

ANYTHING left ON STACK? PRl 

IFILGRAEIKPO?) ,EO.C) GO TO 17 PRl 

NC - hE CCNE? PRi 
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POP 


DATE. = 73178 


lA/44/45 


t^Ai>t uuu<: 


0023 

OCE-? 

0030 

CC31 

0032 

0033 

003A 

0C3'5 

0036 

0037 
0036 

0039 

OCAO 

OOAl 

JG0A2 

0043 

OOAA 

00A5 

CCA6 

00A7 

0CA8 

0CA9 

CC5C 

0051 

0052 

0053 

CC5A 

CCS5 

0056 

0057 

0058 


IF IKST.LE.KSND) GO JO 3 
C YES 

Re TURK 

C SYMBOL VALUE GIVEN 

7 KPOP=NPUSH«KST ) 

C PUT PARENS. IF CIVIDc 

IF (NRAYI KRAY-lI . NE.NSL) GO TO 8 
KPAR=1 

CALL PUTCKR (NLP) 

C CHECK FOR INVERTED ARRAY 

e KMX=MOO(KPOP , ICCCCCO) /lOC'CO 

I F (KMX .£0 . 16, OR .KMX .00 . 26. OR .KMX .EO. 36) GO TO 16 
C FCSM CCRRECT VALUES FROM SUBSCRIPT 

9 IF (KSUB. NE. 0. AND. KPOP .GT. 1 OCOOOO) KP0P=KP0P+K§UB-1000001 

IF(KPGP/iCOCOO.NE , 0. AND. KPOP/ 10 0000. NE, 7) KP0F=KF0P+1 
KST=KST+1 

c PUSH THIS VALUE AND POINT TC NEXT ECUATICN 

IFILSHUVlKPQPt ILOC jKPOPO) ,E0. C) GO TO 12 

C GCT NEXT ECUATICN 

1C KPOS=KOD{KPCPC,ICOO) 

C LOCATION NON-EXISTENT - ZcRC VALUE THEN 

IF (KPCS.SQ.O) CO TO 15 
KTYP = M0C( K PC PC r 1000000 )/ 10000 
C EXTENDED SYMBOL AT BAD LOCATICN 

IF{MOC{KPCPO,lCCCCGO)/lOuCCO.NE.O) GO TO 15 . 

C NCT ENOLGh RCCM Tu PUT SYMBOL - PUT THE LINE 

IF(KRAY->6.GT. 7’) CALL PUTLIN 
C PUT "V/ I" IF NECESSARY 

IF (KTYP. LE.O ) GC TO 11 
CALL PUTCHR (NOVI (KTYP)) 

C PUT NAME 

11 CALL FUTN.iP (NSYMd.KPOS ) ) 

GO TO A 

C NAME FOUND - PUT ELEMENT 

12 IF(KPOP.LT.ACCC ) GO TO 10 
IF(K?CP/1COO.EC.9) GOTO 10 
IFdLCC.SC. C) GO TO A 

C FORCE SGUATION IF NOT EXTENDED SYMBOL 

IF(KEXTC.EO.D) CD TO 1 

IF (KSHUV.EO. G) GC TO 21 

C LOCATE THIS EQUATION SUBSCRIPT IN STACK 

CC 13 IS=1,KSHUV,5 
I T=KSHUV-I S-3 

IF(NSHUV( IT).NE.O) GO TO lA 

13 CONTINUE 

C ■ SUBSCRIPT NCT IN STACK 
GO TO 21 

C GET SUBSCRIPT 


PRl 0773 
PRl C77A 
PRl 0775 
PRl 0776 
PRl 0777 
PRl 0778 
PRl 0779 
PRl 0730 
PRl 0781 
PRi 0782 
PRl 0783 
PRl 078A 
PRl 0785 
PRi 0736 
PRl 0737 
PRl 0786 
PR: 0739 
PRi 0790 
PRl C791 
PRl 0792 
PRl 0793 
PRl 079 A 
PRl 0795 
PRi 0796 
PRl 0797 
PRi 0793 
PRl 0799 
PRl 08 j 0 
PRl C8C1 
PRl 03C2 
PRl 0603 
PRl cacA 
PRl 08CS 
PRL 0806 
PRl 0807 
PRl C8Ce 
PRL 0809 
PRl 0610 
PRl 0311 
PRi 0612 
PRl 0813 
PRl C81A 
PRl 0815 
PRi 0816 
PRl 0817 
PRl 0816 
PRl 0319 
PRl 0820 
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POP 


DATE = 78X78 


CCS9 

. 14 

OCbO 



C 

CC61 

15 

0062 



C 


C 

0063 

16 

0C6<. 



C 

CC65 


CC66 



C 

CC67 


006S 


0Co9 


0C7C’ 



c 

C07X 

17 

Cl72 



C 

0073 


0C7/- 


CC73 


007b 


0077 


CC78 



C 

0079 


CCBC 



c 

oos: 

la 

0032 


CCS3 


00 34 


0C35 

1'9 

GCS6 


CC87 


0083 



C 

ccas 

20 


C 

0C90 


CC91 



C 


c 


c 

0C92 

21 


KPCP=N»USH( KST-1)+NSHUV(I 
GO TO 1 


T) 


PUT ZERO value 
CALL PUTCHR INZER) 

GO TO 4 

START ARRAY MULTIPLY LOOP 
I.SVALin SUBSCRIPT C,N INVERTED- GUY 
IFIKSUR.Nl.GT.l) CO TO 9 
IFIKPCP,.LT,XCCCCCC) GO TO 2X 
SAVE CLC SU8SCS AND GEN NEW ONE 
KSU3N=KSUfi 
KSU9NX=X 

PUT INVERTED ARRAY NAME 
KMXA=KMX/XO 


KSUS=KSt8NX 

CALL PUTMX (KMXA,KSUBN,KSUB ) 

6C TO 9 

CONTINUE MATRIX MULTIPLY 
KHXB=KPOP/XCCCC 

I F ( KMXE. NE. 16 . AND 4 KMX B.NE .26 . AND .KMXB .NE. 36 ) GO TO 4 

T-JCREMENT SUBSCRIPT 2 

1F(KSUBNI.cO.KXM(KHXA> ) GC TC 18 

KSU8N1 = K.SU0NX + 1 

KSUB=f:SLDiM 

KST=NSHLV(KSHUV-3) 

NSHUV {KSHUV-4 )=KSUD 
KFCPC=KPCP 

PUT NEW INV. ARRAY NAME 
CALL PUTMX IKMXA,KSUBNfKSUB) 

GC TO 4 

END OF MATRIX MFY 
X SUB=KSUBN 


NShUV (XSHUV-3 )=NSHUV( KSHUV-2)+-l 
KSUBN1=0 
GO TO 5 
ILOC= ILCSV 

I F (KPAR.2C. X ) KRAY=KRAY-1 
KPAR=KLCCA ( 3,1 LCC) 

GO TO 2 

CMSCK FOR INVALID STATE VARIABLE 
I F IKPCPCVlCCOOO-.LT.X .OR. KFOPC/ 100000 

SIMPLY A ZERO VALUE ARRAY ELEMENT 
CALL PUTCHR INZER) 

RETURN 
E RROjR S 


6T.3) GO TO 22 


INVALID SUBSCRIPT 
CALL ERROR I3X) 


X4/44/45 page 0003 


PRi 

. 0821 


PRl 

0822 


■PRi 

0823 


PRi 

0624 


PRi 

0825 


PRI 

0826 


PRI 

0827 


PRi 

0328 


PRI 

0829 


PRL 

C83C 


PRi 

0831 


PRl 

0832 


PRl 

0833 


PRl 

0834 


PRi 

0835 


PRl 

C336 


PRl 

0837 


PRl 

0838 


PRl 

08 39 


PR 1 

C84C 


PRl 

0841 


PRl 

0842 


PRl 

03 4 3 


PRl 

0844 


PRl 

0845 


PRL 

0846 


PRi 

0647 


PRi 

0848 


PRl 

0849 


PRl 

C85C 


PRl 

0351 


PRl 

0352 

Q 0 

PRl 

0853 


PRl 

0854 


PRi 

PRl 

PRl 

0855 

0856 

0857 


PRi 

0858 


PRl 

0859 

ss ■ 

PRl 

U860 


PRl 

0861 

S a 

PRl 

PRl 

0862 

0863 


PRl 

CB6A 


PRl 

0865 


PRl 

0866 


PRi 

0867 


PR 1 

0868 
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PRi 

0869 


■ C 

ikvalid name to 

POP 



PR 1 

C87C 

CC93 

22 

CALL uRRCR (32) 




PRI 

0871 

CC9A 


END 
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PUSH 


DATS =. 78178 


IA/A4/45 


PAGE 0001 


CCCl 


0C02 

CCC3 

CCCA 

CCC5 

0006 

CCC7 

OCOT 

0C09 


CCIC 

cell 

0012 

0013 

COIA 

0015 


CClo 

0017 

ccie 

ocn 


CCiO 


CC21 


0022 


CC23 

OOCA 

0025 

0026 



SUOROLTINE PUSH (NAPfc,KSUfi) 


PR 1 

0872 

c 



PRl 

0673 

c 

<><><><><> OOOOOOOOOOOOOOOOOO.OOOOOOOOOOPHl 

0874 

c 



PRl 

0375 

c 

PUSH THE EQUATION POINTED TO QY (NAME) 

INTO (NPUSH) 

PRl 

0876 

c 



PRl 

0877 

c 

<><><><><><><><><><><><><><><><><><><><><><><><><><><><>0<><>0<>PR1 

0878 

c 



PRl 

0870 


COPMON /NASA/ N S YM ( 6 , 300 ) , KS YM , NRA Y( 80 ) 

,KKAY,NLEN«8) 

,NPT(8) ,NT(300PRL 

0660 


1) 


PRl 

08 31 


COPMON /PUSHED/ KLCCA(3,AC0) ,KLCC,NPUSH{1000J,KPUSH, 

KPMAX PRl 

0832 


COPMON /FILES/ NPRI NT, NPGH,NOTA ,NEQU 


PRl 

088 3 


CCPMCN /CHARS/ NS P , NLP, NRP , NCOH , N EQ, ND, 

NPLUS,NMINUS, 

NSL,NZER,NXM(3PR1 

C884 


1) ,NUH( 1C) 

< 

PRi 

08S5 


DIMENSION NAM£(6), NAMK6) 


PRi 

0836 

c 



PRl 

0837 


KC0=0 


PRl 

0888 

c 

DECODE NAME - POINT TO SYN3CL TABLE 


PRl 

0669 


CALL CECNAM ( N AME , K SUB , KLOC 1 ,KVAL ) 


PRl 

0893 

c 

e.AC NAPE - I.E. NOT IN SYPBOL TABLE 


PRl 

C891 


IF (KVAL.EC. C) GO TO A 


PRi 

0892 

c 

PUSH VALUE equation 


PRl 

0893 


KLCC=KLCC+1 


PRl 

0894 


KL0CAI1,KL0C)=KVAL 


PRl 

0895 


KLOCAI 2,KL0C)=KPUSH 


PRl 

0896 


KL0CA(3,KL0C)=O 


PRl 

0897 

1 

■ KCO=l 


P.R 1 

0898 


K SHSV = KPUSH 


PRl 

0899 

c 

GET ECUATION SYMBOL 


PRl 

0900 


IF{LPULL(NRAY,KRAY,NAfH,KC0,!<SU31 ) .SQ 

1) RETURN 

PRl 

09 Cl 

c 

FORM SUM CF PRODUCTS FLAG 


PRl 

09C2 


IFIKOO.E0.2.ANO. (KGO.EO. 1 .OR .NP USH ( K SHS V) . EO. NPINUS) 1 KLCCA {3,KL0CPRl 

U9 03 


1 )=1 


PRi 

C9C4 


KG0=1 


PRi 

0905 

c 

DECODE EOLATION SYMBOL 


PRl 

0906 


CALL CECNAM <NAH1,KSU31,KL0C1tKVAL) 


PRl 

0907 

c 

INVALID SYMBOL 


PRl 

C9C8 


IF (KVAL.EQ. C) GO TO 5 


.PRi 

0909 

c 

CHECK TO PUSH FULL NAME 


PRl 

0910 


IF {KVA L. EG.6 j'30 . OR .KV AL . £0 .7000 .OR.KVAL 

.EQ.8U00) GO 

TO 2 PRl 

C911 

c 

'•X=XH" FORMAT - IGNORE 


■ PRl 

0912 


IF(KVAL.EQ.KLOCA(l,KLOCn GC TC 1 


PRl 

0913 

c 

PUSH ECUATION SYMBOL VALUE ONLY 


PRl 

0914 


NPLSH(KPUSH)=KVAL 


PRl 

0915 


KPUSri=KPUSH+l 


PRl 

0916 


IF(KPUSh.GT.KPMAX) GO TO 6 


PRl 

0917 


GC TO 1 


PRl 

0918 

c 

PUSH FULL NAME 


PRl 

0919 
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CC27 

2 

K5=l 

PRi 

C92C 

CC2E 

. 3 

NPUSH('<PUSH)=r{A('UKS) 

PRl 

0921 

0C2'^ 


KPUSH=KPUSH+1 

PRI 

0922 

CC3: 


KS=KS+l 

PRl 

0923 

CC3i 


IFCKPLSH.GT.KPPAX) GO TO 6 

PRI 

092A 


C 

END OP NAHH? 

PRi 

0925 

CC32 


IF {NAy-KKS ) .HO .NSP ) GO TO 1 

PRl 

092 6 



if(ks.gt.61 gc to 1 

PRl 

0927 

C0 3A 


GO TO 3 

PRl 

0928 


C 

ERRORS 

PRl 

0929 


C 


PRl 

093C 


C 

PUSHED F:APE NOT IN SYMBOL TABLE 

PRi 

0931 

0035 

4 

CAUL ERROR (27) 

PRi 

0932 


C 

ECUATICN NAME NOT IN SYMBOL TABLE 

PRl 

0933 

CC36 

5 

CALL ERROR (28) 

PRi 

C93A 


C 

stack TVESFLOh 

PRl 

0935 

0C3"^ 

6 

CALL cRRUR (29) 

PRl 

0936 

CC3' 

» 

5 NO 

PRl 

0937 
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PAGE 0001 


CCCl 


C 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE PUTCFR INCHAR) • . PRi OSBS 

PP.l 0«5 5G 

<><><5<>0<>0<><><><><><><>0<><><><><><><><><><><><><><X><><><>PP'i "0940 

PRl 0941 

PLACE THE CHARACTER (NCHAR) INTO THE OUTPUT ARRAY (NRAY) PRl 0942 

IP (NRAV) LONGER THAN 72, PUT (NRAY) CN FT02FC01 AS FORTRAN PRi 0942 

PRl 0944 

OOOOOOOOOOOOOOOOOOOO <><><> OOOOOOOOOOPRl C945 

PRl 0946 


0002 

C 

common /NASA/ NSYM{ 6,300), 
1 ) 

0003 


IF (KRAY. G7. 72) CALL PUTLIN 

GC04 


■ NP.AYIKRAY ) = NCHAR 

CCC5 


KRAY=KRAY+1 

CCC6 


RETURN 

0007 


ENO 


PRL 
PRl 
PRl 
PRl 
PR 1 
PRl 
PRl 


0947 

0948 
C949 
09 30 
0951 
CS52 
0933 
095* 
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PUTlt 


DATE = 78178 


14/AA/45 


PAGE 0001 


0001 


SUSilOtTINE PUTTT ( KGO , KSTCP , NCHA R , KLEN' , KWT, KP T ) 

PRl 

0955 


c 


PRl 

0956 


c 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRl 

0957 


c . 


PRl 

0958 


c 

PUT THE NAMES FROM ELEMENT LIST STARTING AT TYPE 

IKGC) PRl 

0959 


c 

TC TYPE (KSTOPl IN OUTPUT, PRECEDE BY (NCHARl IF 

PRl 

0960 


c 

KLEN N<^ ZcRC ANC SET NUMBER ELTS IN (KNT) 

PRl 

0961 


c 

AND POINTER OF FIRST NAME IN (KPT) 

PRl 

0962 


c 


PRl 

0963 


c 

<><•><><><><><><><><><><><><><><><><:><><>< ><>< 30000000000P R i 

0964 


c 


PRl 

0965 

0002 


COMMON /CHARS/ NSP , NLP ,NRP , NCCM , NE G , NO , NPLUS , NMI NUS , NS L , NZER , NX M (3PP.1 

0966 



1),NUM(10) 

PRl 

0967 

CCC3 


COMMON /NASA/ NSYM (6,300 ) ,KSYM, NRAY (801 iKRAY, NL EN ( 8 ) , NP T ( 8 ) , NT ( 300PR 1 

0968 



i) 

PRl 

0969 

OCO-t 


COf^MON /HOLCIT/ HEAD('9,li),KHEAO,NOUML,ELTS(300) , 

KELTS, KEL,SOR(50)PR1 

0970 



1 ,NSDR(6,50) ,KSeR,NCUM3,XICt50),KSV,N0UM4,PAR{ 100) 

tNPAR(6,1001,KDP,PR1 

,0971 



2NOUH5,XTA3( 500) ,NTAB(7,50) , KT AB , KTAB2 , K E J 

PRl 

0972 

0C05 


DIMENSION MRAY(6) 

PRl 

0973 


c 


PRl 

0574 


c 

ZERO POINTER AND COUNT 

PRl 

0975 

OGOb 


KPT=0 

PRl 

0976 

0C07 


KNT=0 

PRl 

0-9 77 


c 

LCC? THRU TYPES , 

PR 1 

0578 

CCC^i 


DO A IPT=KG0,KSTCP 

PRl 

0979 


c 

FIND LENGTH OF THIS TYPE 

PRl 

09 30 

CCC4 


MPT=NL*=N( IPT) 

PRl 

0981 


c 

ZERO LENGTH IS NC PUTSKIS 

PR 1 

05 3 2 

00 iO 


IF(MIPT.EO.C) GO TO 4 

PRl 

09 3 3 


c 

FIND THIS GUYS POINTER LESS ONE 

PRl 

093A 

cell 


LFT=N?T(1PT)-1 

PRl 

0935 


c 

LOO? ON THAT CCLNT 

PRl 

0936 

0012 


CC 3 JPT=i,MPT 

PRl 

0987 


c 

PUT THE -EXTRA CFAR IF OME 

PRl 

09E8 

CC13 


NI';PT=0 

PRl 

0905 

0‘j . A 


IF(KL£N.N£. 1) GO TO 1 

PRl 

0990 

C 0 1 3 


KRAY (i )=NCHAR 

PRl 

0991 

cc: 6 


NNPT = 1 

PRl 

C992 


c 

FORM THIS GUYS ABSOLUTE POINTER 

PRl 

0992 

0017 

1 

ljpt=lpt+jpt- 

.PRl 

0994 


c 

BUILD first pointer 

PRl 

0995 

ccis 


IF(KPT.EO.O) KPT=LJPT 

PRl 

0596 


c 

• PUT THIS NAME 

PRl 

0997 

ocn 


MS=6-N.NPT 

PRl 

0998 

CC2: 


DC 2 M = 1',('S 

PRl 

0999 

CC2: 

2 

MRAY( M+NNPT) =NSYK(M,LJPT) 

PRl 

ICCC 

CC22 


CALL PUTNAM (MRAY) 

PRl 

1031 


c 

UPDATE COUNT 

PRl 

1CC2 
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FORTRAN IV Gl RELEASE 2.0 


pum 


CC23 KNT=KNT+1 

• C PUT A COPPA OETWIXT EM 

CC2A CALL PUTCHR (NCCM) 

OC25 3 C0NTI.>JUE 

0C26 4 CONTINUE 

CC27 . RETURN 

0C22 END 


D'ATE =! 78178 


14/44/45 


PAGE 0002 


PRi 1003 
PR-1 1C04 
PRI 1CC5 
PRi 1006 
PRI 1007 
PRI icca 
PRi ICCS 
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FORTRAN IV 
OCCl 


0Cii2 

0003 

CCCA 

000:5 

CC06 

CCC7 


CCC8 

CCC9 

0 C 1 0 
CCil 
0012 
0013 
CCIA 
CC15 

0016 

CC17 

0013 

OCIO 

CC2-: 

0C2I 

0022 


G: release 2.0 PUTLIN date = 78178 IA/AA/A5 .PAGE 0001 


C 

C 

C 

C 

C 

C 

c 

c 


c 

c 

c 


l' 


c 

2 

3 

A 


SIBROLTINE PUTLIN PRX 

PRl 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOppi 


WRITE THE CHARACTER STRING INRAY) IN FORTRAN FORMAT 
AFTER .SLANKING END OF LINE.' CONTINUATIONS ADDED 




PRl 

PRl 

PRl 


0<><><><><><><>00<>01:><><><><><><><><><><><><>'<><><><><><><><>PR1 

PRl 

CCPMON /NASA/ NSYM (6 ,300') , KSYM,NR AY ( 80 ) , KRAYt NLENI 8 ) »NPT( 8 J , NT( 3CCPR1 
. PRl 

CUMKON /FILES/ NPRrN.T,NPGM,NDTA ,N£QU PRl 

COMMCN /CHARS/ NS P , NLP , MRP , NCOM, N EQ, NO , NPLUS, NM INUS ,NSL ,N2£R , NXM( 3PR 1 
1) ,NUMilO> PRX 

dimension NC0MM(6) prx 

•DATA NCOMM/IHC, IHO, IHM, IHM, IHO, IHN/ PRi 

DATA NAST/IHS/ ■ PRi 


NOTHING ON LINE NOW 
IFIKRAY.FO.?) RETURN 
NRAY(o) =NSP 

IF "COMMON " OR «> = " APPEARS, THEN NO CCNTINUATICN 
IF(L0CAT = (NRAY,NC0MM,,7,6).EQ.C) GO .TO 2 
KND=KI.VC(20,KRAY) * • ■ 

00 1 I=a,KND 

IF(MRAYm.EQ.NEO) GO TO 2 

CCNTINUE 

NRAY«6>=NAST 

BLANK TO END OF LINE AND PUT LINE 

DC 3 I=KRAY,73 

NRAYm=NSP 

WRITE INcQU.AI I NRA Y( I ) , I =1 ,72 1 
FORMAT t72Al,8X) ' 

KRAY=7 

RETURN 

END 


PRl 
PRl 
PRl 
PR 1 
PRl 
PRl 
PRl 
PRI 
PRl 
PRl 
PRl 
PRl 
PRi 
PR 1 
PRl 
PRl 
PRl 
PRl 
PRi 


lOlC 

1011 

1012 

1013 

lOlA 

1015 

1016 

1017 

1018 
ICIS 
IC2C 
1021 
1C22 
1023 
102A 
1025 
iC26 

1027 

1028 
1C2S 
1C3C 

1031 

1032 
1032 
103A 
1035 
1C36 
1037 
1039 
lo39 
ICAC 
lOAl 
1062 
1C43 
ICAA 
1045 
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FORTRAN IV G1 RELEASE 2.0 


PUTVX 


DATE = 78178 


1A/A4/A5 


PAGE OCO 


0001 


CCC2 

CC03 

CCC4 


0005 

0C06 

CCC7 

0003 

OCC'J 

CCIO 

cci: 

0012 

•0013 


SUDROLTINt PUTPX ( KMX,XSUBN,KSL'BN1 ) PRl 

C PRL 

C OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRi 

C PRl 

C PUT • WHERE % IS 'L't'R'j OR 'S' ANO PRl 

C ##1 IS KSU8N, ANO hUZ IS KSUBNl . PRl 

C ' , , P^l 

C <><><><><><><><><><><><><><> <><><><><><><><><><><><><><><><><><><>PR1 

C PRl 

CCPRON /CHARS/ NS P , NLP, NRP , NCOM , N £Q » NOf NPIUS, NM iNUSt NSL » N2ER ,NXM( 3PR 1 
DfNUMIlC) PRl 

DI-'^cNSION NAMX(6) PRl 

DATA NAPX/1HX,1HM, Ih ,LHI,1HI,1H / PRl 

C PRl 

NAMX( 3)=NXM{KHX) PRl 

IF(KSUeNl.GT.l) CALL PUTCHR INPLUS) PRl 

CALL PUTNAP (NAPX) PRl 

CALL UCDOLT (KSLON) ■ PRl 

CALL PUTCHR (NCQM ) PRl 

CALL eCCCOT (K5USN1) PRl 

CALL PUTCHR (NRP) PRl 

RETURN PRl 

END . • , PRl 




a 


v’es- 


1046 

1047 
1C48 
iC4<; 
1030 
1051 
1C52 

1053 

1054 
IC55 
1056 
1357 
IC58 

1055 
1060 
1061 
1062 

1063 

1064 
1C65 
1066 

1067 

1068 



FORTRAN IV GI RcLEASS 2.0 


PUTNAM 


DATE = 78178 


1V44/A5 


PAGE 0001 


0001 


SUBROUTINE PUTNAM (NAME) PRl 1069 
C PRi 1C7C 
C <><><><><><><><><><><><><><><><><><><><><>0<><><>0<><><><><><><>PRi 10 7! 
C . ' ' ■ PRl 1072 
C PUT ThP (UP TO) SIX CHARACTER LIST (NAME) IN THE OUTPUT PRl 1073 
C STRING (NRAY). V»RITE (NRAY) IF OVER 72 CHARS. LONG PRl 1C7A 
C PRl 1075 
C <><><><>0<><><><><><><><><><5<><><><><><><>0<><><><><><><><>0<>PR1 1076 
C . ' ■ PRl 1077 


CCC2 


COMMON /NASA/ NS YM( 6 ,300) . KSYM» NRAY (30 ) ,KRAY ,NLEN(8 ) 
1) 

,NPT (0),NT (300PR1 
PRl 

1C7E 

1079 

CCC3 


COMKCN /CHARS/ NS P, NLP, NRP,NCOM,NEQ,ND,NPLUS,NHINUS, 
1) ,NUK(LO) 

NSL,NZER,NXM(3PRi 

PRl 

1C8C 

1031 

OOOA 

C 

DIMENSION NAHE(6), NEW(6) 

PRl 

PRl 

1082 

1033 


C 

RE-FORM MECHANICAL NAME 

PRl 

1C8A 

CCC3 


CALL 3ACK (NAM£,NfcW) 

PRl 

1C05 

0006 


LEN=i 

PRl 

1086 

0007 

' 1 

IF(KRAY.GT.72) CALL PUTLIN 

PRl 

1087 

CCC5 


NRAY(KPAY1=NEW(LEN) 

PRl 

lOSE 

OCO'J 


KRAY=KRAY<-1 

PRl 

1089 


c 

GNC OF LINE? 

PRl 

1099 

GCIO 


LEN=LCN+1' 

PRl 

1C91 

cell 


IF(LEN.GT.6) RETURN 

PRl 

1092 - 

0C12 


IFINEWILEN) .EQ.NSP) RETURN 

PRl 

1U93 

CC13 


GC TO 1 

PRl 

109A 

CC lA 


• END 

PRl 

1095 
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FOBTI^A^ 

CCCl 


0002 

CCt3 

OOOA 

CCC5 


CCC6 

CC07 


0CC5 
CCC9 
0013 
COIL 
CC 12 
0013 
CCi-!* 
CCi 5 
CCli 
CCL7 
CCL6 
CC19 
0020 
0C21 
CG22 
002 3 
CC24 
CC25 
0026 
0027 
CC25 


CC29 


IV G1 BEL6ASE 2.0 SETUP 

StUROLTIN'E SETUP 


date = 78178 


1A/4A/45 


PAGE 0001 


C 

C 

c 

c 

c 

c 

c 


c 

c 

c 


PRl 1096 

. PRl 1097 

- - PRl lice 


G=T ELEPEKT NAKES AND POINTERS INTO (NSYM , INPT 1 , ( NLEN)., AND 
FRQh“ SCEPTRE ROUTINE "SIMUL8" 

OD1 1 1 ns 


PR 

PRl 


1131 

11C2 


PRl 

IMPLICIT REAI,’!'n(A-H,0-Z) ' p^, 

^rrvwnv /r\TRL/ L 1ST , L 1 , T I ME, S 1 ( 2 1 ) , XNOPRQ , S 2 , XNHEAO , S 3 ( 6 ) .XRUNNO.SPRl 

^COPMON-'/KOLCIT/ HE A0( 9 , 1 1 ) , KHS AD , NOUMl ,ELp ( 3C0 ) .KELTS , KcL , SDR ( 50 ) PRl 

l.NSDRt6,5G),KSCR,N0UH3,XIC(5O),KSV,NDUM4,PARUCO).NPAR(6,10O),KDP,PR. 


2N0UM5, XTAD (SCO) ,NTAB( 7 ,50) , KTAB ,KTA02 ,K EJ 


CINENSICN 
DIMENSICN 
D (MEN SION 
DlPSNStCN 
DIPENSICN 
0[ MEN SION 


, NMATB(6) 

, NUERIV(6) 
NCUM(6 ) 


NBA0(13), NELS(8), NL0C(8) * 

OLM(300), NMUT(6,50) 

NSIMTRt 5) , NSRCDR(6), NTABLE(6) 

NXMUTL(51, NNAMES(5)» NPRMTR(6) 

N'I0NPX(5)t NCUTPTC6), NNDPD(3)t 
CTROUIAO) 

ECUIVALENCE ( CTROL ( 1 ) , T IMS ) 

DATA NBAC/34,45,55,71,91,99, 101,102, 111,121,121, 141, 151/ 

i\£LS/96,97 ,9 8,93 ,9^,,95 ,92 ,100/ 

NLQC/7,4,5,6, 1,2,3, 8/ 

NSIMTR/IHS, IHI, IHM, It-T, IPR/ 

NNAMd S/IHN ,lHA,lKM,lhE,lHS/ 
NPRMTR/1HP,1HR,1HM,1HT,1HR,1HS/ . 

ND5RIV/lHC,lhE,’lHR, IHI, IHV.IHS/ 

MCN'BX/lHl ,1HC,1HN,1H , IHX/ 

N0LTPT/1H0,1HU,1HT , IHP , IHT ,1HS/ 

NAST/1H«/ 

NNCPC/IH, ,1I-D,1HP/ 

NSRCDR/iHS vlUR ,1NC ,1H0 ,1HR ,1 HV/ 

NTABLE/IH.T, IHA, 1H8,1HL, IHS.IHS/ 
NXNUTL/lHX,lHM,lHU,ihT, IPU/ 


DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 


GET DATA 
READ TILL 
K.^UT=0 


SUBROUTINE SIMTR FOUND 



PRl 

PRl 

PRl 

PR 1 

PRl 

PRl 

PRl 

PRl 

PRl 

PR 1 

PRl 

PRl 

PRl 

PRi 

PRl 

PRl 

pr: 

PRi 

PRl 

PRl 

PRl 

PPl 

PRl 

PRl 

PRl 

PRl 

PRl 


1105 

11C6 

1107 

1138 

11G9 

me 

1111 

1112 

1113 

1114 

1115 

1116 

1117 

1118 
U19 
112C 
1121 
1122 

1123 

1124 

1125 

1126 
112 7 
1128 
1129 
113C 

1131 

1132 
1123 

1134 

1135 

1136 

1137 

1138 

1139 
114C- 

1141 

1142 

1143 


- 67 - 



FGRTSA^ IV G1 RSLEASE 2.0 


SETLiP 


OATH => 78178 


0C30 

CC21 

0032 

CC33 

CC3A 

OC35 

CC36 

0037 
0C33 
CC3'3 
OCAO 
OCA 1 

CCA2 
OOA3 
c:aa 
CCA 5 
0GA6 
• CCA7 
CCA ri 
OOA9 
0C50 
CC 5 1 

0052 

CCS? 

CC5A 

0055 

0C56 

0057 

0053 

CC5S 

* CC60 
CC=1 
0062 
C063 
CC6A 
CC65 

00 6 ^5 

’CC67 

ccta 

OCoO 

CC70 


1 ■ READ (NPCM, 2,ENC=39 ) NRAV ' • 

2 FORMAT '(80A1) . , 

I F ( locate I NRA Y»NSIMTR *9 i5 ) . EC. 15 GO TO 1 

C CHECK WHICH SCEPTRE IS USED 

READ (NPGF.S) NRAY 
3’ FORMAT ( ////////>80A11 
MAC=70 

1 F (NRAY 16 5 .5Q.NAST ) NAC=160 

c read ccntrols • 

READ ( MDTA r A »EMD=38) { C TRC L (I ) » I -1 i NAC ) 

4 ' format (1P5C1A.7) 

IF (NAC. £0.160) CC TC 6 
DO 5 (=71,160 

5 CFROL ( I )=C .CDC 

C LCCATE certain STRINGS 

6 IF (LREADtNPGM.C) .60.1) GO TG AO 

7 IF(LOCAT6(NRAY,(JSRCOR, 15, 6) .EQ.C) GO TO 20 
IF (LOCATE INRAY ,NT ABLE, 15, 6 ) . EQ .0 ) GO TO 22 
IF ( LOCATE (NRAY , NXPUTL,15 ,5) . EC.O I GO TO 2A 
IrCNAC.60.160) Gn TO 8 

’ IF ILOC ATE (NRAY ,N"i AMES , 15, 5) . EO. 0 ) GO TO 46 
I F (L0CATE(NRAY , APRMTR, 15 ,6 ) . EC.O ) GO TO 28 
IF (LOCAT£(NRAY ,ND£R I V ,15 ,6) . EC.O) GO TG 30 
I F (LOCATE (N'5AY,N lUNBX, 13, 5) . = 0.0) GO TO 32 

8 I F (LCCATE (NRAY,ACUTPT,15 ,6 ) . EC .0 ) CO TO 9 
GO- TO 6 

C CHC-CK FOR INVALID FLAGS IN CONTROLS. 

9 Oc'lO TER=1,13 
KLOC=NBAO (I cR) 

' ■ IF(CTROL{KLOC) .NE.C.COO) GO TC 36 

10 CONTINUE 

C FCRM CFEUG FLAG 

LI ST=XWROBG+XWRTSM 
IF(LIST.GT.C) LIST=1 . 

C "EXECUTE SETUP PHASE ONLY MUST BE FLAGGED 

IF { XE XSTP. N£. IDC) GC TO 37 
C FORMA ELEMENT TYPE POINTERS AND LENGTHS 

NPT(7 )=l 
DC 11 T=l,8 
MPT=NLOCI I ) 

KL0C=‘J=LS(MPT ) 

NLEN(XPT)=CTRQL(KL0C) 

Ir(H?7.£Q.7) GC TO 11 
[•iPTM = NLOC( I-l) 

NPTIMPT )=NPT (MPTM )+NLEN(MPT-M) 

11 CONTINUE 

DO 12 1=1,8 

IF (NLEN( I) .EQ.^ ) NPT(I)=0 
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1A/AA/A5 


PAGE 0002 


PRl llAA 
PRl 11A5 
PRl 1146 
PRl 1147 
PRl 1,148 
PRl 1149 
PRl 1150 
PRl 1151 
PRl 1152 
PRl 1153 
PRl 1154 
PRl 1155 
PRl 1156 
PRl 1157 
PRl 1158 
PRi 1159 
PRl 11^>P 
PRi 1161 
PRl 1162 
PRl 1163 
PRi 1164 
PR: 1165 
PRl 1166 
PRl llo7 
DRl 1168 
PRl 1169 
PRi 1170 
PRl 1171 
PRl 1172 
PRl 1173 
PRl 117^ 
PRi 1175 
PRl 1176 
PRl 1177 
PRl 1178 
PRl 1179 
PRi 1180 
PRi 1131 
PRl 1182 
PRl 1183 
PRl 1184 
PRi 11-85 
PRl 1186 
PRl 1137 
PRl 1138 
PRl 1189 
PP.l 1190 
PRl 1191 



FQRT^A'l 

0071 

CC72 

0073 

007A 

0C75 

0C7o 

CJ77 

0C78 

0079 

ocso 

0031 

0082 

0C83 

QC84 

OCSS 

ocao • 

oce? 

0060 

CC69 

0090 

CC7- 

CC92 

CC93 

0C9A 

CCSD 

0C96 
• CC97 
CC98 


CC99 

010 0 
Cl^I 
C102 

0103 

0104 

cic5 
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12 

CONTINUE 

?R1 

1192 

C 

GIT HEADING CARDS 

PRl 

1193 


KH£AO=XNHEAC 

PRl 

1194 


I F(KH“40.N5.0) READ tNDTA,l3) HEAD 

PRl 

1195 

13 

format {9A0) 

PRl 

1196 

C 

GIT ELEiVENT VALUES AND NAK.ES ’ ' 

PRl 

1197 


KELTS=XN“LTS 

PRl 

119S 


READ (UDTA.i4) ( t L TSd ) i I = 1 t KEL TS ) 

PRl 

1199 

14 

FORNAT (5C14.7) 

PRl 

12CC 


.IFINAC.EO. ISO! BEAD (NDTA.15) KNSV M II t J ) . I=i i 6 ), J = 1,KELTS ) 

PRl 

1201 

15 

FORMAT (72A11 

PRl 

120? 


KSY.M=K=LTS 

PRl 

12-J3 

C 

SKIP CVER ELEMEBT IC VALUES 

PR 1 

12C4 


READ IN0TA,14) t C UKl I) , I =1 , KE LT S ) 

PRl 

1205 


READ (NCTA,14) ( DUM ( I ) , 1= 1 , K£ L TS ) 

PRl 

1206 


I FIKNUT.NE.O) read (N0TA,14) ICUM(I), I=1,KMUT) 

PR 1 

1207 

C 

GIT SCORCc DERV VALLES 

PRl 

12C8 


IFIKSDR.NS.O READ (N0TA.14) ( SDR ( I) , I =I , KSOR ) 

PRi 

1209 

C 

SKIP CVER VMCT- AND' INICT VALUES 

PRl 

121C 


^•5AVE=^JLrM:)•^^L£^^3)•^•NLE^(4)+^LgNI6 ) 

PRl 

1211 


IF(NSAVE.NE.O) READ (M0TA,14) ( DUM II ) t I =i * NSAV E ) 

PRl 

1212 

C 

SKIP CVER C6F PARM DERIVS VALUES 

PRl 

1213 


NCPD=XNDPD 

PR 1 

1214 


IFINOPO.Mc.C) READ (N0TA,i4J <0U« f I) » I =1 , NDPO) 

PRl 

1215 

C 

GIT I3D0T AMD V4D0T VALUES 

PRl 

1216 


KSV = \LEM3)+NLEM4) 

PRl 

1217 


IF (KSV.EQ. 0) GO TO 41 

PRl 

1218 


K2=MLE-'I (3) 

PRl • 

1219 


IF(K2.NE.O) READ (NDTA,14) tX 1C( I > , I=1,K2 » 

PRl 

122C 


■K 1=K2<-1 

PRl 

1221 


K2=NLEN( 4 )+K2 

PRl 

1222 


IFIK2 .NE.NLEN(3n READ (iNCTA,14) ( X IC [ H , I==K1 ,K2 5 

PRl 

1223 

C 

GIT DP NAPES AND VALUES 

PRl 

1224 


KDP-=XNOP 

PRl 

1225 


IFIKCP.NE.O) R = AD IN0TA,14) ( PAR I I) , I = 1,KDP) 

PRl 

1226 


IFIKOP.NS.O.ANC.NAC.EQ.160) READ (NDTA,15J i iNPARIIi J )i 1=1* 6),J=] 

LtPRl 

1227 


IKDP) 

PRl 

1229 

C 

GIT TAGLES 

PRl 

1229 


lF(KTA(!.h-S.O> READ <N-DTA,14) ( XTABI I) , 1= 1 *KTAB ) 

PRi 

123C 

c 

COPY PLTUAL NAPES TC NSYM 

PRi 

1231 


00 17 1=1,KMUT 

PRl 

1232 


KSYM=KSYP+1 

PRl 

•1223 


DO 16 J=l,6 

PRl 

1234 

16 

NSYP.I J,KSYH)=NPLT( J,n 

PRl 

1235 

17 

CCNT INUE 

PRl 

1236 

C 

PUSH EOUATICNS— >STACK 

PRl 

1237 


CALL GETSOS 

PRl 

1238 

C 

FORK cLEHENT TYPE ARRAY 

PRl 

1239 
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01C6 


DO 19 T=1»KELTS 


PRl 

126C 

0iC7 


ip{EL rs: n .Nt. c.cDO) go to is 


PRl 

1261 

0 IC s 

. 

CAUL CECNAN ( NSYM { 1 , I) , 0 , I XUOC, IV AU 


PRl 

1262 

CK9 


NT(1)=1-LCCTST(IXL0C,IVAL) 


. PR 1 

1263 

OliO 


GO TO 15 


PRl 

12 66 

cm 

18 

NT(I)=1 


PRl 

1265 

C112 

19 

CCNTIMJE • 


PRl 

1266 

0U3 


return 


PRl 

1267 


C 

source derivatives exist - GIT THEM 


PRl 

1268 

o::a 

20 

KRAY =22 


PRl 

1245 

0115 


1<SDR = J 


PRl 

125C 

01X6 

21 

IF(LPULL(nRAY,KRAY,NSDR( 1,KS6 r+1) ,0,KSB) . £0.-1 ) 

GO TO 7 

PRl 

1251 

0117 


KSCR=KSCR+1 


PRl 

1252 

oils 


GO TO 21 


PRl 

1253 


C . 

TABLES £XIST~GIT THEM TOO 


PRl 

1256 

CXi'V 

22 

KHAY=22 


PRl 

1255 

C12; 


KTAB=J 


PRl 

1256 

0121 


KTA32=0 


PRl 

1257 

0122 

23 

IFaPULL(.YRAY,KRAY,NTA8(l,KTAB+ll , 0 , NTAB ( 7 , KTAB+1 ) ) . E Q. -1 ) 

GO TO 7PR1 

1258 

0123 


KTAQ=KTAE + NTAfi t7,KTAB2) 


PRl 

1235 

0126 


KTAB2=KTAB2+1 


PRl 

126C 

0^25 


GO TO 23 


PRl 

1261 

• 

C 

C2T MUTUALS 


PRl 

12b2 

0126 

26 . 

KRAY=21 

, 

. PRi 

1263 

01.27 

25 

IF(LPLLL(NRAY,KRAY,NSYMa ,KMUT+l > ,0 ,KS8) .cQ.-l ) 

GO .TO 7 

PRi 

1266 

0128 


KMUT=KMUT+l 


PRl 

12*5 

C129 


GC TO 25 


PRl 

1266 


C 

GET ELEMENT NAMES AND NUMBERS 


PRl 

1267 

0130 

26 

KRAY=21 


PRl 

1263 

0131 

27 

IF (LPULL (NRAY,KRAY,NSYM< 1,KELTS+1) ,0,KSB) .EQ.-l) 

GO TO 7 

PRl 

1269 

01-2 


K£LTS=K£LTS+1 


PRl 

127C 

0133 


XNcLTS=K£LTS 


PRl 

1271 

0136 


GO TO 27 


PRl 

1272 


C 

GGT NUMEER AND NAMES OF DPS 


PRl 

1273 

0135 

26 

KRAY=22 


PRl 

1276 

0136 

29 

IFfLPULL(NRAY,KRAY,NPAR(l,KOP+l),0,KSB).EQ.-i) GO TO 7 

PRl 

1275 

0137 


KCP=KCP+1 


PRl 

1276. 

012S 


NXDP=KDP 


PRl 

1277 

0130 


GO TD 25 


PRl 

1278 


C 

GET NUMBER OF CEFINEO PARAMETER DERIVS 


PRl 

1279 

0160 

30 

DC 31 1=22,70 


PRl 

128C 

0161 


IF ( LOCATE (NRAY,NNDPD,I ,3) .EC.O) XNDPD=XNCPD+1 


PRl 

1281 

0162 

31 

ccntinug 


PRl 

1282 

0163 


GC TO 6 


PRl 

1263 


C 

GST NUMBER CF ELEMENT TYPES 


PRl 

1286 

0166 

32 

IF(MPAY{ 18) ,EQ.NRAy{9>) GO TO 6 


PRl 

1285 

- 0165 


KRAY = 17 


PRl 

1286 

0166 

33 

IF(LPLLL(NRAY,KRAY,NDUM,0,KSB).EQ.-1 ) GO TO 7 


PRl 

1287 
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FORTR;l^ IV Gi RcLEASE 2.0 SL-TUP DATE 


C147 


PC 34 1=2,10 

0148 


IMNDLyi 3).E0,^U^im) GO TO 35 

0149 

' 34 

CONTl.vUc 

0150 


1=2 

0151 

35 

CTR0U\ELS(I-1) )=KSB 

0152 


GO TO 33 


C 

c 

ERRORS 


Vr 

c 

INVALID FLAG ERRORS 

0153 

36 

CALL ERROR (1£P) 


C 

EXECUT": SETUP PHASE ONLY NOT ’SPEC IF lED 

0154 

37 

CALL EPROR 114) 


C 

NO CONTROLS! NO DTASAV) 

Oi55 

38 

CALL ERROR (17) 


C 

NC PRCGRAP SI.VTR(N0 PGNSAV) 

0156 

39 

CALL ERROR (IS) ‘ 


C 

NO TRANSIENT ECUATIONS IN SIMTR 

0157 

40 

call ERROR (19) 


C 

NO STATE VARS 

OiSfi 

41 

CALL ERROR (20) 

Oi'59 


END 
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PR: 12S5 
PRi 1290 
PKi 1291 
P^l 1292 
PRi 1293 
PRI 1294 
PRI 1295 
PRi 1296 
PRI 1297 
PRI 1298 
PRI 1299 
PRI 1300 
PRI 1301 
PRI 13C2 
P.Rl 13C3 
PRI 1304 
PRI 13C5 
PRI 1306 
PRI 1307 
PRI 1308 



FORTRAN IV GI RcLEASE 2.0 


LGRAB 


OATS' = 78178 


1 ^/ 44/45 


PAGc OOOi 


ooo: 


FUNCTION LGRA8 (KPDP) 

PRl 

1309 


c 


PR i 

1.31C 


c 

<><><><><><><><><><><><><><><> <><><><><><><><><><><><><><><><><><>PR1 

1311 


c 


PRl 

1312 


c 

PULL THE PAST LOCATORS FROM (KSHUV) 

PRl 

1313 


c 


PRl 

1314 


c 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRl 

1315 


c 


PRl 

1316 

CCC2 


lyPLICIT REAL=«8 U-H.0-2) 

PR 1 

1317 

0003 


COMMON /SHGR/ N SHUVI 500) , KSHUV , KST ,KEND ,KSUB, KEXTD,KPAR 

PRl 

1318 

0C04 


COf'MO.N /CUTRL/ L I ST , L 1, CTkOL ( 160 ) , KXM ( 3) 

PRl 

1319 

CCC3 


CCPMCN /FILES/ NPRINT,NPGM,NCTA,N£QU 

PRl 

132C 


c 


PRl 

1321 

0006 


IFILIST.EC. 1) WRITE (NPRINT.l) 

PRl 

1322 

CCC7 

1 

FORMAT (9H LGRAEESD) 

PRl 

1323 

cccs 


LGRAB=1 

PRl 

1324 

0000 


IF(KSHUV.EQ.O) RETURN 

PRl 

1325 

OCLO 


LCRAB=0 

PRl 

1326 

CCli 


KPAR=i\SHUVI KSHUV) 

PRl 

1327 

0Q12 


KSHUV=KSMLV-1 

PRl 

1328 

0013 


KPOP=NSHUV IKSHUV) 

PRl 

1329 

CC14 


KSHUV=KSHUV-1 

PRl 

133C 

CCi5 


Kc ND = NSHLV(KSHUV) 

PR i 

1331 

OOU 


KSHUV=KSHUV“1 

PRl 

1332 

CCL7 


KST=N5PUV (KSHUV ) 

PRl 

1333 

CClS 


KSHUV=KSHLV-1 

PRl 

1334 

0010 


KSUQ=NSHUV( KSHUV) 

PRl 

1335 

OCiO 


KShUV=KSKUV-l 

PRl 

1336 

0C21 


RETURN 

PRL 

1337 

0022 


END 

PRl 

1338 
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PAGE 0001 


CCCl 


CC02 

0003 
CC04 
C Vs 0 3 
CCC6 
0007 

cccs 

CCC4 

0010 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c , 


1 


FLMCTION LOCATE ( NR AY , NS T R , KSTK, KLEN ) - PRi 

PRi 

<><><><><><><><><><><><><><> OOOOOOOOOOOOOOOOOOOPRl 

PR 1 

LOCATc STRING (NSTR! IN ARRAY (NRAY) WITH PRI 

LENGTH (KLSN) STARTING AT (KSTR). PRI 

LOCATE=C IF string FOUND . PRI 

LOCATc=l IF STRING NOT FOUND PRI 

PRi 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRi 

PRI 

DIMENSION NRAYiaO), NSTRI80) PRi 

PRI 

LDCATE=1 PRi 

CO 1 '1=1,KLEN PRi 

KSTRIM=KSTR+I-1 PRI 

IF(NRAYtKSTRlM).N£.NSTR!I)) RETURN PRI 

CONTINUE PRi 

LaCATE=0 • . PRI 

RETURN PRI 

END PRI 


133<; 

1340 

1341 

1342 
134 3 

1344 

1345 

1346 

1347 
1345 

1349 

1350 

1351 

1352 

1353 

1354 
13 5 5 

1356 

1357 

1358 

1359 


o o 

^ to 


►rj O 
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LCCTST 


DATE = ,78178 


1A/4A/45 


PACE CCOi 


ooo: 


FUNCTION LOCTST (IL0C,IP0P1 

PRI 

136C 


c 


PRl 

1361 


C 

OOOOOOOOOOOOOOOO <><><><><><>< XXXXXX’XX^OOP^l 

1362 


c 


PRl 

1363 


C, 

DcTcRMTME IF TH5 LOCATOR (IPCF) IS A VALIC SYP80L 

PRl 

1364 


C 

AMO FINC ITS PniNTFR (ILOC) IN (KLOCA) AND SET L0CTST=0 

PRl 

1365 


c 

IF INVALIC SET LCCTST=1. 

PRl 

1366 


c 


PRl 

1367 


c 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRI 

1368 


c 


PRl 

1369 

CCC2 


IMPLICIT REALMS (A-H ,0-Z) 

PRl 

137C 

0003 


COMMON /NASA/ NSY«( 6,3CC) ,KSYP ,NRAY( 80 ,KRAY,ALEN(8 ) ,NPT 18 ) 

,NT (300PR1 

13 71 



1 ) 

PRl 

1372 

COCA 


. CCMHON /PLSHEC/ K LOCA (3 *400 ) , KLOC f NPUSH I iOOO > t KPUSH, KPMAX 

PRl 

1373 

0C05 


COMMON /SHGR/ N£HUV( 500 ) ,KSHUV ,KST,KEKD,KSua,K£XTD,KPAR 

PRl 

1374 

OOOb 


COMMON /CNTRL/ LIST, LUCTROLI 160) ,KXM(3) 

PRl 

1375 

CCC7 


CCMMGN /FILES/ NPRINT ,NPGM, ,\CTA,N EQU 

PRl 

1376 

CCC8 


DIMENSICN MPTI6) 

PRl 

1377 

ccco 


DATA MPT/7,4,5,6,1,2,3,8/ 

PRl 

1378. 


c 


PRl 

1379 


c 

FORM PCINTERS 

PRl 

138C 

0010 


KcXTD=C 

PRl 

1381 

ecu' 


LOCTST=0 

PRl 

1382 

C012 


ILCC=.J 

PRl 

1333_ 

CC13 


HPOP=MQDtIPCP,lCCCCCO) 

PR 1 

1384 

0014 


HLOC=MOD( IPOP, ICCC) 

PRl 

1385 


c 

NC LCCATICN GIVEN 

PRi 

1386 

CC 1 5 


IF(MI.GC,EO.C) GC TO 5 

PRl 

1387 

0016 


MMIO=MPOP-MLUC 

PRl 

1383 


c 

CFiECK FOR INVERTED ARRAYS 

PRl 

1339 

CC17 


MMIDA = f'PIC/1000G 

PRi 

L39C 

CC13 


IF{MMIDA.c0.i6.0R.MMIDA,£C,26,CR.MMI0A.EQ.36) GO TO 8 

PRl 

1391 

0C19 


DO 4 rLOC=l,KLOC 

PRl 

1392 


c 

SEARCH STACK FCR EQUATION POINTERS 

PR 1 

1393 

0C2C 


JPOP=MCD ( KLCCA (1 , I LCC) .lOCOOCO) 

PRl 

1394 

0021 


JL0C = M0D( JPOP 7 ICCO) 

pr: 

1395 


c 

ALL RUT SUBSCRIPT HATCH 

PRl 

1396 

CC22 


IF (MPOP.EQ. JPCR) GO TO 6 

PRl 

13 97 


c 

cxtendfd and type match only 

PRl 

1398 

0023 


IF(MMin.NE.(JPOP-JLOC) ) GO TO 4 

PRl 

1399 


c 

NCT EXTENCEC 

PRl 

14CC 

CC24 


IF( JPO°. LT. ICCCCO) GO TO 1 

PRl 

14C1 


c 

NO POSSIBLE MATCH - TRY AGAIN 

PP.l 

1402 

CC25 


IFIKLCCAIl, iLOr I.LT.IOOOOOO) GO TO 4 

PRl 

14C3 


c 

FOUND A MATCH 

PRl 

1404 

0026 


KEXTD=JPO?/lCCCCO 

PRl 

1405 

0027 


GO TO 7 

PRl 

1406 


c 

DETERMINE SUBSCRIPT POSSIBLE MATCH 

PRi 

1407 
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LOCTSr 


DATE = 78178 


14/4A/A5 


PAGE 0002 


0C28 

CC29 

CC30 
CCBl 
C3 3i 
CS33 
CC2A 


CC 3c 
CC37 
003 3 
CC39 
CC^C 

CO-Vl 
CC A2 

OOA3 

COAA 
00-!-5 
C0-.6 
C C A 7 

OOAS 

CCAS 


1 IFIIPOP.LT.lOOOOOO.AND.KUOCAtIjILOO.LT.iCCCOCC) CO TO 4 - PRl 

IFIMLOC.LT. JLQC) GO TO 4 PRl 

C FI.VO SUBSCRIPT POSSIBLE RANGE PRl 

: CO a J=X, 8 PRl 

JI'PT = RPT(J) PRl 

IFIJLCC.lt. NPT(JMPT)) GO TO 3 PRl 

2 CGNTIWUE PRl 

J=9 PRl 

3 J=J-1 PRl 

C CHECK SUBSCRIPT WITHIN RANGE PRl 

IFU.L'^.O ) go tc 4 PRl 

jpPT=f'OTlJ) PRl 

• IF(NL£N( JMPD.EO.O) GO TO 3 PRl 

IF( (f-'-LOC-JLOC) .LE.NLENI JMPT)“1! 60 TO 6 PRl 

4 CCNTIMJE PRl 

C NO POSSIBLE MA^CH PRl 

5 L0CTST=1 PRl 

GC TO 7 . PRl 

C FORM SUBSCRIPT FROM GIVEN RANGE PRl 

6 IF! IPJP.GT. ICCCCCO.OR.KLOCAI 1,1 LOC) .GT. ICCOOOO.OR.KEXTO.NE.O ) KSUBPRl 

1=PL0C-JL0C+1 PRl 

7 IF(IP'J?.LT.4000»AND.L0CTST.£C.0) GO TO 10 PRl 

8 IFILIST.SO.iJ URITE (NPRINT,9) LOCTST PRi 

9 FORMAT (13H LOCTST VALUE, 12) PRl 

RETURN PRl 

C INVALID VARIABLE ELEMENT ' PRl 

10 CALL ERROR (34) PRl 

END PRl 


1403 

14CS 

1410 

1411 

1412 
14x3 

1414 

1415 

1416 

1417 

1413 
1415 
1420 
l'42l 

1422 

1423 

1424 

1425 

1426 

1427 

1428 
1425 

1430 

1431 

1432 

1433 

1434 

1435 
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LPULL 


DATE = 78178 


14/A4/45 


PAGE 0001 


OCOl 


0002 

CCC3 

CCCA 

DGli 

OCC^ 

CCC7 

0CC8 


OCC-? 

CCIC 

ecu 

CC12 

0012 

OCiA 

CC15 

0016 

C'C17 

OOiS 

CC19 

CC2C 

0021 

CC22 

0023 

CC2A 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


2 


c 

3 

C 


. A 
C 

c 


FUNCTION LP'ULL (NRAY, KSTR, NSYM,KOO,K SUB) . PRl 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPf^l 

PRi 

PULL NEXT SY.XCCL FROM INRAY) BEGINNING AT (KSTR) PRi 

UP- TO COLUMN 72 INTC (NSYM LEFT JUSTIFIED. READ PRI 

NEXT CARD IF END OF LINS FOUND AND CONTINUE PRi 

NEXT CARD -MUST EE CONTINUATION CARO, ELSE LPULL. NE.O PRi 

UPDATE (KSTR) TC NEXT LOCATION. P<ii 

LPULL=C GOT SYMBOL P^i 

LPULL=1 GOT SYMBOL AND GOT NEXT CARD PRi 

LPULL=-1 DID NOT GET SYPBCL OUT GOT NEXT CARO PRi 

PRI 


I r ^ i 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOP^i 

PRI 


IMPLICIT REAL*e(A-H ,0-ZI 

COMMON /PUSFED/ KLOCA(3,AOO),KLOC,NPUSH(1COO) tKPUSHiKPMA 
CCMMO.N /CNTRLX L I ST , LI , CT ROL ( 160 ) , KXM ( 3 ) 

COMMON /FILES/ NPR I NT , NPG M , NOTA , (lE QU , 

COMMON /CFARS/ NSP ,NLP ,NRP ,NCOM ,NEO ,ND ,NPLUS,NMI NUS,NSL , 
1) ,NUM.(10) 

OIMENSICN NRAY(£0) , NSYM(6), NALPHA(36), NXTAS(S) 

DATA NALPHA/iHA, IHO, IHC, IHO, 1HE,LHF,1HG,1HH,1HI 
1 NflHOilF.P ,1HC, IFR, IFS , IHT t IHU, IFV, IHW, IhX, IHY , IHZ, IHO, IH 
21MA,lriS,lH6,lH7,lH8,lH9/ 
data NXTAB/IHT ,1MA, IHB, IHL, IHE/ 

DATA NPER/IF./ 

EQUIVALENCE ( NN ,NALPHA (lA 1 ) , (NX ,NALPHA 12A ) ) 


LPULL=0 

FIND BEGINNING (ALPHA) CHAR OF SYMBOL 
IF(NRAY(KSTR ) .EQ.NSP) GO TOS- 
CO 2 1=1, 26 

I F (NRAY(KSTR) . EC.NAL?HA( I ) ) GC TO 8 

CONTINUE 

NCT A SYMBOL 

IF(KOC.EQ.O) GC TO A 

PUSH ncn-symbol character 
IF(NRAY(KSTR).FQ.NSP) go to a 

N°USH (KPUSH )=NRAY (KSTR) 

KPUSH=KPUSH-»-l 

SET KOO TO 2 IF SUM OF PRODUCTS TERM 

IF (NR AY (KSTR ) .<=O.NPLUS .OR .NRA Y ( KSTR ) . EQ. NMINUS ) 

IF(KPLSH.GT.KPMAX) GO TO 27 

KSTR = KS7RF1 

TRY TC FIND TH= SYMBOL AGAIN 
IF(KSTR.L£.72) GO TO 1 
END OF RECORD-GET ANOTHER 


KD0=2 


PRI 
PR i 
PRI 

NZcR,NXM{3PRI 

PRi 

PRI 

1HL,1HM,1HPR1 

1, 1H2, 1H3,PRI 
PRI 
PRI 
PRI 
PRI 
PRI 
PRI 
PR 1 
PRI 
PRI 
RRl 
PRI 
p'r: 
PRI 
PRI 
PRI 
PRI 
PR 1 
PRi 
PRi 
PRi 
PRI 
PRI 
PRI 
PRI 


1436 

1A37 

1A38 

1439 

lAAC 

iAAl 

1442 

1443 

1444 

1445 

1446 

1447 

1448 

1449 
143C 

1451 

1452 

1452 

1454 

1455 
14 56 
1457 

1453 
1459 
146C 

1461 

1462 

1463 

1464 

1465 

1466 

1467 

1468 

1469 

1470 

1471 

1472 

1473 

1474 

1475 

1476 

1477 

1478 

1479 

1480 
1431 
1482 
1482 
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CG2S 

CC26 

CC27 

0C2B 

OG29 

0C3C 


0031 

0C32 

0033 

CC3A 

0035 

0036 

0037 


CC36 
0C29 
00 AO 

CCAl 

00A2 

CCA3 

CCAA 

OCAS 

CCA6 

0CA7 
C C A 3 


OCAO 

C 


CC5I 

CC52 

CC53 

CC5A 

0055 

CC56 


5 


C 


C 

6 

7 

C 

3 

C 

9 
C 

C 

10 


C 


C 

11 

12 

c 

c 

13 


C 


C 

lA 

C 

15 


NECF=LR£AD(NPGN,0) 

FORMAT llX.eCAl ) 

CHECK CO.’'!TINUATION CHARACTER 
IFtNRAY(6 ).E(3.NSP) GO TO 6 
IFlLIST.EC.l) WRITE (NPRINTt5) NRAY 
FIVO SYMBOL IN NEW RECORD 


KSTR=7 
GO TO 1 
SYMBOL NOT 
LPULL=-1 
KSTR=7 • 
RETURN 
IF SY.MPOL 
IFd.EC.A. 


IN THIS RECORD-NC CCNTINUATICN 


IS «CCC" IGNORE IT 

AND.NRAY (KSTR-1 l.EQ.NPER) -GO TO 


CLEAR SYMBOL 
00 9 1=1,6 
NSYM( n = NSP 
GET SYMBOL 


00 lA 1=1,6 

IGNORE SPACES IN SYMBOL 
1F(NP.AY(KSTR1 .NE.NSP) GO TO 11 
K£TR=KSTRtl 

1FIKSTR.LE.72) GO TO 10 

END OF RECORD - FAKE END- OF SYMBOL 

KSTR=KSTR~1 

GO TO 15 

CHECK ALPhANUMlFRIC CHAR. HERE 
DC 12 J=l ,36 

IF INRAY(KSTR) .EC. NALPHAU) ) GC TO 13 


CONTINUE 
END OF SYMHCL 

GO TO lA ^ 

GOT AN ALPHANUM CHAR. - SAVE IT _ 
NSYM{i)=NRAY(KRTR) 

KSTR=KSTR-H 

CHECK FOR END OF RECORD 
IF(KSTR.LE.72) CO TO lA 
END OF R-ECORD - GET NEW RECQRC 
ME0F=LREAD (NPGM ,C) 

’check FCR CONTINUATION FLAG 
IF (NRAY (6) .EG.NSP 1 GO TO 25 
IF { LI ST. EG. 1) WRITE (NPRIWT,5) NRAY 
OK CONTINUE RECORD-FIND SYMBOL 
KSTR=7 
CONTINUE 

CHECK FOR SUBSCRIPT 
IF(NRAYIKSTR).NE.NLP) RETURN 
KSTR=KSTR-h1 


FCUNC 


3 


PRl 

1A8A 

PRi 

1A3S 

PRl 

1A06 

PRl 

1A87 

PRl 

1A8F 

PRl 

1A89 

PRl 

1A9C 

PR 1 

lAO: 

PRl 

1492 

PRl 

IA93 

PRl 

1A9A 

PRl 

IA95 

PRl 

1A96 

PRl 

IA97 

PRl 

1A98 

PRl 

1A99 

PRl 

1500 

PRl 

1501 

PRl 

1502 

PRl 

1503 

PRl 

15CA 

PRl 

1505 

PRi 

1506 

PRl 

15 07 

PR 1 

iscs' 

PRi 

1509 

PRl 

1510 

PRl 

1511 

PR 1 

1512 


PRl 

1513 


PRl 

151A 


PRl 

1515 


PRl 

1516 


PRl 

PRl 

1517 

1518 


PRl 

1519 - 

.O ni 

PRl 

1520 


PRl 

1521 


PRl 

1522 


PRl 

PRl 

1523 

152A 


PRl 

1525 


PRl 

1526 


PRl 

1527 


PRl 

1528 


PRl 

1529 


PRl 

1530 


PRl 

1531 
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LPULL 
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PAGE 0C03 



C 

FTiKE AN "N" SUBSCRIPT 

PRi 

1532 

CC57 


KSUB=1CCCCC0 

PRl 

1533 

0053 


IF(NRAYIKSTR) .cs.NN} GO TC 19 

PRi 

1534 


C 

IF NOT "N", GET SUBSCRIPT NUMBER 

PRl 

1535 

CC59 


IF{LPUPP2(NRAYt!<STR,KV1).EQ.1) GO TO 22 

PR 1 

1536 

CC65 


IF:nR4Y(KSTR>.EC.NCCM) go TC 16 

PRl 

1537 

0061 


KSU8=KV1 

PRl 

1538 

C C c2 


CC TO 19 

PRl 

1539 

CC53 

16 

DO 17 11=1,3 

PRl 

154C 

0069 


If (NSYMl 3) .EQ.N/MI I I ) ) GO TO 18 

PRi 

1541 

00 3 5 

17 

CCNTINUE 

PRl 

1542 

CC66 


GC TO 22 

PRl 

1543 

CC67 

18 

IF(KXP(II).EO.C) GO TO 22 

PRl 

154 4 

OOoS 


KSTR=KSTR+1 

PRl 

1545 

CC69 


IF(LPURP2(NBAY,KSTR,KV2).£Q.1> GO TO 22 

PRl 

1546 

CC7: 


KSU8=KXR(II)'S’<KV2-i)+KVl 

PRl 

1547 

0071 

19 

RSTR = KSTR + i 

PRl 

1543 

CC7' 


IF{i<SUO.NE.100''000) GO TO 21 

PRi 

1 549 

CC75 

20 

IF(N'>jY(KSTR-1).£G.NRP) GC TC 21 

PRl 

155C 

007i 


KST.R=KS7R + 1 

PR 1 

1551 

OC7n 


GC TO 20 

PRi 

1552 

Cl7<3 

21 

RETURN 

PRI 

1553 


C 

INVALID SUBSCRIPT NUMBER-CHECK EX*PRESSION 

PRl 

1554 

0077 

22 

IF(N$YM( ll.NS.NX) GO TO 26 

PRl 

1555 

CC73 


KSU8=0 

PRl 

1556 


C 

CHECK FOR "XTAOLE" 

P.R 1 

1557 

0079 


IF(L0CATE(NSYM,NXTAB,2,5).EQ.C) GO TO 26 

PRl 

1553 


C 

LCCATE ENC OF ARC STRING IN XPRESSION 

PRi 

1559 

CC8C 


OC 23 I=KSTR,50 

PR 1 

156C 

0031 


IF(NRAYm.EO.NRP) GO TO 24 

PRl 

1561 

0:52 

23 

continue 

PRl 

1562 


C 

NC ENC CF ARG. FOUND. 

PRl 

1563 

CC85 


GC TO 29 

PRl 

1564 

0039 

24 

KSTR=I+l 

PRl- 

1565 

CO 8 5 


RETURN 

PRl 

1566 

OCSi 

25 

LP0LL=1 

PR 1 

1567 

0037 


GO TO 7 

PRl 

i5o8 


C 

"XTADLE" FOUND 

PRl 

1569 

CCS9 

26 

KSTR=K5TR-1 

PRl 

157C 

CC69 


RETURN 

PRl 

1571 


C 

ERRORS 

PRl 

1572 


C 


PRl 

1573 


C 

STACK CVEBFLOW 

PRl 

1574 

CC9C 

27 

CALL ERROR 122> 

PRl 

1575 


C 

BAD SUBSCRIPT 

PRl 

1576 

CC9L 

28 

CALL ERROR (23) 

PRi 

1577 


C 

DAD EXPRESSION ARGUMENT STRING 

PR 1 

157S 

0092 

25 

CALL ERROR (24) 

PRl 

1579 
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LPUMP 


DATS = 78178 


1A/AA/A5 


PAGE 0001 


ccci 


function LPUMP (NRAY,KRAY,KLEN tKVAU) 



PRl 

1531 


c 




PRl 

15S2 


c 

<><><><><><><><>0<><><><><>0<><><><>0<><>< 

50000000000PR1 

1563 


c 




PRl 

1584 


c 

PUMP OUT THE INTEGER VALUE (KVAL) FROM ECO ARRAY (NRAY) 

PRl 

1585 


c 

STARTING AT (KRAY) FOR LENGTH (KLEN). UPDATE 

IKRAY) 

AND 

PRi 

1586 


c 

IGNORE BLANKS. LPUMP=0 IF CKAY CONVERSION 



PR i 

15S7 


c 

LPUMP=l IF ERROR IN CCNVcRSICN 


PRi 

1538 


c 




PRl 

1589 


c 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPR 1 

159C 


c 




PRl 

15S1 

OC02 


COMMON /CHARS/ NSP,NLPfNRP,NCaMfN£Q,ND,NPLUSf 

AMINUS 

t L y 

NZERjNXHOPRl 

1592 



1 ) ,NNUM(10 > 



PRl 

1593 

CCC3 


DIMENSION NRAYISO) 



PRl 

1594 


c 




PRl 

1595 

OCOA 


KVAL=0 



PRl 

1596 

CCCS 


LPUMP=0 



PRl 

1597 

CCCii 


DO 3 1=1, KLEN 



PRl 

1598 

0007 


IF(NSAYfKRAY).FQ.NSP) GO TO 3 



PRl 

1599 

ccce 


DC 1 J=l,10 



PRl 

16CC 

CCCR 


IF (NRAYtKRAY).fiC.NNUMCJ) ) GO TO 2 



PR : 

16C: 

0013 

1 

CONTINUE 



pr: 

1602 

oou 


LPUMP=1 



PRl 

1603 

CC12 


RETURN 



PRl 

16C4 

CC13 

2 

KVAL = KVAL>i=10+U-l) 



PRl 

1605 

00 lA 

3 

KRAY=KRAY+1 



PRl 

1606 

cc:5 


RETURN 



PRl 

1607 

C C 16 


END 



PRl 

16C6 



FCRTR'AN IV Ci RELEASE 2.0 


LF>UMP 2 


DATE = 781.76 


iA/44/45 


PAGE 0001 


CC'Jl 


0002 

0003 

CCC4 

CCC5 

0006 

QOOT 

ccce 

0009 

OCiO 

CCll 

0C12 


c 

c 

c 

c 

c 

c 

c 

c 


c 


FL-NCTICN LPLMP2 ( NR AY . KRAY , KV AL ) PRl 

PRl 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRl 

PRl 

PUKP THE (NEXT SUBSCRIPT SURROUNDED NUMBER PROP NRAY PRi 

RETURN VALUE AS LPUMR PRl 

pr: 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRl 

PRl 

CCMMOiJ /CHARS/ NSP , NLP i NRP ,NCQM , NEQt NO , NPLUSiNH t NUS ,NSL iNIER .NXMIBPRI 


ll.NUKdO) . PRl 

DIMENSION NRAY(aO) PRl 

PRl 

LPUMP2=1 PRl 

LRAY=KRAY PRl 

LRAY=LRAY+1 PRl 

IP (LRAY .GT .72 } RETURN PRl 

Ir(NRAYlLRAYl.N5.NRP.AN0.KRAY(LRAYl.NE.NC0M) GOTO I PRl- 

KXRAY=LR AY-KRAY ‘ PRl 

LPU«P2=LPUMP{.NRAY,KRAY.KXRAY,KVAU . PRl 

RETURN PRl 

END PRl 


16CS 

1610 

1611 

1612 

1613 

1614 

1615 

1616 
16 17 
161B 
1615 
162C 
1621 
1622 

1623 

1624 

X62S 

1626 

1627 

1628 

1625 
1630 



- 81 ' 



FCRTfi^^ 1\/ Gl RGLE/iSi: 2.0 


LRFAD 


date = 78178 


lA/AA/45 


PAGE 0001 


CCCl 


0002 

0C3^ 

CCCA 

0005 


CCC6 

0007 


uCC8 

CCC9 

OCiO 

cell 

0012 

0013 

CCIA 

CC15 

0016 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


c 

c 


c 


1 

2 

c 


3 


FLECTION LRcAO (NFIL,LIST) 


PR 1 

16 3! 



PRl 

1632 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPRl 

1633 



PRl 

L63A 

READ A RE’eoRD FROM FILE :^FIU INTO (NRAY ) 


PRl 

1635 

LIST IF (LIST) IS =1. RETLRN CCDES: 


PRl 

1636 

LRcAD=0 NO END CF FILE AND "END" NOT IN COLS 

7-9 

PRl 

1637 

=1 END OF FILE Oft "END" IN COLS 7-9 


PRl 

1638 



PRl 

1639 

<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><> PRl 




PRl 

16A1 

COMMON /NASA/ N SYM( 6 ,300) , KS YM ,NRAY { 30 ) »KRAY 

,NL£N(8},NPT(8) 

,NT(300PR1 

16A2 

.) ■ ■ ■ 


PRi 

16A3 

CCf'KGN /FILES/ NPRINT,NPGM,NOTA,NEQU 


PRl 

16AA 

DIMENSTCN N6ND(3) 


PRl 

16A5 

DATA NcND/lH6,lHN,lHD/ 


PRl 

16 6 6 



PRl 

1667 

SET PCINTER ANC FLAG 


PRl 

1668 

KRAY=7 


PRL 

1669 

LR£AO=C 


PRi 

1650 

READ NEW LINE 


PRl 

1651 

READ (NFIL,1,ENC=3, ERft=3 ) NRAY 


PRi 

1652 

FORMAT (eCAl) 


PRl 

1653 

IF(L IST.EQ. 1) WRITE (NPRINT,2) NRAY 


PRl 

1656 

FORMAT (IX,80A1) 


PRl 

1655 

"END" APPEARS? 


PR 1 

1656 

IF(L0CATE(NRAY,NEN0,7,3) .EQ. 0) LREAD=1 


PRl 

1657 

RETURN 


PRi 

1658 

LR£AD=1 


PRi 

1659 

RETURN 


PRi 

166C 

END 


PRl 

1661 


■ 82 - 



FDRTRiN IV G1 RELEASS 2.0 


LSHUV 


CATE = 78178 


14/44745 


PAGE OCOl 


0001 


CCC2 

CCu3 

0C04 

CCC5 

CCl.6 

OGCV 

cccs 

ccc? 


OOiO 

pen 

CC12 


0013 

CCIA 

CC15 

0016 

0017 

CClfi 
CC15 
OCiO 
CC21 
CC22 
0023 
0C24 
CC25 
CC2£ 
0027 
002 3 
Gcig 
00 30 
0031 
CC32 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

1 


2 


3 

4 


C 

C 

c 


FUNCTION LSHUV { KPOP , I POs’.KPCPC) pi^i 

PRl 

Ooooooooooooooooooooooooooooooooop:^! 

PRi 

PUSH THE CURRENT SVySOL LCCATCRS INTO (KSHUV) IF THE PR^ 

NEXT SYMBOL IS LOCATED IN THE LIST (KLOCA) OF SYMBOLS PRi 

PRi 

<><>0<><><><><><><><><><><><><><><><><><><>^><><;>^><.>^;> OOOOOPRl 
lyPLICIT P£AL*S (A-H,0-Z) 

COPHCN /PUSHED/ KLOCA ( 3 ,4 CO 1 ,KLOC , NPUS H ( 1000 1 , KPUS H, KPHAX 
CaPHON /ShGR/ NSHUVl 500) , KSHUV iKST,KEND ,XSUB fKEXTO jKPAR 
CCPMCN /CNTRL/ LIST,Ll,CTROL(i60),KXN{3J 
COPMCN /FILES/ NPR] NT,NPGP,MCTA,nEQU 


' PRi 
PR 1 
PRI 
PRI 
PRI 
PRi 


LSFUV^I 

KyGO = Pf;DlKFCPTlOOO ) 

GO IF ECLATIC.N FCUKD 

IF(LOCTST< IPOS,KPOP).EO.O) GO TO 1 

NC ECUiTiCN FOR THIS VALUE 

USE LAST VALUE 

]P0S^K«O0 

KPCPO=KP0P 

OC TO 2 ' 

PUSH sue SCRIPT. POINTERS, VALUE AND SOM-PROCUCT'S FLAG 
LShUV=0 

IF(KSfUV+5.GT.S00 ) CD TO 5 

KSriUV=KSHLV+l 

NSHUV(KSHLV)=KSLB 

KSFUV=KSPUV+i 

NSHUV (KSHUV )=KST 

KSHUV=KSHLV+1 

NSHUV (KShUV)=KEND 

KSHUV=KSHUV+l 

N'SHUV(KSHLV) =KPCPC 

KShUV=i<SHUV+l 

NSFUV IKShUV )=KPAR 

KFC?C=KPCP 

IF(HST.cC.O) RETURN 
KS5=KSHtV-4 

URITE (NPRINT,3) I POS , KPOP ,KSUS , Kc XTD, L SHUV, KSHUV 
IFIKS5.GT.0) WRITE (NPRINT,4) (NSHUV (I), I=KS5, KSHUV) 
FORMAT ( 7H LSHLV ,619) 

FORMAT (1X,5I1C) 

RETURN 

ERRORS 

SHUV STACK OVERFLOW 


PRi 

PRI 

PR I 

PRI 

PRI 

PRi 

PRi 

PRi 

PRI 

PRl 

PRI 

PRl 

PRl 

PR 1 

PPl 

PRl 

PRi 

PRl 

PRl 

PRl 

PRl 

pRl 

PRl 

PR 1 

PRi 

pai 

PRl 

PRl 

PRl 

PRl 

PRl 

PRi 

PRl 

PRl 


1662 

1663 

1664 

1665 

1666 
1667 
1666 
1669 
16 70 
1671 
16 72 

1673 

1674 

1675 

1676 

1677 
16 78 
1679 
160() 
1681 
1622 
1633 
1684 
1635 
1686 
1637 
i6ae 
1639 
169C 

1691 

1692 

1693 

1694 

1695 

1696 

1697 
1693 
1699 
17CC 

1731 

1732 
17C2 

1704 

1705 
17C6 

1707 

1708 

1709 
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CATS = 78178 


14/4V45 


PAGE 0002 


0033 

0C34 


5 CALL ERROR {33) 
END 


PRl 1710 
PRl 1711 
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FC,5r?i^ JV G1 

'bco: 


CCC2 
0003 • . 
OC 0^ 


OOC-5 

0CC6 

CCC7 

ccca 

OCC', 
CC' ■» 
CCl '■ 
0012 
0C13 
CCl'i 
0015 
C01& 
GCi7 
CC IS 
0019 
CC20 
CC21 


RELEASE 2.0 


C 

C 

C 

C 

C 

C 

c 


3 

4 


date = 7S178 l,iV/V-/45 

FUNCi^iOiV NUMB (NRAY ,KST>KEN, I ) ’ '■ 

ooixxxxxxxxxxx;^.^^.^ O <><><><><XXX><><x><xxx> 
DECODE MJKSER AT KST IN NRAY UN ITL KEN ' 

<><XXXXXX>0<XXXXXX><><X><><><XXXX><><XX><><> 

common -/FILES/ NPRiNT,NPGf',NCTA,NEOU 

C'JiMMD-N /CNTRL/ LI ST, L 1 1 C TROL ( 160) , KXM{3 > ;! 

'^Sf’''^>-^'f^2’'P»NC0M,'f£0,N0,NPLUS,NMINUS,NSL,NZER,NXM 

DIMENSION NRAY(8C> 

NUMB=-3 

DC A I=KST,KEN 
IFd.LE. 72) GO TC 2 
READ {NPGMjlJ NRAY 
FORMAT (80AU 
I =7 

KEN=KEN-72-t-6 

CCNTINUe 

I F(NP AY ( I ) . FQ.NSP ) GO TC A 
DO 3 J=l,10 

IF (NRAY( I ) .NS.NUMSIJ >1 GO TO 3 

NLM8 = NU.''e'!‘lO + J-l 

CCNTI NL'E 

CONTINUE 

RETURN 

END 




. "^Rl 

1712 

' . PR 1 

1712 


171A 


1715 


1716 

r R 1 

1717 

<>’-;ri 

1718 

. * • 
,* »■ % \ 

171? 

•' PRl 

172C 

J.^Rl 

1721 

( 3PRi 

1722 

PRX 

1722 

PRl 

1724 

PRl 

1725 

PRl 

1726 

PRl 

1727 

PRl 

1728 

PRl 

172? 

PRl 

173C 

PRl 

1731 

PRl 

1732 

PRl 

1723 

PRl 

1734 

PRl 

1735 

PRl 

1736 

PRl 

1737 

PRl 

17 3S 

PRl 

173? 

PRl 

174C 


PAGE 0001 
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C 

PROSE 2' ' 

PR2 

CCOl 


c 

<<<<<<<<<<< <«<<<<<««««««««>>»»>»»»»»»»»»»»»p:^2 

0002 


c 


PR2 

0003 


c 

PRCSe PART TWO 

PR2 

CCC4 


c 


PR 2 

CC05 


c 

COPYRIGHT JOHN E. O'REILLY, JR. 1977 

PR2 

0006 


c 

10/17/77 

PP.2 

00 07 


c 

CCPYRIGKT JCHN E. C’REILLY, JR. 1978 

PR2 

CCC8 


c 

04/05/78 

PR2 

0009 


c 


PR2 

0010 


c 

<< <<<<<<<<<< <<<<<< <<<<<<<<<<< <<««»» »>»>0 >>>>>>>>>> >>»»»»PR2 

CCll 


c 


PR2 

0012 

OOol 


IMPLICIT REAL^oU-HfO-Z) 

PR2 

0013 


c 

STANDARD SCEPTRE COMMON AREAS 

PR2 

0014 

CCC2 


CCMMON /CNTRLS/ TIM£,S1 (21 ) jXNCPRC, S 2, XNHtAO, S 3 ( 6 ) , XRUNNO, S4( 11) ,X=>R2 

0C15 



IWRTSM , S5( 46) ,XNEY,XN5 7DT,XNBC ,XNBR ,XNBL ,X NLC, XNLR, XNLL, XNJX, XNJ8 DTPR2 

0016 



2,XUJ0,XNJ9, S6(2 ), XNM33, XNM36 , XNM66 , XNDP , XNDPD , XNELTS 

,S7(50) PR2 

0017 

CCC3 


CCht'Ch /CURNTS/ C(l) 

PR 2 

0018 

CCLA . 


COMMON /DERIVS/ D(l) 

PR2 

0019 

CCC5 


COMMON /0PNAM3/ CPS(l) 

PR2 

002 0 

C3C6 


CCMMCN /ELNAMS/ LELS (6,300) 

PR2 

CC21 

CCC7 


COMMON /HDCRDS/ HfcAD(9,ll) 

PR2 

0022 

0C03 ' 


COMMON /NAMES/ XN(1) 

PR2 

0023 



COMMON /CUTPTS/ 0UTR(5,1) 

PR2 

00 24 

CCIC 


COMMON /CtTPUr/ NCUTP(6,50) 

PR2 

0025 

ooi: 


COMMON /PRMTRS/ P(l) 

PR2 

0026 

CC.2 


COMMON /SAVICT/ ClCT(l) 

PR2 

0027 

CC13 


CCMMCN /SRCCRV/ S (1 ) 

PR2 

0C2S 

CCIA 


COMMON /TABLES/ T(l) 

PR 2 

0029 

0015 


COMMON /TAPES/ NPRT ,NOTA ,MPG M 

' PR2 

0030 

C C 10 


COMMON /VCLTGS/ V(l) 

PR2 

00 31 

CC.7 


COMMON /XMUTL/ XM(i) 

PR2 

0032 


c 

PR0SE2 common area 

PR 2 

0053 

0C13 


COMMON /PROSE/ A(50,50),B(50,50),L3R,L3S,L4P,L4S,NSV 

,NEJ,NEL,NIND(PR2 

•JO 3 4 



13C0) ,IPERF(50) 

PR2 

C035 


c 

' 

PR2 

0035 

0C17 


DIMENSION CTR0LS(16C), 5R(5C), EK50), XSR(50)t XEK50) PR2 

00 37 

CC2: 


DIMENSICN NLIN(80), NSRCDR(6), NTABLE(6), N0UTPT(6), 

NUMS(IO) PR2 

CC38 

CC 2 1 


OIMENSICN NLT(6.), NLT2 (6) 

PR2 

CC39 

0C22 


DATA NSRCOR/ IHS, IHR, IHC , IHD , IHR ,1HV/ 

PR2 

0040 

C323 


DATA NTABLE/lHTflHA.lHB, U-L, IHE, IHS/ 

PR2 

CC41 

CC24 


DATA i\riTPT/iHr,iHU,lHT,lHP,lhT,lHS/ 

PR2 

0042 

0025 


DATA NCOM/lH,/,iNBL/lH / 

PR2 

0043 

C 32 


DATA NUM.S/lHO, 1 H 1 , 1H2 , IH 3, 1H4, 1H5, 1H6, 1H7, 1H3 , 1H9/ 

PR2 

0044 

CC27 


DATA M/lhI/,,NV/lHV/ 

PR2 

CC45 

0C28 


EQUIVALENCE ( TI ME ,CTROLS( 1) ) 

PR2 

0046 


c 


PR2 

0047 


c 

BEGIN MAINLINE CODE 

PR2 

CC48 
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FORTRAN' IV GI RELEASE 2.0 


MAIN 


DATE = 7BL7B . 




PAGE 0002 


CC29 

CC30 

CC31 

CC32 

C033 

CC34 

0033 

'CC36 

CC37 

0033 

CC39 

0040 

CC41 
004 3 
0043 

CC44 • 

0045 
CC46 
C C 4 7 
0043 

CC49 
0050 
0C5i 
CC52 
CC 5 3 
CC5<* 

'CC55 
0056 
CC57 
C'0 5S 


0050 

0060 

0061 

0C62 
CC£3 
0064 
C C 65 


C PR2 

•C SET I/C CEVICE CODES PR 2 

NPRT=£ PR2 

Nt)TA=l PR2 

NPG«=2 PR2 

CALL SETLP2 ( N4C , LELSf NSRC, NTA0) PR2 

C PR2 

IF(XWRTSM.EQ.1 .000 ) WRITE (NPRT.l) CTROLS PR2 

1 FORMAT (21H <><><>C0NTR0LS<>OO,/,32(lX. iP5D14.7,/n PR2 

2 FORMAT :i?5Di4.7) . PR2 

C READ FFACiNG CARDS PP2 

NH£AD=XNHtAC PR2 

IF ('MHEAD. AE, 01 READ IN0TA,4) HEAD PR2 

IFUWRTSH.eO.l.CDO.AND.NHEAD.NE.OI WRITE tNPRTi3) ( (HEAD! I , J ) , 1=1, PRl 
19>,J=1,NH£AC) PR2 

3 FORMAT I26H OOOHEAOING C ARCSOOO , / , 11 (IX ,9 A8, / ) » PR2 

4 FORMAT (9A8) PR2 

C READ ELEMENT VALUES PR2 

NcL=XNELTS . PR2 

RLAO INOTA,2) I XfU U , I =1 , ME L ) . PR2 

IFIXWRTSN.EO.I .COO) WRITE (NPRT.S) ( XN( I) ,I =1 rNEU PR2 

5 FORMAT (27H OOOELEMENT V ALU ES <><><> , / , 60 t IX, IP 5D 14 . 7, /)> PR2 

C READ CLEMENTS NAMES , PR2 

IFINAC.eo. 160) READ INDTA,?) ( ( LEL S ( I , J ) , I =1 , 6 ) , J=1 , NE L ) PR2 

IFfXWRTSM.EC.l.OCU ! WRITE (NPRT,6) ( (LELSII ,J I = i,'6) ,J = 1,NE-L) PR2 

6 format (26H OOOcLcMEKT NAMESOOO , /, 30 ( iX . 72 A1 , / ) ) PR2 

7 FORMAT (72A1) PR2 

C READ ELEMENTS VOLTAGES AND CURRENTS {INITIAL CONDITIONS) PR2 

READ (NDTA,2) « V ( I) , I =1 , ImEL ) PR2 

IFIXW'ITSM.L-Q.l. COO) WRITE (NPRT.B) t V ( II , 1 = 1 , NEL ) PR2 

8 FORMAT (25H 'OOOi/C VDL TA GES <> OO , / , 60 UX, I P 5D14.'7 , / )) PR2 

READ iNCTA,2f) ( C ( I ) , 1= 1 ,N EL I ' PR2 

IF tXWRTSM.cC.l.CDO) WRITE. INPRT, 91 I C t II , 1 = 1 , N EL I PR2 

9 FORMAT (25H OOOI/C CURREN TSOO O , / , 60 { IX , 1 P 5014. 7 ,/ ) ) PR2 

C READ f'UTUAL VALUES PR2 

•* NHUT = XKM33+XNM36+XNM66 • PR2 

IFINMUT.NG.C) READ (N0TA,2) ( XM ( I ) , I =1 , NMUT ) PR2 

IFtXWRTSy.EC.l .OOO.ANO.NMUT.NE.OI WRITS INPRT,10) ( XHt I I ,1 =1 ,NMUT) PR2 
1C format (20H <><><>MUTUALS<><>0,/,20 (IX,1P5C14.7,/) I PR2 

C • READ SOLRCS DERIVATIVES PR2 

IFINSRC.NS.CI READ INDTA,2) f S 1 1 ) , 1= 1, NSRCI PR2 

IF(XWRtS.v.£C.l .OC.3.ANO.NSRC.NE.O > WRITE (NPRTrll) ( S ! ! ), 1= l.NSRC I PR2 

11 FORMAT (25H OOOSRC. DERI VS. <><><>,/ ,10 (IX , 1 P5 E14 .7 ,/) i PR2 

C READ STATE VARIA8LS SAVES I /C ■ S PR2 

NIC=XMC'i'XNLL+Xf':£C+XNBL PR2 

IFCNIC.Nc.Ol READ lN0TA,2) ( C I CT (1 ) , I =1 , N IC ) PR2 

IF(XWRTSM.EQ.l.CDO.A.NO.NIC.NE.O) WRITE (NPRT,L2) (CICT ( I ) , 1=1 ,NIC J PR2 

12 FO.RKAT (22H OOOSTATE I/COOO, /, 10 ( IX , IP 50 14 .7 , / ) I PR2 


0049 

CC5C 

0051 

0052 
CC53 
CC54 

0055 

0056 
GC57 
0058 
0J59 
0C6C 
0C61 
0062 
0C63 

0064 

0065 
006 6 
CC67 
006S 
U'369 
CC7C 

0071 

0072 
CC73 
CC74 
0075 
0976 

0077 

0078 

0079 
CC3C 
QC81 
0082 
CC32 
CC84 
0085 
0C86 
0067 

0038 

0039 
CC9C 
OC'91 
0092 
0C93 
GC94 
0095 

■0096. 


- 87 - 


ORIGINAL PAGE I& 
OP POOR QUAury 



FO'^TRA.N IV G1 RELEASE 2.0' 


^^AIN 


DATE =< 70178 


1A/A4/A5 


PAGE 0003 



C • 

COMPUTE NUMBER OF STATE VARIABLES 

PR2 

0097 

CC66 


NSV=X'JEC•^X^LL 

PR2 

• CC98 


C 

READ DEFIAEO PARAMETER DERIVATIVES 

PR2 

0099 

0067 


NCPO=XNOPD 

PR2 

0100 

0063 


IFINDPD.NE.C) READ (N0TA,2) ( D (I -fNSV- 1 ) , I = 1 ,NDPO ) 

PR2 

OlOl 

CC65 


IFUKRTSM.SC.l .OCJ.AND.NDPD.NE.O ) WRITE (NPRT,13) ( OII-fNSV-l ) , 1=1 

,PR2 

0102 



INOPD) 

PR2 

0103 

0C7C 

13 

FORMAT (31H OOODEF .• PARM. 06R1 VS. <><><>,/, 10 t IX t1P5D14. 7, /) ) 

PR2 

0104 


C 

READ STATE VARIABLE I/C'S 

PR 2 

C1C5 

CC7I 


nll=x.':ll 

PR2 

01C6 

0072 


Ir(NLL.Nc.O) READ (NDTA,2) ( C I C T (I > , I =1 , NLL ) 

PR2 

0107 

CC73 


IF(XWRTSM.5Q.1.00w.AND.NLL.NE.O) WRITE (NPRT,14) ( C ICT( I ) , I =1 ,NLLJ PRa 

0103 

CC74 

14 , 

format U9H OOOLL 1/ COOO t / » 10 { IX t IP 5 D14 .7, / »} . 

PR2 

0109 

0075 

. 

NI=NLL-1 

PR 2 

0110 

0076 


NBC=XNPC 

PR2 

0111 

CC77 


IF(NQC.AE.O) READ (N0TA,2) ( C IC7 { UN 1 ) , 1= 1 ,NBC > 

OR2 

0112 

OC78 


IFtXWRTSM.EC.l.CDO.AND.NBC.NE.O) WRITE (NPRT,15) ( C ICT ( UNI ) , 1= 1, 

NPR2 

0113 



ISC) 

PR2 

0114 

c;7s. 

15 

FORMAT <19H OOOBC 1/ COOO , /, 1 C t IX, IP 5D 14 . 7 , / ) ) 

PR2 

0115 


C 

READ DEFINED PARAMETERS 

PR2 

0116 

0000 


NDP=XNDP 

PR2 

0117 

0C8. 


IFI.VDP.NE.O) R = A0 (NDTA,2) ( P ( I) , I = 1 , NDP ) 

PR2 

0118 

C C 2 • 


IF (XWRTSM. EC.: .OC'J . A.NO.NDP.NE.O ) WRITE (NPRT,16) (P ( I ), I=1,NDP ) 

PR2 

0119 

0005 

16 

format I24H OOODEF. PAR MS. <><> O ,/ ,20 ( IX , 1 P5 D1 4 .7 , / ) ) 

PR2 

0120 

C 06-V 


Ir(NDP.NS.O.AND.riAC.5Q.16C) READ (NDTA,17) (0 PS { I ) , I =1 , NDP) 

PR2 

0121 

CC55 


IFIXWRTSM.EC.l .OC d.AND.NOP.NE.O.AND.NAC.£Q.160) WRITE (NPRT,18) (DPR2 

0122 



IPS(I) ,1=1, NDP) 

PR2 

0123 

0056 

17 

FORMAT 112A6) 

PR2 

0124 

COST 

13 

FORMAT (29H OOODEF. PARM. N AM ESOOO, / , 9 ( 1 X , 12A 6 , / ) ) 

PR2 

0125 


C 

READ TABLES 

PR2 

0126 

003 7 


IFINTAB.ME.C) READ (N0TA,2) ( T 1 1 ) , I =l , NTA8 ) 

PR2 

0127 

GCE9 


IFIXW.lTSM.cC.l.CDO. AfsjD.NTAB.NE.O) WRITE (NPRT,19) (T { 1 ) , 1= 1, NTAB ) 

PR2 

0128 

CC90 

19 

FORMAT ( 19H 00<>TABLES<>00,/,50(1X,1P5D14.7,/) ) 

PR2 

0129 

• 

C 

READ OUTPUT RECUESTS 

PR2 

G13C 

ccs: 


NUPRC^XNCFRC 

PR2 

0131 

0C-:»2 


READ |NDTA,17) t ( OUTR ( I , J ) , I =1 , 5 ) , J = 1 , NC PRC) 

PR2 

0132 

OC73 


IF(XWRTSM.EO.I.ODO) WRITE (NPRT,20) ( { 0 UTR < I , J ) , 1=1 , 5 ) , J=1 , NOPRQ) 

PR2 

0133 

CCS4 

2C 

format (28H OOOOUTPUT RE CUESTS <><><>,/, 10 ( IX, 12 A6, / )) 

PR2 

0134 

0095 


N0PRQ=N0PKQ+i 

PR2 

0135 

CC95 


READ (NCTA,7) ( (NOUTP ( I , J ) , I = l , 6 ) , J=1 ,NOPRO) 

PR2 

0136 

CC97 . 


IF(XWRTSM.SC.l.CnO) WRITE (NPRT,21) ( 1 NOUTP ( I , J ) , 1= 1, 6 ) , J= 1, NOPRO ) PR 2 

0137 

CC93 

21 

FORMAT (2CH OOOOUTPUTSOOO ,/ ,1 0 ( 1 X ,72 A1 ,/ ) ) 

PR2 

0138 


C 

ZERO OUT DERIVATIVES 

PR2 

0139 

CCS9 


OC 22 1=1, NSV 

PR2 

014C 

C’lCC 

■22 

Dn)=v. CDC 

PR2 

0141 


C 

CLEAR INDEPENDENT VARIABLE TABLE' 

PR2 

0142 

0101 


DC 2 3 1=1, NSL 

PR2 

0143 

CIC2 

23 

NINO(I)=0 

PR2 

0144 
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C1C3 

0104 

0105 
C1C6 
C1C7 
0103 
C1C9 

one 

on: 

0112 

0113 

0114 


FORM INCEPENEDEWT VARIABLE TABLE 
Ni:J = U 

00 26 I=2,NOPRO 
00 25 J=i,NEL 
OC 24 K=1 ,6 

IF(NOLTP(K,n;A£. LELS£K,JJ) GC to 25 
continue 
MN 0(J) = 1 
■ NEJ=N>;J+i 
. XN(J)=C.COO 
GO TO 26 
CCNTI NUE 
CONTINUE 

COMPUTE OFFSETS TO STATE VARIABLE NAMES 
• L4S=XNE7CT+1 
L4P = >Ij<iRC+L4S-l 

L JS=X JETOT + XNBC + XNflR+XNSL + XNLC+XNLR+l 
L3P=XML + L3S-1 
CLEAR STATE VARIABLES 
IF(L3S.GT.L3PJ go TO 28 
CO 27 •I=L3S,L3R 

c<n=c.o 

IFIL4S.0T.L4P) GO TO 30 
CC 29 I=L4S,L4P 
V(n=J.OC^ 

PcRFOiiy EIGAN ANALYSIS 
CALL EIGAN {SR, El) 

PRINT FEADER 
WRITE {NPKT,3D 

FORMAT { 27HiNASA”LEWI S RESEARCH CENTER, /,27H 

1 ,//,27H prose program ,/,27H ' 

I ^ 1-'/®^ UNIVERSITY OF SOUTH FLORIDA ,/,28H ”IZZ_ 

<<HEAan,0).I=1.9,,0 = l.NHEAD) 

WRITE :nPPT,33) 

t/,21H state variable names,/, 21H 

WRITE 13 STATE VARIABLE NAMSS ' 

I F (L3S.GT.L3P) CO TO 37 

00 35 I=L3S,L3P 

00 34 il=l,5 

NLT2 ( I r+1 )=LELS (11,1) 

NLT2{1)=NI 

CALL BACK fNUT2,NUTl 
WRITE (NPRT,35 ) NUT 
format (1X,6A1) 

WRITE C4 STATE VARIABLE NAMES 


PAGE 0004 


PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
-PR2 
PR2 ' 
-PR2 ' 
PR2 
PR 2 ' 
PR2 < 
PR2 I 
PR2 I 
PR2 ( 
PR2 ( 
PR2 ( 
PR2 { 
PR2 C 
PR2 < 
PR2 C 
PR2 C 
PR2 C 
PR2 C 
PR2 0 


2 0145 
2 0146 
2 0147 
2 0148 
2 C149 
2. 015C 
2 0151 
2 0152 
2 0153 
2 0154 
2 0135 
2 0156 
2 0157 
! 0158 
2 0159 
2 016C 
: 0161 
’ 0162 
! 0163 
^ 0164 
01 o5 
: 0166 
01o7 
0168 
0169 
017C 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 
018C 
0101 
0132 

0183 

0184 
0135 

01d6 

0187 

0186 

0189 

0190 

0191 

0192 
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ClAC 

3.7 

I F IL4S.GT.L4P) GO TO 41 


PR2 

0193 

OiAl 


00 39 I=L4S,L4P 


PR2 

0194 

0142 


CO 38 11=1,5 


PR2 

0195 

:iA3 

38 

NUT2 m + 1 ) = LSLS ( 1 1, I) 


PR2 

0196 

C1A4 


NLT2( 1)=NV 


PR2 

0197 

01A5 


CALL BACK INUT2,.NUT) 


PR2 

0198 

C1A6 

39 

, WRITE (APRT,40 ) NUT 


PR2 

0199 

C1A7 

40 

FORMAT (iX,6All 


PR 2 

02CC 


C 

WRITE 6 MATRIX HEADER 


PR2 

0201 

0148 

41 

WRITE (NPRT,421 


PR2 

0202 

OlA^ 

42 

FCRMAT |/,9H r ,/,9K MATRIX A,/,9H 

j 

PR2 

02 C 3 


C 

WRITE A MAT.RIX 


PR2 

0204 

0150 


DO 43 1=1, NSV 


PR2 

02 J5 

0151 

43 

WRITE (AP'RT,44) ( A { I, J ) , J = 1 , NSV ) ' 


PR2 

02C6 

C152 

44 

FORMAT (1X,1P7014.5) 


PR2 

Q2C7 


C 

WRIT- INOEPENDENT SOURCE HEADER 


PR2 

0208 

0153 


IF (NEJ .EC.O > on TO 50 


PR2 

0239 

015^ 


WRITE (KPRT,45) 


PR2 

02 IC 

0155 

45 

FORMAT I/,2CH ,/,20H INDEPENDENT SOURCES ,/ ,20H 

PR2 

0211 



1 ) 


PR2 

0212 

0156 


DC 47 1=1, NcL 


PR2 

0213 

0157 


IFlNlKOm.NE.l) GO TO 47 


PR2 

0214 

OiST. 


CALL BACK (LELS ( 1, I ) ,NUT) 


PR2 

0215 

0155 


■WRITE (NPRT,46) NUT 


PR2 

0216 

CH5 

46 

FORMAT tlX,6Al) 


PR2 

0217 

0161 

47 

CONTINUE 


PR2 

0218 


C 

WRITE R MATRIX HEADER 


PR2 

v3219 

C162 


WRITE tNPRT,48l 


PRZ 

022C 

0163 

48 

FO.RHAT i/,9H -,/,9H MATRIX B,/,9H 

1 

PR2 

0221 


C 

WRITE 8 MATRIX 


PR2 

0222 

0164 


OC 49 I=i ,NSV 


PR 2 

0223 

0U5 

49 

WRITE tN»RT,44) ( B ( I., J ) , J=1 , NS J ) 


PR2 

0224 


C 

WRITE EIGENVALUES 


PR2 

0225 

0166 

50 

WRIT? INPRT,51) 


PR2 

C226 

0167 

S i. 

FORMAT (/,24M ,/,24H 

EIGENVALUES OF MATRIX 

PR2 

0227 



1A,/,24H ,/,14X,9HREAL 

FART,10X,10HIMAG. PARPR2 

0228 



2T,/,14X,9H ,10X,10H >. 


PR2 

02 2 9 


c 

PERM LARGEST EIGENVALUE AND PRINT THEM 


PR2 

023C 

0163 


XL=O.COC 


PR2 

0231 

OloO 


00 53 I=i,NSV 


PR2 

02 32 

0170 


write :nprt,52I sRm,5im 


PR2 

0233 

0171 

52 

FORMAT 1 i2X,lPEl4.7,5X,lPcI4.7) 


PR2 

0234 

0172 


XV=OS.JRT( ERI n #ER{ n+El (I )X<£H I ) ) 


PR2 

0235 

Cl 73 


IF(0A3StXL).GT.CABS(XVl) CO TO 53 


PR2 

02 36 

0174 


XL=XV' 


PR2 

0237 

0175 


IL = I 


PR2 

0238 

0176 

53 

CCNT LIUS 


PR2 

0239 


C 

DC SENSITIVITY ANALYSIS 


PR2 

024C 
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0177 


WRITE (NPRT,54) IL,XL 



PR2 

0241 

0178 

•54 

FORMAT (/,21H ,/,21H SENSITIVITY A\!AL YS I S, /, 21PR 2 

0242 



iH ,/,21H HAXIHUH EIGENVALUE 

K, 13,1311 

MAGNITUDEPR2 

0243 



2 = ,1PE14.7,//,21H ELEMENT SENS! TI VI TY, / , 21H — 




0244 



3) 



PR2 

0245 

017<5 


EGL=O.CDO 



PR2 

0246 

0180 


IX=0 



PR2 

0247 

0181 


CO 61 I=1,NEL 



PR2 

024 8 


C 

IGNORE IF INOEP. SOURCE 



PR2 

0249. 

0182 


IFiNINOd ) .EO. 1) GO TO 61 



PR2 

025C 


C 

IGNORE IF- ZERO VALUED 



PR2 

0251 

Cl 63 


IF(XN(I) .E0.0.9CO) GO TO 61 



PR2 

0252 


C 

SAVE OLD VALUE 



PR2 

0253 

0154 


SV=XN( n 



PR2 

0259 


C 

. BUMP VALUE BY 10 PER CENT 



PR2 

0255 

0163 


XMI)=XMn*1.0lCO 



PR2 

0256 

0156 


CALL’EIGAN {X£R,XEI) 



PR2 

0257 


c 

CCMPUTC SENSITIVITIES 



PR2 

0258 

0167 


EGU=nSCRT IX£R( TDVXER ( IL ) +X E I ( IL ) ■S'XE 1 1 ID )-XL 



PR2 

C2 5 9 

0U8 


•EG=cGL=!‘iCC. ODC/XL 



PR2 

026C 

0189 


IF(EG.rT£.“1.030C.ANO.£3.LE.-0. 9700> GO TO 56 



PR2 

0261 

CISC 


IF{EG.Ge.0.97D-'>.AND.£G.L£.1.03C0) GO TO 58 



PR2 

0262 

cis: 


CALL DACK < Ld LS (1 , I ) f NUT ) 



PR2 

0263 

0192 


WRITE (NP.RT,55) NLT.EG 



PR2 

0269 

0193 

55 

FORMAT <2X,6A1, 5X,F6.3) 



PR2 

0265 

0194 


IFIDAESieO.LS.CABSIcGL) ) GC TC 60 



PR2 

0266 

0195 


EGL=EG 



PR 2 

0267 

0196 


NL=I 



PR2 

0268 

0197 


GC TO 60 



PR2 

0269 


C 

GOT SENS APPROX. -1 



PR2 

027C 

0198 

56 

CALL BACK ( LELS( 1,1) .NUT) 



PR2 

0271 

0199 


WRITE (NPRT,57 » NUT, EG 



PR2 

0272 

G2C0 

57 

FORMAT (2X.6A1 ,5X,F6.3,5X,46H4** ELEMENT MAY BE 

REMOVED 

BY 

SETTINGPR2 

02 7 3 



1 TO ZERO. ) 



PR2 

0274 

0201 


NIND( I)=-2 



PR2 

0275 

C2C2 


IX=i 



PR2 

0276 

0203 


GO TO. 6C 



PR2 

0277 


C 

GOT SENS. APPROX 1 



PR2 

02 78 

G2C4 

58 

CALL BACK (LtLS(l,I),NUT) 



PR2 

0279 

0205 


WRITE (NPRT.59) NUT, EG 



PR2 

023C 

0206 

59 

FORMAT { 2X,641,£X,F6. 3,5X,5CH4K« ELEMENT HAY 8E 

REMOVED 

8Y 

SSTTINGPR2 

0281 



1 TO infinity.) 



PR2 

0232 

02C7 


NIN0m=2 



PR2 

0283 

0203 


I X=1 



PR2 

0284 


C 

RSSTORP VALUE 



PR2 

0285 

C2C9 

60 

XN(I)=SV 



PR2 

0286 

0210 

61 

CONTINUE 



PR2 

0287 

0211 


IFClX.EQ.l) STOP 



PR2 

0288 
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62 


CALL BACK (LLLSa.NL) 

WRITE (NPRTf62) NUT 
FORMAT (/j35H NO SENSITIVITIES 
IBs REKOVEC (I.S. SET TO ZERC 
29HSET 
STOP 
END 


rNUT) 


TO INFINITY IF SENSITIVITY > ZERO.')) 


PR2 

PR2 

WITHIN 3S OF 1 .,/,1Xt6A1,40H HAY PR2 
IF, 19H SENSITIVITY < ZERO , /, 29X, 3PR 2 

PR2 
PR2 


02S9 

0290 

0291 

0292 

0293 

0294 
C295 
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ccc: 


SU8R0UTINS ATEIG { A, RR f ft I » I AN A, lA ) 

PR2 

02 96 


C 


PR2 

0297 


c 

<<<<<<<<<<<<<<< <<<«««««« <<«<»1'>»>»»»»»»»»»»»»PR2 

0298 


c 


PR2 

02S9 


c 

COMPUTE eigenvalues CF UPPER HESSENBERG KATRI): A 

PR2 

03CC 


c 

M=OftO£R OF f/ATKIX A 

PR2 

0301 


c 

Rft=R£AL eigenvalue RESULTS ARRAY 

PR2 

C3C2 


c 

RI =IPAGINARV EIGENVALUE ARRAY 

PR2 

03C3 


c 

IANA=COMPUTE PETHID ARRAYIC'OR 2=DIRECT) 

PR2 

OJDA 


c 

IA=FIRST DIMENSION OF A MATRIX 

PR2 

03.35 


c 


PR2 

03Ce 


c 

IBM SSP ROUTINE 

PR2 

0307 


c 

■* 

PR2 

0308 


c 

«««<«««««««««««<<««»>> >>>>»»>>>>>>»> >»>»»»P'^Z 

C3C9 


c 

. 

PR2 

031C 

0C02 


IMPLICIT REAL*e( A-H,0-2> 

PR2 

0311 

GCC5' 


OIMENSICN AU). RRlDt Rltl)* PRft(2)t PRH2). lANAll) 

PR2 

0312 

ccc.« 


INTEGSR’i'A P.PlfC 

PR2 

0313 

. 

c 


PR2 

031A 

0C05 


E7=l.G‘=-8 

PR2 

0315 

CCC6 


E6=1.0E-6 

PR2 

C316 

CCC7 ■ 


S 10=l.nE-10 

PR2 

0317 

0C08' 


DELTA=0,5 • 

PR2 

0318 

CCC? 


MAX IT=^0 

PR2 

0319 


c 


PR2 

032C 


c 

ISITIALI2ATICN 

PR2 

0321 


c 


PR2 

0322 

CCIG 


N=K 

PR2 

0323 

CCil 

1 

NI=N-1 

PR2 

032A 

00 1: 


IN =N 1^1 A 

PR2 

0325 

C C i 3 


NN=IN+N 

PR2 

0326 

CC 1 A 


IF(Ni) 2,69,2 

PR2 

0327 

0015 

2 

NP=N + 1 

PR2 

0328 


c 


PR2 

0329 


c 

ITERATION COUNTER 

PR2 

033C 


c 


PR2 

0331 

0016 


IT=0 

PR2 

0’3 3 2 


c 


PR2 

0333 


c 

RCC'^S CF THE 2ND ORDER MAIN SUBMATRIX AT THE PREVIOUS 

PR2 

0334 


c 

ITERATION 

PR2 

0335 


c 


PR2 

•0336 

c:i7 


DC 3 1=1,2 

PR2 

0337 

CCi? 


PRRn)=0.0 

PR2 

0338 

0.210 

3 

PR II I_) = C,C 

PR2 

033'9 


c 


PR2 

0340 


c 

LAST TWO SUPCIAGONAL ELEMENTS AT THE PREVIOUS ITERATION 

PR2 

0341 


0 

, 

PR2 

0342 

CC20 


PAN=0.0 

PR2 

0343 
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00 Zl 


PAN1=C.0 



PR 2 

0344 


. C 




PR2 

0345 


C 

ORIGIN SHIFT 



PR2 

0346 


C 




PR2 

0347 

0022 


R=C.O 



PR2 

0348 

0023 


S=0.0 



PR2 

0349 


C 




PR2 

035C 


C 

ROOTS CF THE LOWER MAIN 2 6Y 

2 SU8MATRIX 

PR 2 

03 51 


C 




PR2 

0352 

COZA 


N2=N1-1 



PR2 

0353 

CC25 


IM = IN-IA 



PR2 

C354 

0026 


NN1=II-I1+N . 



PR2 

0355 

CC27 


N1N= IN-tNl 



PR2 

0356 

CCZS 


N1N1 = IN1+M 



PR2 

C357 

CC29 

A- 

T=A( N1N1>-A<NN) 



PR2 

0358 

0030 


U = T-!<T 



PR2 

0359 

CC31 


V=A.0«A(MN>WA(NN1) 



PR2 

C36C 

CC22 


IF(0A3S( VJ-U«£7) 11,11,5 



PR2 

03 61 

0033 

5 

T=U+V 



PR2 

0362 

0G3A 


IF (DABS (T J-CfU>Cl(U,DABS(V> )>l'E6) 

6, 6, 7 

PR2 

03o3 

CCi'T 

6 

T = C..0 



PR 2 

0364 

0036 

7 

U=(A(.'llNl)+A(NN) ) /2. 0 . 



PR2 

0365 

003 7 


V=CSORT(CABS(T> )/2.C 



PR2 

0366 

C-C35 


IF(T) 15,8,8 



PR2 

0367 

CC29 

8 

IF(U) 10,9,9 



PR2 

0368 

OCAO 

9 

RR(Ni)=U+V 



PR 2 

0359 

CCAl 


•RR (N)=U-V 



PR2 

C37C 

cca; 


GO TO 14 



PR2 

0371 

00 A 3 

1C 

RR(N1 ) = L-V 



PR2 

0372 

OGAA 


RR(N)=U-^V 



PR2 

0373 

CCA3 


GC TO 14 



PR2 

0374 

CCA6 

1 1 

IFITJ 13,12,12 



PR2 

0375 

CCA7 

12 

RR (NU=A( I^INI) 



PR2 

0376 

CCA8 


RR(N)=i(NM 



PR2 

0377 

CCA9 


GO TO 14 



PR2 

0376 

00 50 

13 

RR(N1 )=A(NN) 



PR2 

0379 

cos: 


RR(N)=AtNlNl) 



PR2 

03 50 

CC52 

14 

RI (N)=C.0 



PR2 

0381 

CC53 


RI (N1)=C. C 



PR2 

0332 

0054 


GC TO 16 • 



PR2 

0333 

CC55 

15 

RR (N1 )=U 



PR2 

0334 

CC56 


RR(N) =U 



PR 2 

0385 

0057 


RI(N1» = V 



PR2 

0336 

0053 


RI(N)=-V 



PR2 

0337 

CC59 

16 

IF(N2) 68,68,17 



PR2 

0338 


C 

' 



PR2 

0389 


C 

TESTS OF CONVERGENCE 



PR2 

0390 


C 




PR 2 
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0C6C 

17 

N1N2=M1N1-IA 

PR2 

0392 

V ^ Jk 


RPCO=RR(N’I)’^RR{M )+RI (N1 )*RI{Ni) 

PR2 

0393 

CC62 


£PS=E 1C*DS0RT( RPCD) 

, PR2 

0394 

0063 


IF(0A3S(A(NIV2) )-cPS) 68,68rl8 

PR2 

0395 

CC6A 

18, 

IFIDASS (A<NM ) )-E10’i'DABS t A (NN ) ) ) 69, 69,19 

PR2 

0396 

CC65 

19 

IF<DAf3StPAM-AUlA?l 1 -DASS (A(MN2) )*E6) 67,67,20 

PR2 

0357 

0066 

20 

IF(DA8S(PAN-A( NM) )-0ABS! A(AM) ) ’!'E6> 67,67,21 

PR2 

0398 

0067 

21 

IFdT-’^AXIT » 22,67,67 

PR2 

0399 


C 


PR2 

C4CC 


C 

COMPUTE THE SHIFT 

PR2 

0401 


c 


PR2 

OA32 

CC&8 

22 

J=1 

PR2 

04C3 

cc=? 


DC 24 1=1,2 

PR2 

0404 

0070 


K=\P- I 

PP2 

0405 

0071 


IF(DAaS(RRCK)-PRR(I) )+DABS(RI(K)-PRH n )-0ELTA4 ( DABS( RR( K) )+DABS(RPR2 

04 J6 



IKK) n 1 23,24,24 

PR2 

C407 

0072 

23 

J=J<-I 

PR2 

0408 

CG73. 

24 

CCNTINUS 

PR2 

04 u9 

CC7A 


GC TO (25,26,26,27), J 

PR2 

C41C 

CC75 

25 

R=C, 0 

PR2 

C411 

007-j. 


S = 0.0 

PR2 

0412 

CG77 • 


GC TO 28 

PR2 

0413 

CC72 

26 

J=:H2-J 

PR2 

0414 

0070 


R=RR ( J )v:RR{ J) 

PR2 

0415 

CC30 


S=RR(J)+RR(J) 

PR2 

0416 

cce: 


GC TO 23 

PR2 

0417 

CC52 

27 

R=R^(,\)<c.:iRfM)-RI (N)<,Rl(r<l) 

PR2 

0418 

0 C63 


S=RR(M )+RR (Nl) 

PR2 

0419 


C 


PR2 

C42C 


C 

SAVE THE LAST TWC SUBDIAGCNAL TERMS AND THe'rOOTS OF THE 

PR2 

0421 


c 

SU3MATR IX BEFORE ITERATION 

PR2 

0422 


c 


PR2 

042 3 

CCfcA 

28 

PAN=A(NM) 

PR2 

C424 

OC55 


PAN1=A(MN2> 

PR2 

0425 

0036 


DO 29 1=1,2 

PR2 

0426 

0C87 


K=NP-I 

PR2 

0427 

CC2S 


PRR( n=RR(K> 

PR 2 

0428 

0037 

29 

PRI{ I )=ftI(K) 

PR2 

0429 


C 


PR2 

0430 


C 

SEARCH FCR A PARTITION OF THE MATRIX, DEFINED BY P AND Q 

■ PR2 

04 31 


c 


•PR2 

0432 

CC90 • 


P = N2 

. PR2 

0433 

C CS 1 


IF(N-3) 35,35,30 

PR2 

0434 

CC9Z 

30 

IPI=i'JlN2 

PR2 

C435 

007 J 


DO 34 J=2,N2 

PR2 

043*;. 

0C7* 


IPI=IPI-IA-1 

PR2 

0437 

CC9 = 


IF(OAaS(A(IPn )-EPS) 35,35,31 

PR2 

0438 

0Q9& 

31 

IPIP=IPI+IA 

PR 2 

0439 
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ATEIG 


DATE = 78178 


1^/44/45 


page 0004 


CCS7 


IPIP2=IPI P+I A 

PR2 

044C 

OC93 


0=A( IPI? ) >:=( A( I PIP)-S)+A(I PI P2) «A( IPIP+1) + R 

PR2 

044.1 

CC97 


IF(D) 32,33,32 

PR2 

0442 

ClCO 

32 

I F (OABSI A UFI )*A(IPIP+1 ) l>i‘tCA8S ( At IPIP )+A(IPIP2 + i) 

-$)+DASStAlIPIP2PR2 

0443 



1+2)) )-DABS(D)->EPSl 36,36,33 

PR2 

0444 

ClOl 

33 

P=N1-J 

' PR2 

0445 

c u: 

34 

CCKTIiN'US 

PR2 

0446 

c :c3 

35 

0 = P 

PR 2 

0447 

0L04 


GO TO 39 

PR2 

0448 

0x05 

36 

Pi=P-l 

PR2 

0449 

OK 6 


a=pi 

PR2 

045C 

0107 


IF(Pl-l) 39,39,37 

PR2 

0451 

0x08 

37 

DO 38 1=2, PI 

PR2 

04 52 

G1C9 


I PI = I PI-I A-1 

PR2 

0453 

Clio 


IFtDABStAtrPI) )-EPS) 39,39,38 

PR2 

04 54 

OIU 

38 

Q=Q-1 

PR2 

0455 




PR2 

04 56 


C 

OR 0CU3LE I.TERATIQN 

PR 2 

0457 


c 


PR2 

0458 

0112* 

' 39 

II=(P-1)*IA+P 

PR2 

0459 

U113 


DC 66 I=P,N1 

PR2 

C46G 

0114 


I I 1=II-IA 

PR2 

0461 

0115- 


I i?=n + iA . 

PR2 

0462 

0116 


IFtl-P) 41,40,41 

PR2 

0463 

0117 

4C 

I?I=I I + i 

PR2 

04 64 '* 

0118 


IPIP=IIP + 1 

PR2 

0465 


C 


PR2 

04o 6 


C 

INITIALIZATION OF THE TRANSFORMATION 

PR2 

C467 


c 


PR2 

0468 

0119 


Gx=Al II )»( At 1 1 )-S)+A( IIPj^At I PU+R 

PR2 

0469 

0x20 


G2=At IP I) vt A t IP IP )+At II )-S) 

PR2 

C479 

0121 


G3=At IPI )vA t lPlP+1 ) 

PR2 

0471 

0122 


At IPI + 1)=C. C 

PR2 

0472 

0x23 


GO TO 44 

PR2 

0473 

0124 

41 

Gl=AtIIl) 

PR2 

0474 

0125 


G2=At in+1) 

PR 2 

0475 

0126 


Irtl-''I2) 42,42,43 

PR2 

0476 

6127 

42 

G3=Atni+2) 

PR2 

0477 

0128 


GO TO 44 

PR 2 

0478 

0 129 

43 

G3=0. C 

PR2 

0479 

013o 

44 

CAP=OSORT tGl>?Gl + G2>l'G2+G3>.''G3> 

PR2 

0480 

0131 


IFtCAP) 45,48,45 . . 

PR2 

0431 

0132 

45 

IFlGl) 46,47,4-' 

PR2 

0482 

0133 

46 

CAP=-CAP 

PR2 

0483 

0134 

47 

T=Gl+C4P 

■ PR2 

0434 

0125 


PSI1=G2/T 

PR2 

0485 

0x36 


PSI2=G3/T 

PR2 

0486 

0137 


ALPHA=2 .0/ t 1 .0+PS I 1»PS I 1+PSI 2«PSI 2) 

PR2 

0487 
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0135 

* 

GO TO 49 

C139 

48 

ALPHA=2,.0<. 

C 34C 


PSI1=C.0’ 

0141 


PSI2=C.C . 

0142 

49 

IFll-01 5C,53„5C 

0143 

. 5j 

IFII-P) 51,52,51 

0144 

51 

AUI1)=-CAP 

0>5 


GO rc 53 

0146 

52 

r 

Allli )=-i(IIi) 


U 

c . 

r 

ROW OPERATION • 

0147 

L 

53 

. IJ = II 

0146 


DO 57 J = t ,N 

0149 


T=psi i*a: IJ + 1) 

0150 


IF (I-Nl ) 54,55,55 

c:=i 

54 

IP2J=IJ+2 

0152 


T=T+PST2>^A(IP2J) 

0 1 •: V,. 

55 

£T AsALPHR’^iT + At IJ) ) 

/ • 


A ( I J 1 = A !•: .l-;‘ T. ET A 

r ’.’r 


A(IJtl)=,L<: jd) -PSIL4ETA 

^ • .a 

4. • • 


IF( I-Nl)', 36/57, 57 

OUT, 

• 56 

A(IP2 J) = A( IF2J )-PS I2>:«ETA 

015 3' 

■ 57 

U = I J+IA 


C 

C 

‘COLUMN OPERATION 

0159 

c 

IF(I-Nl) 59,53,58 

0160 

53 

K=N 

Old 


GC TO 60 

Cl 62 

59 

K=I+2 

016 3‘ 

60 

IP=IIP-I 

0164 


CO 64 J=0,K 

0165 


JIP=IP+J 

C166 


JI=JIP-IA 

0167 


T=PSI1'’'A( JIP) 

0163 


Irl'l-Nl) 61,62,62 

Cits 

61 

J1P2=J1P+1A 

0179 


T=T+PSI2»A(JIP2) 

0171 

62 

ETA= ALPHA* (T+A( J I J J 

C172 . 


A(JI)=A(JI)-ETA 

0173 


Af JIP)=A( JIP)-ETA«PSI1 

0 l74 


IF( I-'ll) 63,64, 64 

CL73 

63 

A { JI P2) = A UIP2)-ETA*'PSI 2 

o'; ^ fc 

64 

CONTI 'NUE 

0177 


IF(I-N?) 65,66,66 

0173 

65 

JI=II+3 

0179 


JIP=JI+IA 
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PR2 0488 
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PR2 048C 
PR2 0491 
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PR2 G493 
PR2. 0494 
PR2 0495 
PR2 0496 
PR2 0497 
PR2 049P 
PR2 0499, 

PR2 C5CC 
PR2 ObOl 
PR2 0502 
PR2 0503 
PR2 C5C4 
PR2 0505 
PP.2 05 06 
PR 2 05C7 

P!R2' 'd'509 
PR2 C51C 
PR2 0512 
PR2 0512 
PR2 0513 
PR2 0514 
PR2 0515 
PR2 0516 
PR2 0517 
P.R2 0515 
PR2 0519 
PR2 0520 
PR2 0521 
PR2 0522 
PR2 0523 
PR2 0524- 
, PR2 0525 
PR2 0526 
PR2 0527 
PR2 0526 , 

PR2 0529 
PR2 0530 
PR2 0531 
PR2 0532 
PR2 0532 
PR2 0534 
PR2 0535 
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CISC 

• 

JIP2=JIP+IA 

PR2 

0536 

0 1 e 


tTA=AL'’HAvPSI2*AUlP2) 

PR 2 

0537 

0182 


A(JI)=-ETA 

PR2 

0538 

0.63 


A(JIP )=-=TA=i'?SIl 

PP.2 

0539 

016A 


A( JIP2)=AlJIP2)-cTA“''PSI2 

PR2 

C59C 

0185 

66 

I I =11 P+1 

PR2 

054! 

0186 


IT=IT+1 

PR2 

0542 

0187 


GC TO A 

PR2 

0543 


C 


PR2 

05AA 


C 

ENO OF ITERATION 

PR2 

0545 


C 

» 

PR2 

0546 

oiee 

67 

IF{0A8S(AINM) )-0ABS(A(NlN2) ) ) 69,68,68 

PR2 

0547 


C 


?R2 

0546 


C 

TWO EIGENVALUES HAVE BEEN FOUND 

PR2 

0549 


C 


PR2 

0550 

0185 

68 

IAiNA(,\)=0 

PRZ 

0551 

C15C 


IANA(N1)=2 

PR2 

0552 

01'); 


■ N=N2 

PR2 

0553 

01^2 


IF(N2> 71,71,1 

PR2 

0?54 

• 

C 


PR2 

0555 


C 

ONE EIGENVALUE HAS BEEN FOUND 

PR2 

0556 


C 


PR2 

0557 

0153 • 

69 

RR(NI=A(NM 

PR 2 

0553 

0154' 


RI {.\')=0. 0- 

PR 2 

0 5 59 

0195 


IANA(N )=1 

PR2. 

0560 

0196 


IFINl) 71,71,70 

PR2 

0561 

0197 

7C 

N'!=Nl 

PR2 

0562 

0198 


BO TO 1 

PR2 

0563 

0199 

71 

RETURN 

PR2 

05o4 

02CC 


END 

PR2 

0565 


98 . 



FOSTS^iN IV Gi RELEASE 2,0 


SACK 


DATE = 78178 


1A/A4/45 


page OOOi 


CCC 




CCC2 
CCC3 
0C04 
CCi,5 
CClc 
0C07 
GC J 8 

ccc<? 

OOiO- 

OCil 

c:i2 

CCi3 

0C14 

CC15 

c.: u 

0017 
C0I8 
CC i A 
002J 
0021 
CO 2 2 
CC23 

CCi'r 

c;z 5 

C126 

0027 

0028 
CC24 
00 30 

0031 

0032 
’CC33 

0034 


C 

e 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


I 


2 

3 




5 


6 


7 


subroutine E-ACK (NAHEfNCHAR) • PR2 0566 

PR2 0567 

<><><> OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPR2 0568 

PR2 0569 

CHANCE NAPE TO NCHAR IN MECHANICAL FORMAT . • PR2 057C 

CHANGES PR2 0571 

6P=T? DXXX=ZXXK AND CHANGE XXX PR2 0572 

I=F EXCEPT ICW=TJ PR2 0573 


V=V EXCEPT VCW=KJ 
ELEMENTS 

CU=M EU=U JU=R 

C-A=J EW=G JH-C 


PR2 0574 
PR2 0575 

LU=K RU=D PR2 0576 

PR2 0577 
PR2 057£ 

<><><><><><><><><><><><>< ^<><><><><><><><><><><><><><><><><><><><>PR2 0579 


DIPENSICN NSL(5), NELFIS), NST(2), NSTF(4). NCHAR(6)» NAME(6) 
EQUIVALENCE ( NO ,NELF ( 7) ) , ( NC ,NELF (5 ) ) , (NT,NSTF(3)) 
EQUIVALENCE (NU,NSLFC 2) 1 , (NC,NELF(en, lNVI,NSTF(4n 
DATA N5L/ lKC,U-e; J.HJ, 1HL» IHR/ 

DATA NELF/lHM,lHU,lHR,lHKaHD,lHJ,lHG,lHC/ 

DATA NST/IHI,1HV/ * ■ 

DATA NSTF/11-F,1HV,1HT,1HW/ 

DATA i\2/lhZ/,NP/lHP/,NBL/lH / 

KST=1 

KN=0 

DC 1 1=1,6 

IFtNAPSlU.NE.'NBU KN=KN+1 
NCHAR I n = NAHE{ I ) 

IFUN.LE.l) RETURN 

lF{MCHNfttl>.NS.NG.0R.NCHARl2).Ne.NP) GC TO 3 

NCHAR ( 1)=NT 

RETURN 

DO 4 KTNS=1 ,4 

KL0C=KN-KINS+1 

IF IKLUC.LS. 1 ) RETURN 

IF (NCHARIKLCC) .EC.NU.OR.NCHARtKLOO.EQ.NW J GO TO 5 

CONTINUE 

RETURN 


PR2 05S0 
PR2 0581 
PR2 0582 
PR2 0533 
PR2 C584 
PR2 0SS5 
PR2 0586 
PR2 0587 
PR2 0588 
PR2 0585 
PR2 0590 
PR 2 C5 91 
PR2 C592 
PR2 C593 
PR2 C554 
PR2 0595 
PR2 05 J6 
PR2 0597 
PR2 0598 
PR2 0599 
PR2 06‘t0 
PR2 C6C1 
PR2 0602 
PR2 06J3 


IFINCFARU I.NE.NC) GO TO 6 

NCHARIl) =NZ 

KST=2 

CG 7 Kl=l,2 

IF{NCH4R(KST».N£.NST(K1)) GG to 7 
NCriAR(KST)=NSTr(Ki) 

KST=KST+1 

IFINCHARIKSTI .EC.NC.AND.NCHARCKLQCl. 

GO TO 8 

CONTINUE 


PR2 C6C4 
PR2 0605 
PR2 0606 
PR2 0607 
PR2 06C8 
PR2 0609 
PR2 0610 

,NW) NCHAR(KST“1) = NSTF(KH-2)PR2 0611 

PR2 0612 
PR2 06 i 3 
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0035 

a 

K£L = 5 

P.R2 

06X4 

0036 


K5t.F=:0 

PR2 

0615 

0C37 


IF(NCHAR(KLOC).EC.NU) GO TO 9 

PR2 

0616 

CC38 


KEL=3 

PR2 

0617 

003V 


KELF=5 

PR2 

0618 

CCA*: 

9 

CO 10 K1=1,KEL 

PR2 

0619 



IF(N:hARIKST).AE.,N6UK11 ) go to 10 

PR2 

062C 

00A2 


'lCHARIKST)=NELFtKl+K£LF-l 

PR 2 

0621 

0043 


' GO TO 11 

PP.2 

0622 

CCAA 

10 

CCNTINUe 

PRZ 

0623 

CCA5 


RiTUR.V 

PR2 

0624 

CCA 6 

11 

IFIKLCr .EO.KN) GO TO 13 

PR2 

0625 

CC47 


K£LF=KN-KLCC 

PR2 

0626 

CCA8 


■ DO 12 Kl = i ,KSLF 

PR 2 

0627 

0049. 

12 

MCHARIKLOC+Kl-l>=NCHARfKLCC + J<l) 

PR2 

062S 

C C5 0 

13 

NChAR (KN J = N6L 

PR2 

0629 

c<;3: 


GC TO 2 

PR2 

C63C 

005.2 


END 

PR2 

0631 
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PAGE CCOi 


0001 


CCC2 

CCC3 


CCCA 

0005 

0CC5 

CCC7 

0003 

COCO 

ceil 


ecu 


CC12 

CC13 

CCIA 


0015 
C C 1 
0017 

ccie 

ooiy 

CC2o 

CC21 

0022 

CC23 

002A 

0025 

CC26 

CC27 


C 

C 

C 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


1 

c 

c 


c 

c 


2 

3 
C 

4 


5 


SUBROUTINE EIGAN (ERjEI) PR2 0632 

PR2 0633 

<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< «»>>>>>>»>)>»»>>»>> »»»»p^ 2 0634 

PR2 0635 

SUPERVISE COKPLTATIDN OF EIGENVALUES PR2 0636 

FROM SCEPTRE ROUTINE SIMTR PR2 0637 

PR2 0638 

<<<<<<<<<<<<<<<<<<<<<<<<<<< <<<««<»»»»»»»»»»»»»»>»P22 06 39 

PR2 0640 

IMPLICIT REAL*ft ( A-H,0-Z) PR2 0641 

COMMON /CNTRLS/ TJ NS , SI ( 2 1 ) , XNCPRQ,S2 , X NHfcAC , S3 16 ) , X RUNNO, S4 ( 11 ) , XP.R 2 0642 

IWRTSF., S5( 46) , XNEY, XN5 7DT, XNBC , XNBP. , XNB L , XNLC ,XNLR,XNLL. XNJX ,XNJ8DTPR2 0643 
2 »XNJ3 ,XNJ9,S6 ( 2 ) , XNN3 3r XN!^36, XNM66, XNDP , XNOPO ,XNElTS , S7{ 50) PR2 0644 

CGNMON /CLRNTS/ C(l) ' PR2 064= 

COMMON /OEREVS/ 0(1) PR2 0646 

COMMON /N4M5S/ XNtl) PR2 0647 

CCMMCN /CLTFUT/ GUTP(l) PR2 0648 

COMMON /TAPES/ NPRT ,NDTA , NPG N PR2 0649 

COMMON /VCLTGS/ V(l) ‘ 

COMMON /PkCSE/ A(50,50), 0(5O,5O).L3PtL3S,L4P,L4S,NSVf NEJfNEL»NlND(PR2 
13C0) rIPERM(50) 

CIMENSICN ER(50), EK50), XA(5C,50), SSV(50) 

BEGIN CODE 

SET UP INITIAL CONSTANT OFFSETS 
CALL SIMTR (KOLM.C) 

DO 1 1=1, NSV 
SSV( I )=C< I ) 

LOOP THRU STATS VARIABLES' 

THRU TYPE 3*S FIRST 
IFtL3S.GT.L3P) GO TO 4 
OC 3 I=L3S,L3P 
C(I) = i.COC 

COMPUTE DERIVATIVES , 

CALL SIMTR (KDUM,0) 

FORM AN A MATRIX ROW 
DO 2 J=1,NSV 

A( J, I-L3S-H ) = 0( J )-SSV (J ) 

0( D^C.CDO 

THRU TYPE 4 STATE VARIABLES 
IF (LAS. GT. LAP) GO TO 7 
DC 6 I=L4S,L4P 
V( l') =1. CDC 
CALL SIMTR (KDUM,C) 

DC 5 J=1,NSV 

A( J,I+L3P-L3S-L4S+2) =0( J)-SSV( J) 


ii 

n} 2 , 


PR2 

;pR2 

0650 

C651 


PR 2 

0632 

©nJ 

PR2 

0653 

Cj ^ 

PR2 

0654 


PR2 

0655 

n si 

PR2 

0656 

3 ^ 

PR2 

0657 

PR 2 

0658 


PR2 

0659 


PR2 

0660 


PR2 

0661 


PR2 

0662 


PR2 

066 3 


PP.2 

0664 


PR2 

066 = 


PR2 

0666 


PR2 

06a7 


PR 2 

0668 


PR2 

0669 


PR2 

0679 


PR2 

0671 


OR2 

0672 


PR2 

0673 


PR2 

06 74 


PR2 

0675 


PR2 

0676 


PR2 

0677 


PR2' 

0678 


PR2 

0679 
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CC28 

6 

V(I) =c.cco 

PR 2 

0680 


C 

FORM B MATRIX IF iNECcSSARY 

PR2 

0631 

OCiO 

7 

IFINEJ.eO.O) GO TO 10 

PR2 

0682 

CC3C 


11=1 

PR2 

0683 

0C31 


DO 9 1=1, ^EL 

PR 2 

0684 

0032 


IFIN IMD( I) .NS. 1) GO TO 9 

PR2 

0685 

0C33 


XN ( I )=1 .OEO 

PR2 

0o3 6 

CC34 


CALL SIMTR (KCUM.O) 

PR2 

0687 

00 34 


DO 8 J=l,iNSV 

PR2 

0688 

0C36 

8 

e tJ, I I)=C(J )-SSVU ) 

PR2 

06 8 9 

CC27 


11=11+1 

PR2 

C69C 

CC38 


XN(I J =C. COO 

PR2 

0651 

0C3V 

9 

CONT IMUt 

PR2 

0692 


C 

FCRM EICEN VALUES 

PR2 

0693 


C 

USE ANOTHER ARRAY (SINCE A IS DESTROYED BY HSBG) 

PR2 

0694 

OOAO 

10 

DO 11 1=1, NSV 

PR2 

0695 

C04i 


CO 11 J=1,NSV 

PR2 

0696 

CCA2 

n 

XA(I , J)=A(1,3) 

PR 2 

0697 

00A3 


call HSBG (NSV,XA,50) 

PR 2 

0698 

CC*44 


CALL ATEIG (NSV, XA, ER , El , IPERF , 50) 

PR2 

0699 

CC45 


RETURN 

PR2 

07CC 

CC4£, 


END 

PR2 

07C1 
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CCCl 


SLRROITINc HSBG (K,A, 

I A) 

PR2 0702 



c • 



PR2 0703 



C 

<<<<<<<<<<<<<< <<<<<<<<<<<<<<< <<<«<>>> >>>>»>>»»>»»»»»»»> 

PR2 07 



c 



PR2 07C5 



c 

RcDUOE A PATRIX TC UPPER HESSEM3ERG FORP 

PR2 0706 



c 

N=ORDER OF A. IA=FIRST DIMENSION OF A 

PR2 0707 



c 

IBM SSP ROUTIN'^ 


PR2 07C8 



c 



PR2 07C9 



c 

<< <<<<<<<<<<<<<<<<<<<<<<<<< ««««»»»»»»»»»>»»»»»»PR 2 0710 



c 



PR2 0711 


ccca 


If-PLICn REAUXsa U-H,0 

-Z) 

PR2 0712 


CCC3 


DIMENSrCN Ad) 


PR2 0713 


0C04 


t=N 


PR2 0714 


CCC5 


NIA=L'»I A 

* . 

PR2 0715 


CCC6 


1IA=NIA-IA 


PR2 0716 



c 



PR2 0717 



c 

L IS THE ROW INDEX 

OF THE £LIM-INATIOM 

PR2 0718 



c 


. 

PR2 C719 


CCC7 

1 

IF { L-3) 19,2»2 


PR2 072C 


0003, 

2 

LIA=LIA-1A 

' 

PR2 0721 


CCC9 


L1=L-I . 


PR2 0722 


CCIC 


L2=Ll-l 


PR2 0722 



c 



PR2 0724 


• 

c 

. SEARCH FOR THE PIVOTAL ELEMENT IN THE LTH R£W 

PR2 0725 



c 



PR2 0726 


CCil 


I SLS = LTA+L 


PR2 0727 


03X2 


IPIV=1SL3“IA 


PR2 0728 


0C13 


PIV=OAPS(A( IPIV) ) 


PR2 0729 


CCI^ 


lFlL-3) 6,6,3 


PR2 073C 


CCI5 

3 

M=IPI V-I A 


PR2 0731 


0016 


DO 5 I = L,.«,IA 


PR2 0732 


CC17 


T = OASS (Ad)) 


PR2 0733 

2 o 

0 C Ic 


IFtT-PtV) 5,5,4 


PR 2 0734 


oon 

4 

IP1V=I 


PRZ 0755 


CC2C 


PIV=T 


PR2 0736 

SS 

ccz: 

5 

CONTINUE 


PRZ 0737 


0022 

6 

IF(PIV) 7, 18, ‘7 


PR2 0738 


0023 

7 

IF(PIV-CA8S(A(ISU8) ) ) 

11,11,8, 

PRZ 0739 

2 hM 


C 



PR2 074C 



c 

INTERCHANGE THE COLUMNS 

PR2 0741 

H Q 


c 



PR2 0742 

S ^ 

CC2A 

8 

M=IPIV-L 


PR2 0743 


OC23 


OG 9 1=1, L 


PR2 0744 


0026 


J=H+I 


PR2 0745 


0027 


T = A(J ) 


PR2 0746 


0‘C28 


K = LIA'+I 


PR 2- 0747 


0C29 


A{ J)=A(K) 


PR2 0742 


0030 

9 

AU) = T 


PR2 0749 
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C • 


PR2 

075G 


C 

IVTERCHANGE THE ROWS 

PR2 

0751 


C 


PR2 

0752 

CC3I 


t'=L2-f^/IA 

PR2 

•0753 

00 32 


DO 10 I=Li,NIA,IA. 

PR2 

0754 

CC33 


T=A(I ) 

PR2 

0755 

CC34 ■ 


J = I-N< 

PR2 

0756 

CC25 


Am=A(J) 

PR2 

0757 

0S3& , 

1C 

A(J>=T 

PR2 

0753 


C 


PR2 

0759 


C 

TcRI'S CF THE ELEMENTARY TRANSFORMATION 

PR2 

076C 


C 


PR2 

0761 

0037 

11 

•DO 12 r=L,LIA,!A 

PR2- 

0762 

CC36 

12 

A(I)=A(I)/A(1SUE) 

PR2 

0763 


C 


PR2 

0764 


C 

RIGHT TRANSFORMATION 

PR2 

0765 


C 


PR2 

076 6 

CC39 


J = -IA 

PR2 

0767 

cca:. 


DO 14 1=1, L2 

PR2 

0768 

00^1 


J=J+IA 

PR2 

0769 

CC42 


LJ=L+J 

PR2 

077C 

CC43 


DC 13 K=l,Ll 

PR2 

0771 

004* • 

♦ 

KJ=K+J 

PR2 

0772 

OC-5 


KL=K+L lA 

PR2 

0773 

COifi 

13 

A(KJ) =A(KJ)-A(U)*A(KL J 

PR2 

0774 

CC47 

14 

CONTI NUE 

PR 2 

0775 


C 


PR2 

0776 


C 

LEFT TRANSFORMATION 

PR2 

0777 


C 


PR 2 

0778 

CC46 


K=-IA 

PR 2 

0779 

0045 


DO 16 1=1, N 

PR2 

0780 

0C5 j 


K=K+IA 

PR2 

0781 

CC5I 


LK=K+L1 

PR2 

0782 

oos; 


S=ACLK) 

PR2 

0733 

OC53 


LJ=L-IA 

PR2 

0784 

CC 54 


DC 15 J=1,L2 

. PR2 

0735 

CC55' 


JK = K4-J 

PR 2 

0736- 

0056 


LJ=LJtIA 

. PR2 

0787 

0057 

15 

S=S + A(LJ )t>AUl<>«1.0D0 

PR2 

0708 

CC3£ 

16 

A(LK) =S 

PR2 

0789 


c 


PR2 

0790 


c 

SET THE LOWER PART OF THE MATRIX TO ZERO 

PR2 

0791 


c 


PR2 

0792 

CCS? 


DO 17 I=L,LIA,IA 

PR2 

0793 

0060 

17 

At n = o.c 

PR2 

0794 

OCtl 

la 

L=L1 

PR2 

0795 

C C 62 


GO TO 1 

PR2 

0796 

0065 

19 

R ETURN 

PR2 

0797 
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hSBG 


CC6A 


END 


DATE = 78178 
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FORTRAN 

CCCl 

CCC2 

CCC3 


IV G1 RELEASE 2.0 INTGR8 DATE = 78178 1A/4A/A5 PAGE 0001 


C 

C 

C 

C 

C 


SUBROUTINE INTGR8 PR2 

<<<<<< <<<<<<<<<<<<<<<«««««««»»»»>»»»»»» »»»»»PR2 

PR2 

OUPMY INTGR8 RCUTINE PR2 

■ PR2 

<<<<<<<<<<<<<<< <<<<«««<<<<<<<«<»»»»»»>»»»»»»»»» PR2 
RETURN PR2 

END PR2 


OTS? 

0800 

06CI 

C8C2 

0803 

D304 

03G5 

08C6 


o 9 

if 

A >t 

r-< 

p ry 
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FORTRAN IV SI RSL5ASE 2.0 


INJVERS 


DATE = 78178 


IA/A4/45 


PAse cool 


CCCl 

0002 


0C03 

CCCA 

0005 


0006 

0CC7 

cccs 

OCO'^ 

CCIO 

ecu 

0012 

CC13 

ecu 

0015 
CO 16 
CC17 

0C1% 

CC15 

0020 

C021 

CC22 

0023 

QCiA 

CC25 

0026 

C027 

CC28 

CC29 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROLTINc INVERS ( N AMMAT » AH AT i MR t NC , URRGR) ^^2 

IMPLICIT REAL’i'lMA-M.O-Z) ,INTEGER<'4( I-N) PR2 

PR2 

PR2 


NAME OF MATRIX TO Be OPERATED ON PR2 

THE MATRIX TO BE SOLVED OR INVERTED PR2 

THE NUPBER OF RCWS IN AMAT PR2 

THF NUMOcR OF COLLHNS IN AMAT PR2 

AI NC = MR + 1 FOR EQUATION SOULUlON PR2 

NC = 2 * AR FCR MATRIX INVERSION PR2 

IF THE MATRIX IS SX AGUUAR* OTHERWISE ZER0PR2 

PR2 

<<<<<<<<<<<<<<<<< <<<<<< «<«<<«<»»>»»>» 

PR2 


INNAMMAT 
AKAT 
NR 
. NC 


lERRDR = ONE 


Re At* A NAMMAT 

COMMON /TAPES/ NCUTTP , I NT A P £ , L IBT P, INOUTP , NED I TP , N S AVTP 
COMMON /CNTRLS/ T I M.E , X STO PT , X MXERT ,X I S , XT I SSS ,X MN ISS t XMX ISS 
1,XMXAI-,XIC,XR=RM0, XH XP A S , XHNR IE , XHXR I S , XHNI SC , XMXl C P , XI CRE 
2P. .XMiXCTP.XXCTMS ,XYDIHS.XIC?AS,XNOPRQf XNTNOXf XNHcADi XNOFEQ) X 
3£RN,XCNTNL,XSTPNC,XPASN0,XRUNNC,XSTPSZ 
DIMEMSTON AMAT!NR,NC) 

CCUIVaLcNCE (NCUTTP.OUTAPE) 

INTE3EP OLTAPE 

IFIXICPAS.LE.XMXICP) GO TO 2 
WRITS tNO'JTTP,!) XMXICP 

FORMAT (34H1THE NE WT CN-RAPHS ON PASS LIMIT OF rF5.0»48H HAS BEEN 
ICESOEO WITHOUT ATTAINING CONVERGENCE! 

I ERB0R=1C 
GC TO 20 
ItRROR=C 
NRX2»2>>NR 
INK=:NR 

IF(MRX2.N£.NCJ INK=NC 
CHECK FOR 1X1 MATRIX 
IFINR.NS.i ) GO TO 3 
AMATUt2> = AMATU»2>/AMAT11U) 

GO TO 2C 
NB2=2wNR 

IF MATRIX INVERSION DESIRED. 

IF(NC.LT.NR’) GC TO 6 
J INC=.MR + 1 
CC 5 UiCR=l,NR 
INCRl3W=INCR + NR 
00 5 INCOLMsJINC ,NC 
IFIINCGLH. SO. INCROW) GO TO 4 
AMATIIN'CR, INCGLM1=0. 

GO TO 5 


PR2 

PR2 

,XHNAIEPR2 
R,XICAEPR2 
RT. XNRPR2 
PR 2 


AND SOLVE IF FOUND 


GENERATE IDENTITY MATRIX 


PR2 
PR2 
PR2 
PR2 
PP.2 
EXPR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR 2 
PR2 
PR2 
PR 2 
PR 2 
PR2 
PR2 
PR 2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 
PR2 


OSOT 

OOOR 

08CG 

081C 

0811 

0812 

0813 

0814 

0815 

0816 
0317 
0813 

0819 

0820 
0821 
0822 

0823 

0824 

0825 
0626 

0827 

0828 
0629 
063C 
083 1 

0832 

0833 

0834 

0835 
C836 
0837 

0836 
0839 
084C 
0841 
0342 
C843 
0844 
uS45 

0846 

0847 
034P. 
0849 
C35C 
0851 
0832 
0853 
0654 
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LVVEftS 


DATE = 78178 


1A/AA/A5 


PAGE 00D2 


CC2C 

4 ‘ 

APAT( INXRt I^COU^')=^ . 


PR2 

0855 


0031 

5 

CONTINUE 


PR2 

0856 


0C32 

6 

CO 19 I=i,NR 


PR2 

0857 


C033 


P=I + 1 


PR2 

0858 


0C3A 


J = I 


PR2 

0859 



C ■ 

CHECK FOR LAST DIAGONAL ELEFENT 


PR2 

0860 


CC35 


IF II. 50. NR) GO TO 9 


PR2 

0861 



c 

FIND THE LARGEST VALUE BELOW DIAGONAL ELEMENT IN THIS 

column 

PR2 

C862 


00 36 


DO 8 N=K,NR 


PR2 

0863 


0037 


I F (DABS (AM AT (J, I ) )-DAB S ( AKA T ( N , I ) ) ) 7r8,8 


PR2 

0864 


CC3S 

7 

J = N 


PR2 

0865 


CC29 

8 

CONTINUE 


PR2 

0866 



C 

SEE IF THE LARGEST ELEKEM IK ROW IS ZERO 


<’R2 

0867 


OCAO 

9 

IF (AKAT U, I) ) 12,10,12 


PR2 

0363 



C 

INDICATES THAT KATRIX IS S’lKGULAR 


PR2 

0869 


OOAl 

1C 

I=RRDR = n 


PR2 

0870 


CCA2 


WRITE (NOUTTPill) NAMMAT 


PR2 

0871 


CC43 

11 

FORf'AT 18H MATRIX ,A4,12H IS SINGULAR) 


PR2 

0372 


CC4A • 


GO TO 2C 


PR2 

0872 



c 

SEE IF ROW INTERCHA'JGE IS NEEDED TO GET LARGEST ELEMENT ON DIAG, 

PR2 

0874 


CCA5 

12 

IF(J-I) 15,15,13 


PR2 

0875 



C 

INTERCHANGE ROWS OF OIAGCNAL AND THOSE TO THE RIGHT OF 

IT 

PR2 

0876 


OOA6 

13 

DO 14 N=I,NC 


PR2 

0877 


OC A 7 


SAVE=AWAT ( I ,N) 


PR2 

03 78 


CCA6 


AMAT( I,N)=AKATf J,N) 


PR2 

0879 


00A9 

lA 

AMATIJ,N )=SAVE 


PR2 

0380 


0C5C 

15 

diag=a«at ( I, I) 


PR2 

OSSl 



C 

CIVICS ALL ELEMENTS TO RIGHT OF DIAG. IN DIAG. ROW BY 

DIAG. ELEM 

.PR2 

0882 


CC51 


DO 16 J=M,NC 


PR2 

0833 


0C52 

16 

AMAT( I,J)=AMAT(I,J)/DIAG 


PR2 

0884 


:c?3 


AKAT( I , I )=0 . 


PR2 

0885 


CC 3A 


DO 18 K=i,NR 


PR2 

0686 


00 55 


IF(AVAT(K,1 ).EQ.O.) GO TO 18 


PR2 

0887 

o9- 

0C55 


I3=I+i 


PR2 

0838 

^ Q 

CC57 


DC 17 J=IP,KC 


PR2 

0889 


C 15 S 

17 

AKAT(K,Jl=AKAT{K,J)-AKAT(K,nKAKAT( I,J) 


PR 2 

089C 


0C59 

U 

CCNTINUE 


PR2 

0891 

Is 

CC4.0 

19 

CCKTIKUE 


PR2 

C892 

CC61 

2C 

CONTI NLE 


PR 2 

0893 

iD H(j 

0062 


RETURN 


PR2 

0894 


CC63 


END 


PR2 

0895 








^ si 
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R5LEAS5 2.0 


CUTP 


DATE = 78173 1A/4V45 PAGE 000 


OOOi 


CCC2 

0003 


C 

C 

c 

c 

c 


subroutine OUTp ■ PR2 

<<<<<<<<<<<<<< <<<<<«««<<<<<<<«<»>>> >»>>>5>»»»>>>?>>»»»PR2 

PR2 

DUMMY nUTP ROUTINE ■ PR2 

PR2 

«<««««<«<<<«««««<<<<<<«>>>»>»»>5>>>»>>>>>»>»»»PR2 
RETURN PR2 

END PR2 


0896 

0897 
0098 
0899 
G900 
09C1 
09C2 
0903 
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FOST'^AN IV G1 RELEASE 2.0 


SAVEPO 


DATE = 78178 


14/A^/A5 


PAGE OOOJ 


OCOl 


CCC2 
000 5 


C 

c 

c 

c 

c 


=>R2 OOO'i 

r<f<<<<r<<<<<r<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>^ 

PR2 0907 

DUMMY SAVEPD ROUTINE _ PR2 0903 

RETURN PR2 0911 

END 
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FO'^TRAN !V Gi i^ELEASE 2.0 


SETUP2 


OATS = 73173 


14 / 44/45 


PAGE OCO 


0001 


SII8I10LTINE SETUP2 { NAC, LEVS ,NSRC ,NTA30 

PR2 

0912 


c ■ 


PR2 

0913 


C 

<<<<<<<<<<<<<<<<<<<<<<<<<<< ««««»»»»»»»»»»>»»»»»P^Z 

0914 


c 


PR 2 

0915 


c 

SETUP ROUTINE 

PR2 

0916 


c 


PR2 

0917 


c 

<<««««««<««««««««««>>»»»»>»>»>»»»»»»»P^2 

0918 


c 


PR2 

0919 

0C02 


IMPLICIT RSAL-fSlA-HfO-Z ) 

PR2 

0920 

CCC3 


• CCMMGN /CM9LS/ T IMS, SI (2 1 ) .XACPRCtS 2» XNHEAG, S3 t 6 W XRUNNOf S4U1) 

, XPR2 

0921 



IWRTS.M, S5f 4 6) ,XMEY,XNc7DT,XN8C ,X.M8R,XNBL,XNLC,XNLR,XNLLf XNJX,XNJ3 0TPR2 

09 2 2 



2f XlNJ0,Xi^lJ97S6( 2), X.'iMSS, XMM36 , XjNM66 . XND P , XNO PD r XKE LT S » S7 ( 50! 

PR2 

0923 

CCC4 


CCMPC-N /TAPES/ XPRT ,,SiOT A, NPGM 

PR 2 

0924 

CCC5 


. DIMcflSICN LcLS{6f3C0) 

PR 2 

092 5 

.0006 


OIMENSIOU NSIHTRl5!f NSRCDRI6), NTABLE(6), NNDP0(3) 

PR2 

0926 

CCC7 


CIHE.--ISICN KCUTPTI6), NXMUTLtS). .M.MAMES(5)i NPRMTR(5I 

PR2 

0927 

CCC3 


DIMSNSICN MCNnx(5), ,NDcRIV(6!) NRATISOlt CTROLSU60) 

PR2 

0928 

0CC9 


DIMENSION NLMSIlC), LOG SI 8) 

PR2 

0929 

CCiO 


data ns IHTR/IHS, IHI, IHM, IHT, IHR/ > 

PR2 

09 30 

ecu. 


DATA NSRCCR/IHS jIFRtIHC. IPC, IHR, IHV/ 

PR 2 

C93! 

CC12 


DATA NTABLE/iHljlHA ,1HB t IHL ,1HE ,1 HS/ 

PR2 

0932 

0.-13 


data XMCPC/IH, rlhOilHP/ 

PR2 

0933 

CCiA . 


DATA NCUrPT/lHCilhU,lHT, IhPtlPT, IHS/ 

PR2 

0934 

CC15' 


DATA NXMUTL/lHX,lHM,iHUilHT ,IHL/ 

PR2 

0935 

0016 


•DATA NMAM£S/lHN,iHA,lHn,lHE,lHS/ 

PR2 

09 3 6 

CC17 


DATA APRyTR/lnF, ih.R, IHM, IhT, IhR, IHS/ 

PR2 

C937 

CCIP 


DATA MCi\BX/ 1HT >iHC,iHN,iH ,iHX/ 

PR 2 

0938 

0010 


DATA JDERIV/iHDjiHS ,1HR,1HI jiHVflHS/' 

PR2 

0939 

GC20 


DATA LnCS/96,97,93,93,94, 95,92,100/ 

PR2 

0940 

ccc: 


DATA M-MS/lHO,i (-1 , 1 H2 , 1 H3 , 1 (-4 , 1H5 , IH6, 1H7, 1H8, 1H9/ 

PR2 

0941 

0022 


DATA NAST/1H<=/ 

PR2 

0942 

0023 


ECUIVAL-2.NCE (T IM = ,CTROL S 1 1 ) ) 

PR2 

0943 

CC 2 ^ 


ECUrVJLENCS |iNCOM,NNDPD(I ) ) 

PR 2 

0944 

n ^ - c 
4 - 


SOLI VALENCE ( NR LK ,NI 0N8X (4)') 

PR 2 

0945 

00 26 


NTAB*0 

PR2 

0946 

C027 


NSRC=J 

PR2 

0947 

00 23 

L 

READ (.\PGM,2,£ND=33 ) NRAY 

PR2 

0948 

0029 

2 

FORMAT UCAl) 

PR2 

0949 

0050 


IF (L0CATE(NRAY,NSIMTR,9, 5).EQ.l) GO TO 1 

PR2 

09 50 

CC31 


READ (NPGM.3! NRAY 

PR2 

C95I 

003? 

3 

FORMAT (////////,80Ai) 

PR2 

0752 

0033 


NAC=70 

0R2 

0953 

C':34 


IF{NRAY(6).EC.NASn NAC=160 

PR2 

0954 

CC35 


READ INDIA, 4rEAD=35) (CTRCLS ( I ) , 1=1 ,NAC ) 

PR2 

0955 

0036 

4 

FORMAT UP 50 14. 7) 

PR2 

0956 

CC3 7 


I F (NAC.EG.160 ) GO TO 6 

PR2 

0957 

CC3 3 


00 5 1=71,160 

PR 2 

0958 

0039 

5 

CTROLS( I)=C.COC 

PR2 

0959 
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FCftTS^^ IV G1 RELEASE 2.0 SETUPS 


0C4C 

‘6 

READ (NPGP.a) NRAY 

0041 

• 7 

IF(LDC4TE(NRAY,NSRC0R ,15,6) .EC. 0) 

0C42 


I F (LOG AT 5 (NRAY, NT ABLE, 15, 6) .CQ.O) 

CC43 


IF(NAC.EC.160) CC TC 8 

0044 


IF(LDCATE<NRAY,NXMUTL,15 ,5) . EC.O) 

0045 


IF(L0CAT6(NRAY,NNAI1£S, 15, 51 .£0.0) 

CC46 


I F : LOCATE {NR AY, NPRHTR, 15,6 ) .60.0 ) 

CC47 


IF<LOGAT£ <NRAY,.\D£RI V,15,6) . EC.O) 

0043 


IF{LDCATe(NRAY,NIC.NBX,13,5) .EQ.O) 

CC49 

8 

IF (LOCATE (NRAY.NCUTPT, 15,6 ) .£0.0 ) 

CCSC 


GO TO 5 

0051 

9 

K = 22 

0052 

10 

Nl-OL=NrNT (NRAY.K ) 

CC53 


NSRC=NSRC+NHOL 

0054 


IF(NRAY( 72) .EQ. M3LK) GO TC 6 

0055 


READ ('‘lPGf‘1,2) NRAY 

CC56 


IF(NRAY(6 l.NE.NAST) GO TO 7 

0057 


K = 7 

0053 


GO TO 1C 

C'C59 

U 

NT AB= NUMB (NRAY, 25,27,K ) 

CCtO 


GC TO 6 ' 

CCcI 

12 

K = 21 

0062 

13 

VH0L=NCNT(NRAY,K) 

0063 


XNY33=XNP33+NHCL 

CC64 


IF(\RAY(72) .EQ. NBLK) GO TC 6 

0065 


READ (NPGK,2) NRAY 

0066 


■iF(NRAY(6).N£.NAST) GO TO 7 

CC67 


K = 7 

00o3 


GO TO 13 

006 9 

14 

NELTS=0 

CC70 


1=21 

CC7i 


N = 1 

0072 

15 

DO 16 J=l,6 

CC73 

16 

LELS ( J,NELTS+1 )=NBLK 

CC74 

17 

IF(NRAY(I l.EC.NCCM) GO TO 18 

0075 


LELS(N,NELTS+1)=NRAY(I ) 

007o 


IF(NRAY(I).£0.N8LK) GO TO 19 

CC77 


N = N+1 

C07S 


1 = 1 + 1 

0079 


IF( I.LE.72) GO TO 17 

CC80 


READ (NPGM,2) NRAY 

ccei 


IF(NRAY(6).N£.NAST) GO TO 20 

0032 


1 = 7 

0083 


GC TO 17 

CC64 

18 

NSLTS=N£LTS+i 

0035 


XNcLTS=N£LTS 

0086 


N=1 

CC£7 


1 = 1+1 


DATE =, 7S178 


I4/4A/45 


PAGE 0C02 


GO TO 9 
GO TO U 

GO TO IZ 
GO TO 14 
GO TO 21 
GO TO 23 
GO TO 25 
GO TO 32 


PR2 096C 
PR2 0961 
PR2 0962 
P5L2 C963 
PR2 0964 
PR2 09o5 
PR2 C966 
PR2 0967 
PR2 0968 
PR2 0969 
PR2 097C 
PR2 0971 
PR2 0972 
P^2 0973 
PR2 0974 
PR2 0975 
PR2 0976 
PR2 0977 
PR2 0976 
PR2 0979 
PR2 C93C 
PR2 0931 
PR2 0982 
PR2 0933 
PR2 0934 
PR2 0935 
PR2 0936 
PR2 0937 
PR2 0938 
PR2 0939 
PR2 099C 
PR2 0991 
PP.2 0992 
PR2 0993 
PR2 0994 
PR2 0995 
PR2 0996 
PR2 0997 
PR2 0996 
PR2 0999 
PR2 ICOO 
PR2 ICCl 
PR2 1002 
PR2 1003 
PR2 1CC4 
PR2 1005 
PR2 1006 
PR2 1007 
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FO‘?T^A.\ IV G1 RELEASE Z.O 


S£TL'P2 


DATE = 78178 


OCR? 

• 

IF( I.Le.72) GO TO 15 

CC£5 


RFAO (NPGP,25 NBAY 

CC9C 


IF (NRAY(6) . ^E. ^AST) GO TO 7 

0C91 


1 = 7 

OC32 


GO TO 15 

CC93 

19 

N£LTS=N'ELTS+l 

0094 


XVELTS=NELTS 

0C95 


GO TO 6 

CCSo 

20 

NELTS=N£LTS+1 

CC97 


X.')£LTS=MELTS 

0C93 


GO TO 7 

CC97 

21 

K=22 

01 CO 

22 

NHCL=AnNTlNRAYtK') 

0101 


XNOP=XNDP+NHOL 

0102 


I F (MRAY ' 72 J -EG.NBLK ) GO TO 6' 

G1C3 


READ (NPGP,2) KRAY 

010 4 


IF :mRAY( 6) . Nc. AAST) GO TO 7 

0105 


K = 7 

C1C6 • 


GC TO 22 

C1C7 

23 

DO 24 I=22t70 

0103 . 


IF(L3CATE(NRAY,^NCPD tl ,3).EQ.C) XNDP0 = XNDPD+1 

C1C9 

24 

CQNT INUE 

0110 


GCT06- 

on: 

25 

IF (NRAY( 18) .EQ.NRAY( 9)) GC TO 6 

C112 

'26 

K= 19 

0113 

IF (NRAY (K ) .EC.NeLK) GO TO 6 

0114 


00 27 I=2t10 

0113 


IFlNRAYfK l.EO.NUKSUn GO TO 29 

0116 

27 

CCNTIllUE 

0117 


WRITE (.NPRT,28) KRAY,K 

0113 

28 

FORMAT { 42H ERROR 4*4 - - INVALID NUMBER IN 

10L=,I3) 

0119 


STOP 3 

0120 

29 

J=I-1 ■ 

0121 


K = K+2 

C122 

30 

IF(K.GT.72) RE4C <NPGM, 2 ) NRAY 

0122 


IF(K.GT,72) K=K-72+6 

0124 


IF(NRAY{K).E0.^'8LK) GO TO 31 

0125 


K=K f2 

0126 


GC TO 30 

0127 

31 

XST=K 

•0'-24 


K :n=k+2 

0129 


CTROLS (LGCS U) 1=NUMe{NRAY,KSTtKEN,K ) 

C1.2C 


K = KEN+5 

0131 


IF(K.LE.72) go TO 26 

0122 


READ (NPGM»2) NRAY 

0123 


IF(NRAY(6),NE. NAST) GO TO 7 

0134 


K=K”72+6 
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PR2 1029 
PR2 1C3C 
PR2 1031 
PR2 1032 
PR2 1C33 
PR2 1C3A 
PR2 1035 
PR2 1036 
. • PR2 1037 

S1KTR,/,80AI,6H CPR2 1033 
PR2 1039 
' PR2 ICAC- 
PR2 1041 
PR2 1042 
PR2 1C43 
PR2 1044 
PR2 1045 
PR2 1046 
PR2 1047 
PR2 1048 
PR2 1U49 
PR2 1C5C 
PR2 1C51 
PR2 1052 
PR2 1C53. 
PR2 1054 
PR2 1055 



SETUP2 
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0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 


32 

33 

34 

35 

36 


GO TO 26 
ft tTURN 

WRITE (NPftT,34) 
format 139H ERROR 
STOP a 

WRITE (NPRTi361 

format 140H *** ERROR *** 


END 


DATE = 73178 


14/44745 


page 0004 


end of file in SIMTR) 

END OF FILE ON OTASAVJ 
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FORTRAN IV Gl" 


RELEASE 2«0 


■ 0001 


CCC2 

0003 


C 

C 

C 

C 

C 


SUBROUTINE TPHEAO . 

<<<<<<<<<<<<<<<<< <<<<<<<<<<<< <<<<<<>>>> >>>>>>>>> >>>»>>>>PR2 

PR 2 


DUMMY TPHSAO ROUTINE' 


PR2 

PR2 


<<««««««<«««««««<<««»> 

RETURN • 

END 


PR2 

PR2 
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1066 

1067 

1068 
1069 
107C 

1071 

1072 
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XLOAO 
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PAGE COOL 


0001 


CCC2 

0003 


C 

C 

c 

c 

c 


SUOROLTINl XLOAD , PR2 

<<<<<<<<<<<<<<< <<<«<««<<<<<<<<«»>>»>>»>>>»»>»>»»>»>»PR2 

PR2 

DUMMY XLOAD ROUTINE PR2 

P R2 

<«««««<«««««««<<<<<««»>>>>»»>>>»>>>»>»»»»»PR2 
RETURN PR2 

END . PR2 
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1076 

1077 
1076 
1C7S 
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locate 


date t 781.78 


14/4A/45 


PAGE OCOl 


OCOl 


OC02 

0003 

00C4 

CCC5 

CCC6 

0007 

ccca 

OC.C-; 

0010 


c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 

•c 


1 


FUNCTION LOCATE { NRA Y , NSTR t K STR r KLcN ) PR2 

PR2 

<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><; ;><><; 2 

PR2 

LOCATE STRING (NSTR) IN ARRAY (NRAY) WITH PR2 

LENGTH (KLSM STARTING AT (KSTR). PR2 

LaCATS-C IF STRING FOUND PR2 

LGCATS=1 IF STRING NOT FOUND PR2 

PR2 

<><><><><><><><><><><><><><><><><><><><><><><><><><><>0<>0<;><>.(;v,p:^2 


DINEWS ION NRAV{80I» WSTRI80) 

L0CAT==1 

CO 1 I=1,KLEN 

KSTRIR=KSTR+I-1 

IFlWRAVIKSTRlNI.WE.NSTRan RETURN 
CONT I\'U£ 

L0CATE=0 

RETURN 

END 


PR2 

PR2 

PR2 

PR2 

PR2 

PR2 

PR2 

PR2 

PR2 

PR2 

PR2 


1081 

1082 

1C88 

1084 

1065 

1C86 

1C37 

1038 

1C8S 

1C9C 

1091 

1092 
1C93 

1094 

1095 
1C96 
1097 
1099 
1C9S 
IICC 
1101 
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0001 
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CCC3 
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0007 
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C 

c 

c 

c 

c 

c 

c 


NCNT 


DATE = 78178 
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PAGE COOi 


FUMCTION NCNT (NRAY.NST) 


PR2 

PR2 


P R 2 
PR2 
PR2 
PR2 


COUNT NtMSER OF COKRAS+1 IN ARRAY NRAY 
STARTING AT NST' 




DIPcNSICK NRAY (60) 

DATA NBLK/IH /,NCG)'/1H,/ 

NCNT=1 

DC 1 I=NST,72 

IF(NRAY(I).EQ. KCCH) NCNT=NCNT+1 

IrtNRAYi I ) .cQ.NBLK) RETURN 

CONTINUE 

RETURN 

END 


PR2 

PR2 

PR2 

PR2 

PR2 

PR2 

PR2 

PR2 

PR2 


1102 

1103 

1104 
11C5 
1106 
1107 
11C8 

1109 

1110 
1111 
1112 

1113 

1114 

1115 

1116 

1117 

1118 
1119 
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CCOl 


0002 

0CC3 

CCCA 

0005 

0CC6 

CCC7 

CCCS 

0007 

CCIO 

ecu 

0C12 
0C13 
CCIA 
CC15 
CC16 
CO 17' 

ccia 

0019 

CC20 

CC21 


C 

C 

C 

C 

C 

C 

C 


1 

2 

3 

A 


FUNCTION NUPB t NRAY , KST » KEN, I ) PR2 112C 

PR2 1121 

<<<<<<<<<<<<<<< <<<<<<<<<< <<<<<<«<<»>> >>>>>>> »»»»>»»»»»> PR2 1122 

PR 2 1123 

DECODE NUMBER AT KST IN NRAY UMTL KEN PR2 112A 

PR2 1125 

<<<<<<<<<<<<<<< <<<<<<«««<<«««»>>» >>»>5>>»>»»»»»»»PR2 ii26 

PR2 1127 

COMMON /TAPES/ NPRT,NOTA ,NPGN PR2 1128 

DIMENSION NRAYISO), NUMSIlO) PR2 1129 

DATA NUPS/iF0,lHTlF2ilH3»lH4'»iK5,lH6,lH7,lM8j'lH9/ PR2 113C 

DATA NBLK/IH / PR2 1131 

NUHG=0 PR2 1132 


DC 4 I=KST,K£N PR2 1133 

IFII,LE.72) GO TC 2 PR2 1134 


READ (MPGMrl) NRAY PR2 1135 

FORMAT (SuAU PR2 1136 

1=7 PR2 1137 

KEN=KEN- 72-t-6 PR2 U38 

continue PR2 1139 

IFINRAY (I l.EQ.NELKJ GO TO 4 PR2 1I4C 

OC 3 J=i*iC PR2 1141 

IFI.NRAYt I )..N£.NUMSIJ) ) GO TO 3 PR2 1142 

NUPB=M!Pe*10+J>l PR2 1143 

CONTINUE ■ PR2 11^4 

continue PR2 1145 

RETURN PR2 1146 

END PS2 1147 
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CCOI 


FUNCriCN XTABLE (TABLE, X) 

PR2 

1148 

CCC2 

C 

IMPLICIT REALM'S (A-H,D-2) , lNTEGER>i‘4 ( I-N) 

PR2 

PR2 

1141 

U5C 


C 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPR2 

U5l 


C • 

FUNCTION TO PERFORM STRAIGHT LINE EXTRAPOLAT I CN OF TABULAR OATA 

PR2 

1152 


C 

PR2 

115 3 


C 


PP2 

il54 


C 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPR2 

1155 


C 


PR2 

1156 

0003 


DIMENSION TABLE! 2) 

PR2 

1157 

CC04 


LCDP=TABLE(1)-2. 

PR2 

1158 


c 

LCCK FOR INTERVAL OR EXACT VALUE OF INDEPENDENT VARIABLE 

PR2 

1159 

CCC5 


DO 4 INDX=3,LOOP,2 

PR2 

1160 


c 

- TEST LOWcR COORCINATE 

PR2 

1161 

CCCft 


IF(TAeLElINCX)-X) 1,2,5 

PR 2 

1162 


c 

TEST LPPES CGCRDINATE 

PR2 

1163 

0007 

1 

IF(TA8LE( INOX+?)-X) 4,3,5 

PR2 

1164 


c 

EXACT MATCH TO LOW COORDINATE - EQUAT AND RETURN 

PR2 

■1165 

occe 

2 

Y=TABLE(INDX+1 j 

PR2 

1166 

ooov 

. 

GO TO 6 • ■ 

PR2 

1167 


c 

EXACT MATCH TO UPPER COORDINATE - EQUAT AND- RETURN 

PR2 

1168 

cci: 

3 

Y = TABLE(UCX+3 ) 

PR2 

1169 

ooi: 


GO TO 6 

PR2 

1170 

0012 

4 

CCNTINU6 

PR2 

1171 

CC13 


INDX=LCCP-1 • 

PR 2 

1172 


C 

PtFORM STRAIGHT LINE EXTRAPOLATION 

PR2 

1173 


C 

CALCULATE SLOPE 

PR2 

117A 

CCl^ 

5 

S LOPE=(T ABLE (I NCX43 t -TABLE! INCX-H ))/( TABLE ! INDX-42 )-TABL£ ( INOXI ) 

PR2 

1175 


C 

CALCULATE Y INTERCEPT 

PR2 

1176 

0015 


YINTR = TABL6! 1NDX+ 1) -SL0PE*TA3 LE ! I NDX) 

PR2 

1177 


C 

calculate depencent variable 

PR2 

1178 

CC16 


Y=SLCF£*XtYINTR 

PR2 

1179 


C 

OUTPUT PERTINANT VARIABLES 

PR2 

1180 

0017 

6 

XTA6LE=Y' 

PR2 

1181 

CC18 


return- 

PR 2 

1182 

CC15 


end 
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